Stabilizer states and 
local realism 



oo 

o 

O 



by 



l> 



Matthew B. Elliott 






B.A., Mathematics & Physics, Carleton Cohege, 2003 



> 

00 
(N 

l> 

O 
OO 

o 
> 

X 



DISSERTATION 

Submitted in Partial Fulfillment of the 
Requirements for the Degree of 

Doctor of Philosophy 
Physics 



The University of New Mexico 
Albuquerque, New Mexico 



July 2008 



Ill 



Acknowledgments 



There is an exceptional body of researchers at UNM I would like to thank, includ- 
ing my advisor Carl Caves, close collaborator Bryan Eastin, as well as collaborators 
Steve Flammia, Andrew Landahl, and Seth Merkel. Thanks is also deserved by 
others here at UNM, Sergio Boixo, Aniniesh Datta, Ivan Deutsch, and Anil Shaji. 

I also have publications involving researchers not at UNM. My gratitude goes out 
to Jon Barrett, Dan Browne, Akimasa Miyake, Stefano Pironio, and Tony Short. I 
should also mention those not at UNM with whom I have had interesting physics 
discussions, namely Joonwoo Bae, Howard Barnum, Jaewan Kim, Alex Monras, and 
Patrick Rice. 

The quantum circuits in this document were drawn using a wonderful ETgXpackage 
called Q-circuit [39]. 



Stabilizer states and 
local realism 



by 
Matthew B. Elliott 

ABSTRACT OF DISSERTATION 

Submitted in Partial Fulfillment of the 
Requirements for the Degree of 

Doctor of Philosophy 
Physics 

The University of New Mexico 
Albuquerque, New Mexico 

July 2008 



Stabilizer states and 
local realism 

by 

Matthew B. Elliott 

B.A., Mathematics & Physics, Carleton College, 2003 
Ph.D., Physics, University of New Mexico, 2008 



Abstract 

The central focus of this work is to make progress towards understanding entan- 
glement as a resource for computation by examining the quantum correlations that 
can be extracted from stabilizer states. As such, we focus on the stabilizer formalism, 
local realism, and the convergence of the two. 

We contribute to the understanding and manipulation of stabilizer states by pro- 
viding a simple and compact graphical description of them. Using our stabilizer-state 
graphs we demonstrate how Clifford operations and Pauli measurements on stabilizer 
states can be illustrated graphically. Moreover, we demonstrate how our graphical 
description of stabilizer states can actually be extended to encompass all stabilizer 
codes. 

The inability to explain quantum measurement results in local realistic terms 
is the most striking divergence of quantum mechanics from classical physics. We 
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consider measurements on spatially separated quantum systems and learn when the 
measurement correlations are consistent with local realism. If they are, we provide 
an algorithm that generates a local realistic model for them. When the measurement 
correlations are not so describable, our algorithm generates a Bell inequality violated 
by them. Unfortunately, this program becomes computationally infeasible when the 
number of measurements becomes large. For this reason, we also explicitly give a 
large class of Bell inequalities for an arbitrary number of systems and dichotomic 
observables. We also examine local realism through probabilistic theories, which 
include local realistic models, quantum mechanics, and more. 

Finally, we study the stabilizer-state violation of local realism. We begin by exam- 
ining local realistic models for Pauli measurements on stabilizer states which satisfy 
a single highly-motivated property. Beginning with these models, we then examine 
how much stronger quantum correlations are than classical ones. This is done by 
allowing some classical communication. Under this paradigm, we introduce a simple 
communication model that reproduces all global quantum correlations. Studying 
the successes and failures of this simple model uncovers two properties essential for 
any successful communication model. Furthermore, we construct a universal model 
with these properties that reproduces all quantum measurement results for Pauli 
measurements on stabilizer states. 
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Chapter 1 



Preliminaries 



1.1 Introduction and overview 

The field of quantum information has taught us that information-processing tasks are 
more appropriately discussed in the context of the physical apparatus implementing 
them. We cannot separate the physical system from the computational task. More- 
over, our physical reality seems to be most fundamentally quantum mechanical, so it 
is reasonable to approach information theory assuming a physical system governed 
by the laws of quantum mechanics. Moreover, processing information encoded in 
states of quantum systems seems to be advantageous [3, 44, 48, 52, 88], a fact that 
has even changed our notion of computational complexity [14, 96]. However, despite 
these results of quantum information, there still remains the deep question of why 
quantum mechanics makes such feats possible. 

From the perspective of measurement-based quantum computation [79] the full 
power of quantum computation can be consigned to quantum correlations. This 
means studying the power of quantum computation can be done by studying out- 
comes of measurements on an appropriate quantum state. This remarkable discovery 
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says that studying objects such as Bell inequalities [13] is akin to studying the power 
of quantum information itself. Most research in this direction concerns itself with 
the amount by which quantum states can violate the local realistic bound of a Bell 
inequality. However, one interesting related perspective [71], which was first explored 
in Ref. [fil], is to learn how much quantum correlations violate local realism by giv- 
ing the parties access to local realistic correlations and asking how much additional 
information, in the form of classical communication, is needed to reproduce quantum 
correlations. Notice that the original shared randomness is essential, since, without 
it, results diverge considerably [ ]. This is the approach we choose to follow, mean- 
ing we study the power inherent in quantum information by examining local realistic 
models assisted with classical communication. 

We do not study arbitrary measurements on arbitrary quantum states. Instead we 
restrict our attention to Pauli measurements on stabilizer states, which is actually not 
as restrictive as it initially seems. Most quantum states which are used as resources 
for measurement-based quantum computation are examples of stabilizer states. More 
than this, the set of measurements usually includes Pauli measurements plus only 
one additional non-Pauli measurement. Thus the situation we consider is shy of full 
measurement-based quantum computation by only a single non-Pauli measurement. 
Admittedly, Pauli measurements on stabilizer states can be efficiently simulated and 
are by themselves probably not even universal for classical computation [1]. On the 
other hand, extending our results for Pauli measurements beyond stabilizer states 
opens the possibility of universal quantum computation as well [18, 31]. Hopefully, 
the progress made here will inspire results beyond stabilizer states, or on stabilizer 
states with this extra measurement included. 

We break up the material that follows into three main areas of study. First, 
in Chapter 2, we review the stabilizer formalism in great detail. This material is 
essential since we want to eventually study quantum correlations resulting from out- 
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comes of Pauli measurements on stabilizer states. Even readers familiar with the 
stabilizer formalism should read the background material in Chapter 2, as much of 
it cannot be found in standard texts. This chapter culminates with the presentation 
of stabilizer-state graphs, a novel and useful graphical representation of the stabi- 
lizer formalism introduced in Ref. [12]. These graphs are inspired by graph states, 
which are a class of stabilizer states privileged with a natural graphical representa- 
tion, and the fact that any stabilizer state is local-Clifford equivalent to some graph 
state. In fact, a stabilizer state is usually local-Clifford equivalent to several graph 
states. Representing these local Clifford operations by decorations on the nodes of 
the graph, then, produces a graphical representation of all stabilizer states. Since 
a given stabilizer state is typically local-Clifford equivalent to several graph states, 
there are generally several graphs that represent the state. However, we show that 
only two graph transformations are sufficient to generate all the graphs representing 
a stabilizer state. As a result, these graph transformations also yield all graph states 
that are local-Clifford equivalent to a given stabilizer state. 

Certain transformations, namely Clifford operations and Pauli measurements, 
take stabilizer states to other stabilizer states. Since both of these stabilizer states 
can be represented by a graph, it follows that the transformation can be represented 
by a graph transformation. Aside from decoration changes, it turns out that all local 
Clifford operations are generated by a single kind of graphical transformation called 
local complementation. With regard to Pauli measurements, the post-measurement 
state is described by a sequence of steps which take the initial stabilizer graph to the 
post-measurement stabilizer graph. In the case of individual Pauli measurements, 
meaning those that act on only one node of the graph, the transformation always 
results in that node being disconnected from the graph, possibly preceded by a graph 
transformation. 

The next main subject matter comes in Chapter 3. This chapter introduces the 
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concepts of locality and realism, which say that all measurement results have pre- 
established values that do not change once particles are spatially separated. Quan- 
tum mechanical correlations violate local realism. Local realistic correlations are 
exactly those correlations achievable classically and, in our case, they serve as the 
basis for communication strategies. Communication is then used to quantify the 
difference between classical and quantum correlations. 

Chapter 3 contains some novel results in the study of local realism, particularly in 
the area of Bell inequalities. These are correlation inequalities that must be satisfied 
by any local and realistic model. In an experimental setup to test local realism, par- 
ticles are put into an entangled quantum state and then projective measurements are 
made on the them. Over many trials, the measurements are found to be correlated 
in such a way that cannot be reproduced under the conditions of locality and real- 
ism. Given a set of correlations, thought of as being obtained from an experimental 
setup, we describe a linear program that accomplishes one of two tasks. First, if 
the correlations can be described in a local and realistic model, the program outputs 
such a model. Second, if the correlations violate local realism, the program outputs 
a Bell inequality which the correlations violate. 

Unfortunately, the runtime of the aforementioned linear program is exponential 
in the number of measurements used in the experiment. We use this as motivation 
to derive a novel technique for generating a large class of Bell inequalities applicable 
to arbitrary numbers of parties and measurements. These Bell inequalities result 
from the unique perspective that local realistic models can be mocked up by using 
random variables. We can derive the allowed correlations for random variables, and 
Bell inequalities arise from testing if the correlations are within this allowed range. 

The third and last subject area we study, in Chapter 4, is a melding of the 
previous two chapters to study the stabilizer state violation of local realism. First, 
we study local realistic descriptions of Pauli measurements on stabilizer states. We 
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do not want to include all local realistic descriptions of Pauli measurements, but 
we also do not want to limit our study to those local realistic descriptions that 
reproduce the correlations of quantum mechanics. Therefore, we end up looking at 
probability-preserving local realistic descriptions of Pauli measurements on stabilizer 
states. These are local realistic descriptions which yield correlations that are the 
same as quantum mechanics, up to a sign. The idea is then that these tables can be 
fixed rather simply. We introduce a standard method for constructing probability- 
preserving descriptions and show that any probability-preserving local realistic model 
can be constructed in the standard way. 

With probability-preserving local realistic descriptions serving as a basis, we can 
then supplement them with classical communication. Following the pioneering work 
of Toner and Bacon [92], and Tessier et al. [90], the intent is to use classical com- 
munication to reproduce the observed correlations present in Pauli measurements 
on stabilizer states exactly. The simplest kind of communication model we consider 
restricts communication between neighbors in the stabilizer-state graph and insists 
that nodes in identical local environments behave identically. The first restriction 
is called nearest-neighbor communication while the second is called site invariance. 
Surprisingly, such a model is capable of reproducing all global measurement corre- 
lations, global measurement correlations being those correlations involving all the 
individual Pauli measurements present in a given Pauli product. It is interesting to 
note that in the ubiquitous two party scenario, global measurement correlations are 
the only correlations considered. Therefore, it seems that the interesting features of 
quantum correlations come from submeasurements of a global measurement. These 
are correlations arising from the consideration of restricted sets of individual Pauli 
operators in the global Pauli product. 

To reproduce all submeasurement correlations, it is necessary to relieve the con- 
straints we placed on our model. In particular, we look at what happens to models 
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for two special classes of graph states, one-dimensional and two-dimensional cluster 
states, when these restrictions are lifted. Cluster states are simply graph states where 
the underlying graph is a square lattice. We find that nearest-neighbor communica- 
tion is too strong a restriction to reproduce submeasurement correlations for either 
one-dimensional or two-dimensional cluster states. Site invariance, on the other hand, 
seems to be more revealing. Site invariance is too restrictive for two-dimensional clus- 
ter states, but we explicitly provide a site-invariant model for the one-dimensional 
cluster state. This is interesting because two-dimensional cluster states are universal 
for measurement-based quantum computation while one-dimensional cluster states 
are not. 

In general, we show that the communication distance, defined as the number of 
successive edges through which communication can be sent, must scale linearly with 
the girth of the graph state in order for a communication model to reproduce all 
submeasurement correlations on all graph states. With this in hand, we know a 
model that successfully reproduces all submeasurement correlations must not be site 
invariant and must let communication span the entire graph. Indeed, under these 
conditions, we do generate a successful model. 

As previously mentioned, much of the material in Chapter 2 preceding the discus- 
sion of stabilizer graphs cannot be found in standard texts or lecture notes. Also, the 
results regarding stabilizer-state graphs is entirely new. Stabilizer graphs, and the 
action of Clifford operations on them, are introduced in Ref . [ 1 2] . The material on 
the graphical description of Pauli measurements is presented in Ref. [11]. Material 
involving the graphical representation of stabilizer codes has not yet been written 
up for publication. In addition, both the linear program and random-variable Bell 
inequalities covered in Chapter 3 comprise novel unpublished work. The unique ap- 
proach to generalized probabilistic theories in that chapter is unpublished as well. 
Chapter 4 contains entirely new results. The material covering LHV tables for Pauli 



Chapter 1. Preliminaries 

measurements on stabilizer states has not been written up for publication. However, 
much of the later content in the chapter is an expanded exposition of the published 
material in Ref. [11]. 

In the rest of this chapter we talk about some conventions and notation we use, 
as well as briefly reviewing some of the central background in linear algebra and 
quantum mechanics. After all, understanding of the content in this thesis requires 
first an understanding of vectors and matrices, and then an understanding of basic 
quantum theory. Those readers already familiar with linear algebra and quantum 
mechanics can safely disregard most of the rest of this chapter, excepting the next 
section on some conventions that we choose to follow. Another purpose of this 
chapter is to introduce some notation we use later, since much of the notation, 
although standard, may be unfamiliar to some. The next section, below, includes 
all of the notation we use which is not considered standard practice. After that, the 
rest of this chapter is designed to highlight aspects of basic quantum theory central 
to the chapters that follow. It is therefore advised that even readers experienced in 
quantum information skim the material to make sure it is fresh. 



1.2 Conventions 

Our first point of business is to discuss the conventions we use throughout this work. 
All of these conventions have been set up to help clarify the exposition of the material. 

First, it is our hope that readers are able to go through each chapter without 
getting bogged down by overly complicated or tangential results. For this reason 
complicated proofs, or proofs of tangential topics, are put into an appendix. This 
allows us to state results and focus on their application rather than dwelling on their 
truth value. 
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We also find that sometimes the result of the calculations in a section is ob- 
scured by the calculations themselves. When this is the case, the result itself will be 
highlighted as such. 

Result. This is the result of our calculations. 

This way, the purpose of all our calculations is absolutely clear. Other important 
results are in the form of theorems, lemmas, and corollaries and so are already 
highlighted. 

Now, a word on some notation we use which, though not standard, is both useful 
and natural for presenting the material herein. We find ourselves often performing 
algebra using modular arithmetic. There are even times when modular arithmetic 
using two different values appears in the same equation. To help clarify what sort of 
modular arithmetic is being used we use a few techniques. 

First of all, the symbol "=" means that quantities are identical without using 
modular arithmetic. When we do use modular arithmetic, the two types we encounter 
are modulo 2 arithmetic and modulo 4 arithmetic. When an equation is meant to 
hold modulo 2 we write the symbol " =2" . In other words, 

X = y (mod 2) <^ x =2 y. (1.1) 

Similarly for equations that hold modulo 4 we have 

X = y (mod 4) -^ x =4 y. (1.2) 

For example, we have that =2 2 and 3 =4 —1. The final situation we encounter 
is that modulo 2 quantities sometimes appear in equations that involve arithmetic 
that is otherwise modulo 4. Our notation when this situation arises is thus: any 
expression with an overbar is meant to be taken modulo 2. For example we have 



that 1 + 1 + 5 + 2 =4 -1, since 1 + 1 + 5 (mod 2) = 1 and 1 + 2 = 3 =4 -1. In 
contrast, notice that 1 + 1 + 5 + 2 7^4— 1. 
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Column Vector 


Row Vector 


Binary Vector 


a'' 


a 


Real Valued Vector 


a 


a'^' 


Pure Quantum State 


\a) 


(a 



Table 1.1: There are three contexts in which vectors are used and so different vector 
notation is used for each context. Binary vectors are row vectors and they are rep- 
resented by an unadorned letter. Regular real valued vectors are given the standard 
column vector notation with an arrow over the letter. Quantum states are column 
vectors written in Dirac notation, otherwise known as bra and ket notation. 



Finally, a word on vector conventions. Much of the relevant modular arithmetic 
is formulated in terms of binary vectors, i.e., vectors whose entries are either 1 or 
0. We would like for these vectors to not be confused with ordinary real valued 
vectors, which we also use. To make things worse, quantum states are also typically 
represented by vectors. To keep these vectors straight, we employ specific notation 
in each case. A binary vector, a, is represented by an unadorned letter. To make 
the distinction even more clear, binary vectors are always written as a row vector; 
in order to get the column vector, one must write a^ . Normal real valued vectors 
are taken by default to be column vectors, and we write them with an arrow over 
the letter such as a. Finally vectors representing quantum states are also column 
vectors, and for these vectors we use Dirac notation. Table 1.1 summarizes the 
different notation for vectors in these three contexts. 



1.3 Quantum theory of qubits 



The quantum states we consider are those of qubits, meaning two-dimensional quan- 
tum systems. Two-dimensional quantum states encompass a very special class of 
quantum systems, yet there is enough richness in these states to occupy this entire 
thesis and more. Quantum information is also often expressed in terms of qubits be- 
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cause of the analogy with classical computers which operate using two-dimensional 
classical systems called bits. Indeed, the term qubit derives from a shortened form 
of quantum bit. 

We review only the basics of qubits below; it is assumed that any reader is 
familiar with both quantum mechanics and linear algebra. For those interested, a 
more thorough overview of the material here can be found in Ref. [70] . The purpose of 
the following subsections is only to highlight and review those aspects of the quantum 
theory of qubits particularly relevant to the following chapters of this thesis. 

1.3.1 Single qubits 

If a qubit is in a pure state, then we can represent its state with a two-dimensional 
vector with complex entries. Because of the connection to classical bits, we write our 
basis for the two-dimensional vector space as 

-|1)- (1-3) 

Any pure state of a qubit can be written as a column vector, |^), which must be a 
linear combination of our basis vectors, i.e.. 



a\0)+P\l). (1.4) 



The state of a qubit is really just an object that tells us what outcomes are pro- 
duced when a measurement is made on the system. We also encounter other ways 
to represent states, with the purpose of all these representations being to predict 
measurement outcomes on the system. 

A projective measurement of any 2x2 Hermitian operator. A, can be made on this 
state. According to the postulates of quantum mechanics, the measurement outcome 
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is one of two possibilities corresponding to the two real eigenvalues of A. Once the 
measurement is carried out, the post-measurement state is one of the eigenvectors 
of A, which one depending on the observed measurement outcome. Throughout this 
thesis, we take the nondegenerate measurement outcomes to be +1 and —1. If this 
seems strange, just think of multiplying A by a constant and adding to it a multiple 
of the identity, /, so as to make the eigenvalues equal to ±1. This does not change 
the post-measurement state, meaning the eigenvectors of A, and so it can be thought 
of as simply relabeling the two different measurement outcomes as ±1. 

The probability, p±, for a measurement on a system in the state {ip) to produce an 
outcome ±1 can be inferred from the average value of the measurement, (iplAlip) = 
(A), since p± = (1/2)(1 ± (A)). This is how the state gives rise to measurement 
outcome probabilities. If we do not know the state exactly but only know that 
it is some state {ipj) with probability pj, then the average value is {iIjj\A\iIjj) with 
probability Pj. Thus, 



(A) = J]p,(V^,|A|^,) = Tr 



Tr [Ap] . (1.5) 



j 

The 2x2 matrix p allows us to assign a state to systems about which we have incom- 
plete knowledge, assigning probabilities, p± = (1/2)(1 ± Tr[Ap]), to measurement 
outcomes as best they can be. 

Finally, the state must assign probabilities such that all the measurement outcome 
probabilities sum to 1. The outcome probabilities summing to 1 can be expressed as 
(/) = 1, since both outcomes in this case are 1 and hence (/) = p+ + P-- For a pure 
state vector this means that (iplip) = 1, so that 

mi^)={a* /^* ) I " J =\a\' + m' = l, (1.6) 

and more generally we also have that Tr [Ip] = Tr [p] = 1. 



11 



Chapter 1. Preliminaries 

1.3.2 Projective measurements 

The measurements we talked about in the previous subsection, which are actually 
the only kind of measurements we ever consider, are projective measurements. They 
are called projective measurements because the state is projected onto an eigenstate 
of the observable, A. Thus the post-measurement state is an eigenstate of A, and so 
any further measurement of A must yield the same outcome previously obtained. 

A projection operator P, which in general can be of any finite dimension, is 
such that one application may change the state, but subsequent applications do not, 
i.e., P^ = P. By rewriting this condition as P{P — /) = we see that projection 
operators can only have eigenvalues 1 and 0. The number of entries with a 1 is the 
number of linearly independent vectors unaffected by the projection operator. These 
vectors form a basis for the subspace onto which P projects, since P\ip) is in the 
subspace for any lip). The dimension of this subspace is simply calculated by taking 
Tr [P], as this sums up all the eigenvalues of P. 

If a projective measurement of A gives an outcome of ±1 then it means the 
post-measurement state is given by 

-7^1^- , ' -liiimh (1.7) 

where P± = (1/2) (/ ± A) is a two-dimensional projection operator. If we want to 
measure another operator, WAU, unitarily equivalent to A, then we get a post- 
measurement state 

^ \l±U^AU)\^) = U^-^^ ^ ■-(/±A)|0), (1.8) 



2' " ' /aiUt^ ama\ 2 



'mi±WAUm ' V(</'li(/±^)l0) 

where |0) = U\iIj). This is a fundamental result. In words, measuring U'^ AU on \il)) 
is equivalent to measuring A on U\il)) and then transforming the post-measurement 

state by [/"'". 
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1.3.3 Pauli matrices 

Any 2x2 matrix can be written as a linear combination of the identity and the Pauli 
matrices. These are written in the |0), |1) basis as, 

A' = I ° M , y = ( ° "' I , and Z = P " . (1,9) 

The Pauh matrices are traceless, Hermitian, and unitary. They are sometimes written 
as X = (Ji, Y = a2, Z = a^, and may include / = ctq- The three non-identity Pauli 
matrices are sometimes combined into a column vector, <? = ( cri o"2 cts ) ? whose 
entries are operators. There are many important properties of the Pauli matrices, 
for instance see Ref. [2 ], but we restrict our discussion to those aspects relevant to 
our studies. 

The Pauli matrices are exceptionally useful in the study of qubits. The fact 
that Pauli matrices form a basis for 2 x 2 matrices means that quantum states and 
operations on quantum states can both be expressed in terms of them. For a quantum 
state p we have 

p=^(/ + a-a), (1.10) 

with \a\ < 1, and where the vector a is the Bloch vector defining the state. Pure 
states are those with \a\ = 1. Likewise, a Hermitian operator. A, with eigenvalues 
±1 can be written as 

A = a-a, (1.11) 

where this time |a| = 1. In both cases, Hermiticity ensures that a has real entries. 

1.3.4 Multiple qubits 

Interesting features arise when we consider a system consisting of multiple, say n, 
two-dimensional subsystems. The state of the overall system, if in a pure state, is 
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represented by a vector with complex entries. Since each subsystem can be in one 
of two orthogonal states there are at very least 2" orthogonal states of the overall 
system, one for each of the subsystem configurations. Thus the overall system can 
be taken to be 2" dimensional. The general state of a system of multiple qubits is a 
linear combination of these orthogonal basis states. 

Given this, the way to combine the vector spaces of multiple systems mathemat- 
ically is quite intuitive. The operation that combines the vector spaces, called the 
tensor product, is denoted by the symbol "(8>." It combines vector spaces of dimen- 
sions (Ia and ds to yield a space of dimension dAds- The essential property it must 
satisfy is that for subsystem operations A and B and subsystem states {ip) and |0), 
it should not matter whether A and B act before or after we combine the systems, 
i.e., 



(A®5)(|^)®|0))=A| 



B\ 



:i.i2) 



Given this requirement, there are multiple ways of representing the tensor product 
of two matrices, of which vectors are a special case. Our choice is to multiply each 
entry of A by all the entries of B. Thus we have. 



A^B 



ciiiB CI12B 
0'2iB a22-B 



an 



021 



^11 &12 

^21 ^22 

hu hi2 

^21 ^22 



ai2 



022 



&11 ^12 

&21 ^22 

&11 hi2 

^21 ^22 



/ aiihii aiihi2 012^11 012612 \ 

011621 0'llb22 012621 O12622 

021611 021612 022611 O22612 

\ O21621 021622 O2262I O22622 / 



;i.i3) 
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Repeated application of this rule can be used to combine an arbitrary number of 

systems. 

We use various notation for the tensor product. A subscript on a matrix, such 
as Aj, generally means the tensor product of all identity matrices, /, with a single A 
matrix in the jth spot. Mathematically, we have Aj = I^---^I^A^I^---^I. 
For instance we can write A^ B = AiB2. It is clear when this notation is not being 
used, as terms will appear such as Ai ® B2. When this happens, we know that Ai 
does not mean A^ I. There is some special notation for tensor products of vectors 
as well. For the standard basis, we almost never write the tensor-product symbol 
preferring the shorthand |00) = |0) ® |0). Lastly, we indicate n copies of a state {ip) 
by writing [ip)'^"- = [ip) ^ ■ ■ ■ ^ {ip), and likewise for matrices. 



1.3.5 The Pauli group 

The multi-qubit extension of Pauli matrices leads to a very important class of unitary 
operators known as the Pauli group, which is as useful for states of multiple qubits 
as Pauli matrices are for single qubits. More specifically, we talk about the Pauli 
group on n qubits and denote it as P„. This group is of central importance to the 
stabilizer formalism, and therefore merits an in depth discussion here. The Pauli 
group on n qubits is defined as all n-fold tensor products of the Pauli matrices X, 
Y, and Z, along with the identity / and phases ±1, ±i. An element of V5 is, for 
example, X ® il ® —Z ® iY ® —il. Of course, one can always pull off the overall 
phase and write —iX ® I ® Z ®Y ® I instead. 

The Pauli group is, first of all, a group which means it is a set of elements, along 
with an associative multiplication operation, that satisfies three properties. The first 
is that it contains the identity operator. The second is that the product of any two 
elements of the group is also a element of the group. Finally the third property is 
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that there must exist, for each element of the group, an inverse within the group. 
For any group there exists a subset of elements that generate the group. That is, any 
group element can be written as a product of these generators. Moreover, elements 
in the generating set can be chosen to be independent, in that no element in the 
generating set is the product of other elements in the set. This is simply because a 
dependent element in the generating set can be thrown out without losing the ability 
to generate the whole group. A subgroup of a group is a subset of group elements 
that by themselves constitute a group. 

We will now prove that the Pauli group is, indeed, a group. First of all it contains 
/ ■ • ■ (S> / by definition, which proves the first group property. As for closure under 
multiplication, two elements oiVn, g = gi^- ■ -^g-n and h = hi®---®hn multiply to 
gh = gihi (S> ■ ■ -^gnhn- The fact that this is also in the Pauli group then follows from 
from the case n = 1 which is easily verified. Finally, given g = gi® ■ ■ ■ ® gn & Vn, 
it follows, again from the n = 1 case, that g"^ = gl ■ ■ ■ gj^ ^ 'Pn- Moreover, 
d'^g = 9i9i (8 ■ ■ ■ ® gl^gn = I also follows from the n = 1 case and so there exists an 
inverse in the group for all elements. 

First, we mention some basic properties satisfied by members of the Pauli group. 
Since Pauli matrices are traceless, this implies that all Pauli group elements are also 
traceless. From the unitarity of Pauli matrices, we also conclude that all members 
of the Pauli group are unitary as well. Elements of the Pauli group that have overall 
phase ±1 are Hermitian, while those with overall phase ±i are anti Hermitian. As 
a result, the Pauli group splits into two parts. Those elements, g G Vn, with overall 
phase ±1 satisfy g^ = I and those with overall phase ±i satisfy g'^ = —I. One 
consequence is that elements with ail overall phase have eigenvalues ±1, while the 
other Pauli group elements have eigenvalues ±i. 

A remarkably useful property of the Pauli group is that any two elements either 
commute or anticommute. Verification of the n = 1 case is left to the reader. Once 
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this is established we can define the following function that takes elements of Vi as 
input. 

{0 if Oj and hj commute 
(1.14) 
1 if gj and hj anticommute. 

This function allows us to reverse the order of multiplication through the identity 
9j^j ~ (— l)'^''*'^^'*/^jfi'j- Now we consider two elements of Vn, g = Qi® ■ ■ ■ ® Qn and 
h = hi ® • • • ® hn- We calculate, 



\g^ h]= gh-hg = hg f{-l)^l=if(9j>hi) _ A ^^^^ 
{g, h} = gh + hg = hg ({-l)^l-if(9,A,) + i 



:i.i5) 



We now see that g and h either commute or anticommute and they commute if and 
only if Yl^=ifi9j'^j) i^ even, which is to say they contain an even number of Vi 
elements that anticommute. 
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The stabilizer formalism 



Representing quantum states as a linear combination of basis states can quickly 
become futile as the dimension or number of systems grows. However, if we limited 
ourselves to studying states that can be written as a linear combination of just a 
few basis vectors, we would miss most quantum states and therefore many of the 
interesting features of quantum mechanics. Many methods have been developed 
to study more complicated quantum systems [68, 73, 87], with one very important 
example being the stabilizer formalism [19]. The stabilizer formalism has found uses 
in quantum error correction [49], measurement-based quantum computation [, ,,], the 
study of quantum entanglement [35, 43, 95], as well as in other areas [27, 34]. 

We review the stabilizer formalism with special emphasis on those topics relevant 
to our needs. Much of what follows is focused on developing mathematical tools 
used to augment the stabilizer formalism to better understand its power. Probably 
the most significant result in this regard is the development of stabilizer graphs in 
Sec. 2.3. Much of the work before this section is simply background on the stabilizer 
formalism, with a few novel results interspersed. 
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2.1 Stabilizer states 

If a state {ip) is a +1 eigenstate of some operator A, i.e., A\iIj) = {ip), we say that 
\^) is stabilized by A. For example, |+) = (l/v2)(|0) + |1)) is stabilized by the 
X operator. In fact, |+) is the only state stabilized by X. If we consider a tensor 
product of two Pauli matrices, then we find that two linearly independent states 
are stabilized by the operator. As an example, we find that X ^ X stabilizes both 
(l/v^)(|00) + |11)) and (1/V2)(|01) + |10)). Stabilizer states are so called because 
they are states that are stabilized by tensor products of Pauli matrices. Thus we 
need to explore a bit the stabilizing properties of these operators. 

We can generalize the previous examples to an n-fold tensor product of Pauli 
matrices, g, by considering (1/2) (/ + g). This operator has the property that, 

iW + g)Y = l{I + 2g + g^) = \{I + 2g + I) = \{I + g), (2.1) 

and so is a projection operator. The eigenvectors of (1/2) (/ + g) are the same as 
those of g, with |(/ + 5')|V') = IV^) for those \i/j) stabilized by g and |(/ + g)]^^) = 
for the other eigenvectors of g. These properties together imply that (1/2) (/ + g) 
is the projector onto the states stabilized by g. Using this projector we determine 
the dimension of the subspace stabilized by g to be Tr [(1/2) (/ + g)] = 2"^^, which 
follows from Tr [(?] = 0. Thus an n-fold tensor product of Pauli matrices stabilizes 
2n-i linearly independent states. 

Let us consider once again the states stabilized by X X. If we also impose the 
condition that the state be stabilized hj Z ^ Z, then the only state we are left with 
is (l/\/2)(|00) + |11)). This is the unique state stabilized by both X 0X and Z Z. 
Generally, let g and h be two commuting n-fold tensor products of Pauli matrices. 
They, of course, need to commute in order to have common eigenvectors. As before, 
we can determine the number of linearly independent states stabilized by g and h by 
considering the projector onto the subspace stabilized by both operators. The trace 
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of the projector yields the subspace dimension 

'1 



Tr 



1(1 + 9)1(1 + h) 



Tr 



^{I + g + h + gh) 



2"-^ (2.2) 



since gh is traceless as long as g ^ h. Thus 2""^ states are stabilized by both g and h. 
We see that this same projection argument will give that k independent, commuting 
n-fold tensor products of Pauli matrices stabilize 2^~^ linearly independent states, 
with each new independent operator cutting the dimension in half. The independence 
property means that no generator can be written as a product of other generators, 
up to a sign, and this is equivalent to products of operators being traceless. Thus 
the subspace dimension ends up being (l/2^)Tr [/] = 2""^. 

Everything still works the same way if we choose our operators to be n-fold 
tensor products of Pauli matrices with overall phases ±1. For example, the state 
(l/v^)(|00) - |11)) is the unique state stabilized by -X ® X and Z ® Z. Other 
elements of P„, however, are not valid as stabilizer operators because an overall 
phase of ±i means that there is no +1 eigenvalue of the operator. On a more 
fundamental level, a product of Pauli matrices with overall phase ±1 is Hermitian, 
and is therefore an observable, while those with phase ±i are not. The property 
that stabilizer operators are also observables is crucial for many applications of the 
stabilizer formalism. That being said, a subspace stabilized by tensor products of 
Pauli matrices, with overall phases ±1, is called a stabilizer subspace. The case where 
the number of operators equals the number of qubits is special. 

Definition 1. A stabilizer state on n qubits is the simultaneous +1 eigenstate of n 
independent, commuting elements of the Pauli group. 

It is now appropriate to introduce some terminology related to stabilizer states. 
These terms arise from a connection between stabilizer operators and group theory. 
The key observation is that for two Hermitian Pauli group elements, g and h, which 
stabilize some state, the product gh also stabilizes the same state. This is simply 
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the closure under multiplication property of groups. Clearly / stabilizes all states, 
and also Hermitian Pauli group elements are their own inverses. Thus the Pauli 
group elements that stabilize a stabilizer state form a group. Once n independent, 
commuting Pauli group elements that stabilize a state are specified, products of those 
operators generate the group. For this reason these n operators are called stabilizer 
generators. The whole group is called the stabilizer group, and an element of the 
stabilizer group is called a stabilizer element, stabilizer operator, or simply stabilizer. 
This is an abelian group by construction, meaning all elements in the group mutually 
commute. 



Let us consider an example of a stabilizer state. We stated that (l/\/2)(|00) — 111)) 
is the simultaneous +1 eigenstate of —X ® X and Z ® Z. It is simple to verify that 
—X X and Z ^ Z commute, since an even number of individual Pauli matrices 
in the tensor products anticommute. Thus (l/\/2)(|00) — |11)) is a stabilizer state 
with stabilizer generators —X Cg> X and Z ^ Z. The stabilizer group for this state 
consists of the operators / (S> /, —X ^X, Y ^Y, and Z ^ Z, all of which stabilize the 
state. The operator Y ^Y is a stabilizer element, but not one of the listed stabilizer 
generators. 



We generally specify a stabilizer state by listing n stabilizer generators for its 
stabilizer group. This set is not unique, however, since given two generators, g and 
h, replacing g by gh generates the same group. It is possible to define a canonical 
set of generators for any stabilizer state, although doing so requires a more powerful 
description of the stabilizer formalism. The rest of this chapter is devoted to devel- 
oping the mathematical machinery used to better understand and manipulate this 
important class of states. 
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2.1.1 Binary representation 

The binary representation of Pauli matrices associates a two-dimensional binary 
vector r{aj) with each Pauh matrix aj, where r(/) = ( ), r{X) = ( 1 ), 
r{Y) = ( 1 1 ), and r{Z) = ( 1 )• The vertical line is simply a means to help 
visualize the two halves of the vector, a schematic that will be helpful later. The 
first element of r{aj) can be thought of as specifying the X part of aj, since it is 1 
only for X and Y = iXZ. Likewise, the second element specifies the Z part, since 
it is 1 only for Z and Y = iXZ. In the sense that Y = iXZ, Y has both an X part 
and a Z part, while I has neither. By taking ri{aj) and ^2(0"^) to mean the first and 
second entries of r{aj), we can reconstruct the Pauli matrix from the binary vector, 

aj = fii''jy2{<rj)xM'^,)z''^^''^\ (2.3) 

The binary representation of an element of Vi specifies that element up to the 
overall phase of ±1 or ±i. We will worry about this phase later, but ignoring the 
phase does have a few benefits. For example let us consider the product of two Pauli 
matrices, using the binary vector representation. 

_ ^ri(CTj)r2(o-j)+T-i(o-fc)r2(o-fc)/_-|^\ri(crfc)r2(crj)j^ri(o-j)j^ri((Tfe)2-r2(crj)2'»'2(crfc) (2.4) 
_ ^ri(o-j)r2(crj)+ri((Tfc)r2(o-fe)+2ri((Tfc)r2{o-j)j^ri(o-j)+ri((Tfc)2'r2{o-j)+»'2(o-fe) 

Since the binary representation ignores the overall phase, we find from Eq. (2.4) that 



r{ajak) =2 r{aj) + r{ak). (2.5) 

Multiple applications of this formula yields a similar formula for the product of an 
arbitrary number of Pauli matrices. 
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The binary representation is sometimes called the symplectic representation for 
the following reason. Define a symmetric 2n x 2n matrix, A, as 





A = I ^— I , (2.6) 

where, as before, the lines split the matrix in half. The symplectic inner product of 
two binary vectors is determined by the n = 1 version of this matrix as, 

f if [aj,ak] = 
r(crj)Ar (ak) = ri(crj)r2(o-fc) + r2{aj)ri{ak) =2 < (2.7) 

[^ 1 if {crj,o-fc} = 0. 

The last equality can be checked simply by plugging in all the binary vectors. Thus 

the symplectic inner product is a natural one for binary vectors representing Pauli 

matrices because it evaluates whether the matrices commute or anticommute. Notice 

that r(o"j)Ar"^(crfc) is simply the function f{aj,ak) introduced in Eq. (1-14). 

The binary representation of Pauli matrices is generalized to an arbitrary n- 
fold tensor product, g, as follows. Since we can represent each Pauli matrix in 
the tensor product as a two-dimensional binary vector, we can represent the entire 
n-fold tensor product as a 2?7,-dimensional binary vector, r{g). It is conceptually 
helpful to break up r{g) into its first n components and last n components so that 
^(9) = ( ^1(5') T"2ig) )■ The binary vectors ri{g) and r2{g) generalize the notion 
of ri{aj) and r2{crj) for Pauli matrices. The jth component of ri{g) is 1 only if the 
jth Pauli matrix in gf is X or Y, while the jth component of r2{g) is 1 only if the 
jth Pauli matrix in g is Z or Y. In other words, the jth and (n + j)th entries of 
r{g) give the binary representation of the jth Pauli matrix in the tensor- product 
decomposition of g. 

Let us consider a few examples. For two qubits, the binary vector has dimen- 





sion 2x2 = 4. Some examples in this case are r{X X) = ( 1 1 

r{X0Y)=(l 1 1 y and r(X O Z) = (" 1 

amples the first Pauli matrix is always X, that means the first entry is always 1 



1 ) . Since in these ex- 
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and the third entry is always 0. In the n-qubit case, we have r(/) as the zero vec- 



1 ), and 



tor, r(X®") = (^ 1 ... l|o ■■■ o),r(Z«")=(o ■■■ 
rCK®") as a vector of all I's. 

The way we construct r{g) may, at first, seem a bit peculiar. Perhaps the obvious 
way to generalize the binary vector for a tensor product of Pauli matrices would 
have been to simply place the binary vector for each Pauli matrix side by side. The 
reason for our construction, however, becomes more obvious when we consider the 
symplectic inner product of r{g) and r{h), 

r{g)Ar'^{h) = ri(5f)r^(/i) +r2{g)r'({h) 

n (n ii\ 

= Y.{^n{gUr,{h)], + [r,{g)Un{h)],). 

From Eq. (2.7), we see that the quantity [ri{g)]j[r2{h)]j + [r2{g)]j[ri{h)]j is or 1, 
modulo 2, depending on whether the jth Pauli matrices of g and h commute or not. 
Therefore, following the discussion after Eq. (1.15), we see that g and h commute if 
r{g)Ar'^{h) is even, and anticommute otherwise. This gives the following result. 

Result 1. For two Pauli group elements, g and h, we have that 

, if \g,h] = 
r{g)Ar' [h) =, <^ ^^' ^ (2.9) 

1 if{^,M = 0. 

Eq. (2.5) also holds for arbitrary n-fold tensor products of Pauli matrices g and 
h. To see this we write, 

n 

g = ^E;=i[n(s)b-[r2{9)], K>) J^[ri(g)],^[r2(9)], 

(2.10) 
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Now we can generalize Eq. (2.5) by considering the product, gh, 



(2.11) 



X (g)"_ x[''i(»)1^Z[''2(9)]jX['"i('')1jZ[''2W]j 

^ gh = zE"=i{[n(9)],[r2(g)], + [ri(h)],[r-2(h)],+2[r-2(3)],[ri(/,)],} 

X (g)"_ X[''l(»)l-'+t''l('')]jZt'"2(s)lj+I'"2('')l^ 

This leads us to our next result about binary vectors. 

Result 2. For two Pauli group elements, g and h, we have that 

r{gh)=,r{g)+r{h). (2.12) 

This formula is important to the discussion of stabilizer states that follows. 

2.1.2 Generator matrix 

So far we have discussed how to represent an n-fold tensor product of Pauli matrices 
as a 2n dimensional binary vector. Our task now is to use this formalism to alge- 
braically describe stabihzer states. Since a stabilizer state is specified by listing n 
stabilizer generators, we can describe the stabilizer state by writing the 2n dimen- 
sional binary vector for each generator and then using these vectors as the rows of 
a n X 2n matrix. This matrix is called the generator matrix for the stabilizer state, 
and is given the symbol G, 



( r{g,) \ 



G 



(2.13) 
\ — r{gn) — j 

where we have labeled the n stabilizer generators as gi,. . . ,gn- This matrix only 
specifies the generators up to an overall phase of ±1, a deficiency addressed in Sub- 
sec. 2.1.3. 
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Let us once again consider the stabilizer state (1/a/2)(|00) — |11)), with stabihzer 
generators gi = —X ® X and g2 = Z ® Z . These generators have binary represen- 



tations r{gi) = (l 1 ) and r(^2) = ( 
matrix for this state is, 



1 1 
G ' 







1 1 



1 1 ) . Therefore a generator 



(2.14) 



Not every n x2n binary matrix is a vahd generator matrix; there are two impor- 
tant properties generator matrices must satisfy. First, the stabihzer generators for a 
stabihzer state must be independent, meaning that no product of them can produce 
±/. Since Eq. (2.11) teUs us that taking a product of stabihzer generators corre- 
sponds to adding their binary vectors, no rows of G can sum to zero and therefore 
G must have fuU rank. The second important property comes from the fact that 
stabihzer generators must mutuaUy commute. We can calculate commutativity by 
taking the symplectic product of stabilizer generators, giving us that 

GAG^ =2 0. (2.15) 

Any n X 2n matrix satisfying Eq. (2.15) and having full rank is a valid generator 
matrix. 

As previously noted, we are free to replace any stabilizer generator by a product 
of itself and another generator. This means that we are allowed to row reduce 
the generator matrix, performing addition modulo 2, without changing the stabilizer 
state it represents. Furthermore, swapping qubits has the effect of swapping columns 
of the generator matrix; swapping qubits j and k is equivalent to swapping column 
j with column k and swapping column n + j with column n + k. These operations 
are sufficient to put the generator matrix into standard form. 

Result 3. To any stabilizer state one can assign a generator matrix in standard 
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form, 

/ A 
G 





B . 

(2.16) 



This standard form was obtained by performing row reduction, including swapping 
qubits, to the left half of G. The left half of G need not have full rank, and so it is 
possible that this process terminates with rows of zeroes as shown. After this, we 
row reduced the right half of G making the lower right corner the identity matrix, 
which is possible since G has full rank. Finally, Eq. (2.15) implies that A and A'^ 
appear as indicated and that B = B^ , i.e., i? is a symmetric matrix. As an example, 
we see that the generator matrix in Eq. (2.14) is already in standard form. 

Any set of stabilizer generators that produce a generator matrix in standard 
form are called canonical stabilizer generators. It turns out, and this is surprisingly 
fortunate, that the standard form of a generator matrix is not unique, and hence 
neither are canonical stabilizer generators. However, in Sec. 2.3 we will learn not 
only how to generate all sets of canonical stabilizer generators for a given stabilizer 
state, but also how to efficiently test whether two generator matrices in standard 
form correspond to the same stabilizer state. Before turning to this, though, we 
finish up the binary representation of stabilizer states by finally dealing with the 
sign ambiguity inherent in the formalism. 

2.1.3 Products of Pauli group elements 

The reason for this subsection is a bit subtle, so we begin with an example to illustrate 
the problem. Suppose we had been given the stabilizer generators for (l/-\/2)(|00) — 
\11)) a.s gi = Y ® Y and g2 = Z ® Z . The resulting generator matrix looks like, 

,1111, 
G=\ , (2.17) 

11 
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and is not in standard form. This is not a problem because we can row reduce this 
matrix by adding the second row to the first. The resuh is the generator matrix 
in Eq. (2.14), which is in standard form. By converting the two rows of the matrix 
in Eq. (2.14) into tensor products of Pauh matrices, we get canonical generators 
g[ = X ^ X and g2 = Z ^ Z. Of course, the correct canonical generators are 
g[ = —X (g) X and g2 = Z ® Z, which means we lost the overall phase on the first 
generator. The purpose of this subsection is therefore to determine the overall phase 
of the canonical generators g' as a function of r{g') and the given generator matrix 
G. 



In order to determine the exact form of a canonical stabilizer generator, g', all 
we need now is the overall phase of g', henceforth called the sign of g' and denoted 
sign{g'). Since this sign is either +1 or —1 we will also sometimes find it convenient 
to write sign(f7) = (—1)*^, with the sign of generators being simplified to sign{gj) = 
(—1)*^ . To determine the sign of g' we will assume that during row reduction we kept 
track of which stabilizer generators were multiplied together to obtain g'. In other 
words, we know the values aj G {0, 1} such that g' = gl^ ■ ■ ■ (?."", where gi, . . . ,gn are 
the given stabilizer generators. We know the canonical generator g' must be of that 
form because g^ = I and the given generators commute. Thus we will determine 
sign(5f') in terms of the binary vector a whose jth entry is aj. 



So let us assume we know the values of the aj such that g' = g'l' ■ ■ ■ gn" ■ Iii that 
case we want to generalize the product in Eq. (2.11) to our situation. Let's begin by 
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just writing down the product, recognizing the dot product that appears. 

. . . fjT.]=i[ri{9n)]j[r2{gu)h (g)"_^ X[''i(9n)b-^[''2(9n)]jy" 
= ^E"=i{«i[''i(9i)]j[''2(3l)]j+---+a„[ri(g„)]j[r2(5i„)]j} 

X (^"_ X"i[''"i*^^i)]jZ"^[''2(^i)]j' ■ ■ . X""f''l(^"^l■'Z""l''2(3n)]J 
_ jairi{gi)r^{gi)^ ^a„ri{g„)r'^ {g„) 



(2.18) 



In the above formula we do not include a possible sign of the generators. However, 
if the generators do have a sign, we can simply put those signs in at the end. Now 
we want all the X's together and Z's together. It is clear that X'*'''''^^^'')!-' needs to 
commute through the k — 1 Z operators that come before it, giving us 

„/ ^ jT,]=i'^jn{gj)r^(gj)ti\T,]^iT,i<k"-kai[r-i{gk)h[r2{gi)]j 

X (^"'_ X'^i[''i(^i)l-'+'"+""[''i(^"^l^Z"i[''2(fi'i)]j+'''+""[''"2(9")]j (2.19) 

= jT,]=ia-jri{gj)r^(gj)+'2J2k>l^kairi{gk)rJ(gi) ^^_ X^^^^'^ Z^'^^^'^ 

If we make the simple observation that 

n 

g' = sign(^') X f^ia'^^ia') (g) x'-^^a') z'-^ia') ^ (2.20) 

then we can immediately pull off sign(5f') from Eq. (2.19) to get 

Sign(o') = i^l=i'^3'^'^^3j)r^(gj)+'2T,k>l"-kairiigk)r^(gi)j^-ri{g')r'^(g')^ (2.21) 

Ostensibly we are now done. However, Eq. (2.21) can be simplified considerably with 
some basic, albeit tricky, algebra. Given the unsightly nature of Eq. (2.21), we now 
carry out this simplification to obtain the final result of this subsection. 

Using the notation, introduced in Sec. 1.2, that quantities underneath an overbar 
are understood to be evaluated modulo 2, our first task is to simplify 2 X]fc>/ ^k^'i^i {dk) 
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Tg ((?;). This only needs to be known modulo 2, and since ^ . f{j) =2 ^- f{j), it is 
sufficient to simplify the quantity 2Y2f^^iakairi{gk)r2{gi). This allows us to use the 
fact that the generators commute, i.e., from Eq. (2.9) 



rii9kMi9i)=ri{9i)rU9k), (2-22) 

to write, 

ri{9k)rU9i) = \ [n{9k)rl{9i) + n{9i)rl{9k)) (2.23) 



Using Eq. (2.23) we can simplify '2.Ylik>i'^kO'iri{9k)r2{9i) because 



2 Efc>« akain{gk)r^{gi) = 'Ek>i{(^kairi{gk)r^{gi) + akairi{gi)r'^{gk)} 



= E"fc=i akairi{gk)r^{gi) - XlLi akriigk)r2i9k)- 
Plugging Eq.(2.24) into Eq. (2.21) yields the moderately nicer looking formula. 



(2.24) 



sign(g') = ^SJ'=laj{'•lfe)'■i'fe)-n{^i)r^((;j)}^EJ,fe«j«fe''l(9J)'•i'(9fc)^-'■l(^')'-|'(9') (2.25) 

The best way to write this, though, is to recognize that we can split up the given 
generator matrix G as 

G = ( Gi I Gs ) . (2.26) 

The combination T = G1G2 is really the quantity of interest in these formulas for 
sign(5f'). The reason for this is that 

n n 

Tjk = $^(Gi),7(G2)« = Yy^^3i)W2{9k)]i = r,{g,)rl{gu) (2.27) 

1=1 1=1 

So in terms of this T matrix we have that 

sienfoO = i^l='i-°-i^'^n~'^n\^j,k"-3"'kTjkl-'n{g')r'^{g') 

- -2 (2.28) 

where Ta is the T matrix restricted to those rows and columns corresponding to 
nonzero entries of a. In that last step we also used the fact that, for the binary 
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=;2 



matrix T^, we have ^jj^{Ta)jk = Tr T^ . To truly complete the formula, we should 
note that if the generators themselves have a sign then the signs of the generators 
needed to obtain g' should be multiplied in front of Eq. (2.28). This is the final result 
of this subsection, and it is a quite useful result of which we often make use. 

Just to make sure we did not get lost in all the algebra, let us use this formula for 
its originally intended purpose. The example we began with was the generator matrix 
given in Eq. (2.17). We wished to find the canonical generators for the stabilizer state 
specified by this generator matrix. After row reducing this matrix, we found that 
the canonical generator matrix for the state was the one given in Eq. (2.14), with 
canonical generators g'^ = gig2 and g2 = g2- In this case it is simple enough to 
multiply out gig2 to determine the first canonical generator. Since this is not always 
so simple, we could alternatively deduce this sign from Eq. (2.28). We note that 
a = ( 1 1 ) and we calculate 



Ta= \ = and T, = . (2.29) 






Hence we get, using ri{g[)r2 (g'l) = 0, that sign((7']^) = {^lTa-o+o]j^-o ^ j2 ^ _-^ 
While calculating the formula in this simple case was overkill, for more complicated 
cases it has to potential to save considerable effort. 

2.1.4 Canonical generator sets 

Whenever we deal with stabilizer states we always describe them in terms of canonical 
generators, and the last few sections were devoted to showing how one can find 
canonical generators. This subsection will illustrate what one can do once canonical 
generators have been found. 

Given a generator matrix, G, in canonical form, we define the canonical dual 
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generator matrix to be 
/ 



H 




(2.30) 

■ 



It is simple to check that the dual generator matrix H satisfies the constraints in 
Subsec. 2.1.2 to be a generator matrix. It is also straightforward to check that 
GAH^ =2 I. This means that a dual generator, hj, commutes with all but one of the 
generators, Qj, with which it anticommutes. Dual generators with this property exist 
for any set of stabilizer generators, but for our canonical generators we have the dual 
generators taking a particularly simple form. The importance of dual generators is 
that they reveal the connection between a and r{g), where g = g^^ ■ ■ ■ g'^\ Before, we 
assumed that we knew a for a given g, but now with dual generators we can derive 
a simple relation. 



The fact that g = g^^ ■ ■ ■ g'^" gives us r{g) once we know a, since 

G 

This final equality turns out to be invertible, since from GAH^ =2 I we have 



is) =2 Y^o-A'Ji) = oG = ( a ) . (2,31) 






(2.32) 
i G \ f \ 



where A is the matrix in Eq. (2.6). Hence we can find a once we are given ri^g^ by 

(a ) =2 T{g)k ( G^ if^ ) A, (2.33) 

producing the following result. 

Result 4. For any stabilizer element g = g1^ ■ ■ ■ g'^", where gi, . . . ,gn are canonical 
generators, we have that 

a = r{g)AH^. (2.34) 
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We should note that this equation holds without using modulo 2 arithmetic. The 
reason is because all entries of r{g)AH^ are already either 1 or due the form of H. 
The second condition from Eq. (2.33), that r{g)AG^ =2 0, simply states that g must 
commute with all the stabilizer generators in order to be a stabilizer element. 

Now the only missing piece of information about g is its sign. However, this is 
remedied through the application of Eq. (2.28), with sign(5f) = {—lY^ and sign(f7j) = 
(— 1)**^ We calculate that 



T= 1= \^T= . (2.35) 






This leads us to our next result. 

Result 5. Let g = g^^ ■ ■ ■ (?"", where gi, . . . ,gn are canonical generators whose signs 
are specified by the entries of a binary vector s. Then, 

R n 

2sg=4 2as^ + a( \ a^ - ri{g)r2 (g) 





/ 5 

2rig)isHAf + nig) 

■ 



(2.36) 



In the final equality we simply used Eq. (2.34) to eliminate a in favor of r[g). 

2.2 Unitary operations and measurements 

Unitary operations and measurements are the typical actions on quantum states. 
Indeed any quantum operation be thought of as a unitary operation followed by a 
projective measurement whose outcome is ignored [71]. For stabilizer states, we want 
to consider those operations and measurements that are natural in the context of 
the formalism we have developed. We begin this section by discussing the unitary 
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evolution of stabilizer states and then we end the section by learning about projective 
measurements. 



2.2.1 Clifford operations 

This subsection concerns itself with the issue of unitary operations on stabilizer 
states. Given a stabilizer state {ip) and a unitary operator U, the state resulting from 
the application of the unitary to the state, i.e., U\iIj), is generally not a stabilizer 
state. In other words, general unitary transformations take us out of the stabilizer 
formalism. In order to stay in the stabilizer formalism, we can only consider unitary 
transformations that take any stabilizer state to another stabilizer state. The group, 
and it is a group, of unitary operations that does this is called the Clifford group. 
Elements of the Clifford group are called Clifford operations , Clifford elements, or 
Clifford unitaries. 

We now work towards a better understanding of the Clifford group. Since stabi- 
lizer states are described via stabilizer generators, we must know how stabilizer gen- 
erators transform under Clifford operations. So let \il)) be a stabilizer state, f/ be a 
Clifford unitary, and (?i, ...,(?„ be stabilizer generators for \%l)). From these generators 
we can construct the stabilizer generators which stabilize U\il)) by conjugating the 
stabilizer generators of \il)), yielding the operators UgiU\ . . . , UgnW^ . It is simple to 
verify that these are stabilizers of the new state since {UgjW)U\ilj) = Ugj\ip) = U\iIj). 
The generators UgjW commute since [UgjU\UgjW] = U[gj,gj]W = 0. They are 
also independent since if YYi=ii^9j^'')"'^ ~ ^-^ then it follows that n?=i fi*/ ~ ^-^ 
as well. Thus we have indeed constructed a set of stabihzer generators for the new 
state. 

Since we want U\iIj) to be a stabilizer state, we require that UgjU'^ be in the 
Pauli group. This means that the unitary operators U map elements of the Pauli 
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group, Qj, to elements of the Pauli group under conjugation. In group theory terms 
we say that the operators, U, normalize the Pauh group. Indeed, the Clifford group 
is defined to be the normalizer of the Pauli group, meaning it consists of all unitaries 
that normalize the Pauli group. This group of unitary operators is a natural one to 
consider in the context of stabilizer states. 

As an example, we point out that the Pauli group elements themselves are in the 
Clifford group. For instance, since two Pauli group elements g and h either commute 
or anticommute, we have that ghg'' = ±h. There are, not surprisingly, other elements 
of the Clifford group as explained below. 

2.2.2 Generators of the Clifford group 

The Clifford group is generated by the generators of the local Clifford group plus 
one additional two-qubit Clifford operator [ ]. Local Clifford operations are those 
elements of the Clifford group that act on a single qubit. Put another way, the 
local Clifford group is the normalizer of Vi, and hence maps Pauli matrices to Pauli 
matrices under conjugation. Local Clifford operations are generated by, 

(i o\ 

and S = , (2.37) 

meaning that any local Clifford operation can be written as a product of H and S 
operators. We give these operators the names Hadamard operator, for the H matrix, 
and phase operator, for the 5* matrix. 

The H and S operators are indeed elements of the local Clifford group. We can 
verify this by calculating, 

HIH = /, HXH = Z, HYH = -Y, HZH = X, 

(2.38) 
SIS^ = L SXS^ = r, SYS^ = -X, and SZS^ = Z, 
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where we use the fact that W = H to ehminate W . In words, H interchanges X 
and Z, giving Y a phase, and S interchanges X and Y, possibly adding a phase, and 
leaves Z unchanged. In the binary representation, Eq. (2.38) becomes 



ri(aj) r2{aj) 
ri{(Tj) r2{(Tj) 



S 






ri a,- 



ri((Tj)+r2((Tj) 



(2.39) 



for any Pauli matrix aj with j = 0, 1, 2, 3. Eqs. (2.38) and (2.39) verify that H and 
S are in the local Clifford group. 

We will not prove that H and 5* generate the local Clifford group. For a proof of 
this fact the interested reader is referred to Refs. [38, 49, 50]. The basic idea behind 
the proof is to show that the local Clifford group contains 24 elements and that H 
and 5* generate them. The reason we omit the details is that there is nothing new 
this work would contribute to the proof; it would simply be a restatement of the 
referenced papers. 

Refs. [38, 49, 50] also prove that only one more unitary operation is needed to 
generate the entire Clifford group. There is some freedom in which extra unitary to 
add, but we choose this extra unitary to be the two-qubit controlled-Z operation, 

/ 1 \ 
10 



Cr 




(2.40) 



10 
\ -1 / 

The two qubits acted on by this transformation can be split up into a control qubit 
and a target qubit. If the control qubit is in the |0) state, / is applied to the target 
qubit. Only when the control qubit is in the |1) state is the Z operation applied 
to the target. General controlled operations act in the same way, with a controUed- 
U operation being written as U. The controlled-Z operation is special in that it 
happens to be symmetric in which qubit is the control and which is the target. 
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2.2.3 Pauli measurements 

When a measurement described by a Hermitian matrix M is made on a stabilizer 
state we want the post-measurement state to also be a stabilizer state. We will see 
that the measurements satisfying this property correspond to tensor products of the 
identity and the Pauli matrices. Certainly, since the post-measurement state is an 
eigenstate of the operator, the operator we measure must be a tensor product of 
the identity and the Pauli matrices in order for the post-measurement state to be a 
stabilizer state. What we show below is that, in fact, the post-measurement state 
is a stabilizer state every time such a Pauli product is measured. We could also 
consider the negatives of tensor products of the identity and the Pauli matrices, but 
we neglect this case because the results are the same as the case we consider only 
with the probabilities of the two outcomes interchanged. 

So, suppose we want to measure a Hermitian Pauli group element, M, with sign 
+1. If we suppose that all the stabilizer generators for the initial state {ip) com- 
mute with M then ±M must be a stabilizer element, since adding M to the list 
of generators yields n + 1 generators for an n-qubit state. We know that adding 
independent generators cuts the stabilized dimension in half, so M cannot be inde- 
pendent of the existing generators. Thus M = (— l)™5f™^ ■ ■ ■ (7™", with m, rrij G {0, 1} 
determining how the generators multiply to produce M. This gives us that M\iIj) = 
(— 1)™(?™^ ■ ■ ■ Qn"]^) = (— 1)™!'?/'). A measurement of M therefore gives result (—I)'" 
with probability (^|(l/2)(/ + (-1)'"M)|^) = (1/2)(1 + (-1)™(-1)'") = 1, and the 
post-measurement state is (l/2)(/ + (— 1)'"M)|-?/') = lip). 

The case we want to consider now is when M anticommutes with at least one 
generator g. The first point of business is then to find an appropriate set of generators 
for the stabilizer state. For each given stabilizer generator, Qj, M either commutes 
with Qj or anticommutes with it. If M anticommutes with a generator Qj, then 
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MqjQ = —QjMg = QjqM and so M commutes with the product QjQ. Therefore if 
we replace Qj with QjQ in our generator hst, we end up replacing a generator that 
anticommutes with M with a generator that commutes with M. In this way we can 
come up with a set of generators such that only one anticommutes with M. 

Now we claim that these generators, Qj ^ g, that commute with M are valid 
stabilizer generators of the post-measurement state. For, consider the action of Qj 
on the un-normalized post-measurement state, gj{l/2){I + {—1)"^M)\iIj) = (l/2)(/ + 
(—l)'^M)gj\ip) = (1/2)(I +(—l)'^M)\ip), where (—1)"^ is the measurement outcome. 
Now we only need one more independent stabilizer generator to completely specify 
the post-measurement state. It should be clear now that M performs this role. 
M must certainly be independent of all the generators because it anticommutes 
with g. We also know that (— 1)'"M stabilizes the post-measurement state because 
(-l)™M(l/2)(/ + (-1)™M)|V^) = (l/2)((-l)'™M + I)\ij) = (l/2)(/ + (-1)™M)|V^). 
Thus if we obtain the measurement result (—1)'", we find independent, commuting 
stabilizer generators for the post-measurement state to be (— 1)™M and all gj ^ g. 
Furthermore, it follows that {iplMlip) = {ip\Mg\i;) = {ip] - gM\ip) = -{iplMlip) => 
(iplMltp) = 0, so the measurement outcomes are ±1 with equal probability. 

To summarize, the action of measuring a Hermitian Pauli group element, M , on 
a stabilizer state is as follows. Find a set of generators for the state such that at most 
one anticommutes with M. If all the stabilizer generators commute with M, then the 
outcome is certain and the state is unchanged. If one generator anticommutes with 
M, then replacing that generator by ±M gives a new set of generators that stabilize 
the post-measurement state. In this case the measurement outcome is completely 
random. 

An important special case is to consider the action of individual Pauli measure- 
ments on stabilizer states. These are precisely the kinds of measurements considered 
in Chapter 4. If we measure a Pauli operator cr = X,Y, or Z on qubit j the mea- 
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surement operator looks like M = I^---^I®a^I®---®I, with a occurring in 
the jth spot. This operator commutes with exactly those generators whose tensor 
product decomposition have either a cr or an / in the jth spot, and anticommutes 
with all the others. The generators for the post-measurement state, except for ±M, 
can all be chosen to have an I in the jth spot. 



2.2.4 Measurement correlations 

The remarkable thing about quantum measurements is that they exhibit correlations 
that are stronger than those allowed classically. We will see why this is true in Chap- 
ter 3, but for now we simply examine some basic properties of quantum measurement 
correlations. 

We can only speak of correlations for two commuting Hermitian operators, A and 
B. In this case the correlation between A and B comes from calculating the average 
of the product AB, which is Hermitian from our assumptions. The average value 
obtained by measuring the product AB can be calculated by measuring A and B 
separately. This value is equal to measuring A and measuring B, which can be done 
since A and B commute, and then multiplying together the results and averaging 
over many trials. This is not generally the same as multiplying the averages of A and 
B, i.e., mathematically (AB) ^ {A){B). In the case that they are the same, then A 
and B outcomes are independent of each other and so A and B are not correlated. 
The generic case allows for A and B to be correlated. 

To illustrate this, we show how individual Pauli measurements can build up the 
statistics of any other measurement that is a tensor product of Pauli matrices. Con- 
sider, for example, a measurement of a ® cr' with a, a' = X, Y, or Z. The probability 
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that the outcome is +1 is 

') UIj) = (ilAl (16^1 + a6^I)(T6^T + T6^cr') Uh) 

(2.41) 



II (J ® / + a ® a') IV^) = (V^ll (/ ® / + a ® /) (/ ® / + / ® a') 



+(v^l i (/ ® / - a ® /) (/ ® / - / ® or') iv^; 

The quantity on the right is expressed in terms of various outcomes for measurements 
of the individual Pauh operators a®/ and I®a' . Specifically, we took the correlation 
between both outcomes being +1 and both outcomes being —1, since these are the 
cases in which their product is +1. Thus the statistics oi a ® a' can be obtained by 
looking at the average of the product of the individual Pauli measurements involved. 
However, we stress that while measuring a ® I and I ® a' gives the same statistics 
as measuring a ® a', the post-measurement states in these two scenarios are quite 
different. 

To review, the essential point of this subsection is that one can look at the corre- 
lations between two commuting measurements. This is done by averaging over many 
trials the product of the measurement outcomes, rather than taking the product of 
the average of each measurement by itself. Mathematically, we account for corre- 
lations by taking the average value after multiplying the operators together, rather 
than multiplying together the average values of the individual operators. We also 
found that correlations between measurement outcomes are equivalent to correlations 
resulting from measuring a new operator which is the product of the two commuting 
operators whose outcome correlations are being studied. 

2.3 Stabilizer-state graphs 

The traditional representation of stabilizer states is to list n commuting, indepen- 
dent stabilizer generators for the state. In this section we introduce an alternative 
description in which stabilizer states are represented by graphs [42]. We will see 
that this graphical representation of stabilizer states is immensely powerful, not only 
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calling on results from the binary representation of stabilizer states but also drawing 
strength from the circuit model of quantum computation. 

2.3.1 Graph theory 

There is a considerable collection of terminology associated with graph theory [ ] . 
Fortunately, most of this terminology is intuitive, with only a few exceptions. All 
that aside, learning this terminology is central for being able to manipulate stabilizer 
states through their graphs. Hopefully Fig. 2.1 will help in this regard. 

The simplest kind of graph is portrayed pictorially as points, called nodes, con- 
nected by lines, called edges. Two nodes that are connected by an edge are called 
neighbors, while the neighborhood of node j, denoted Af{j), is a set including all of 
that node's neighbors. The adjacency matrix, so called because its elements are de- 
termined by which nodes are neighbors, and therefore are adjacent to each other, is 
capable of completely describing any graph. The adjacency matrix, F, for an n-node 
graph is a.n n X n binary matrix defined as 

{1 if nodes ?' and k are connected 
(2.42) 
otherwise. 

A 1 on the diagonal of an adjacency matrix is drawn as a self loop, but in such cases 
we still do not consider that node as its own neighbor. Also note that adjacency 
matrices are symmetric, meaning F-^ = F. It is clear that graphs and adjacency 
matrices contain the same amount of information. 

There are a few important operations on graphs that we need to know about. 
The most basic operation is called complementation. Complementing the edge be- 
tween two nodes removes the edge if one is present and adds one otherwise. Note 
that an edge need not exist in order for it to be complemented. From this operation 
we construct two types of local complementation, local complementation on a node 
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Node 
Self loop 
5**=== — / Neighbor 

Edge 




(b) i 







Figure 2.1: The graph in (a) is used to illustrate the basic terms used to describe 
a graph. Each black dot is a node, with node number 2 having a self loop. We 
indicate in (a) one of the neighbors of node 2. Finally, we represent edges simply 
with a line that connects two nodes. In (b) we take the graph in (a) and perform 
local complementation on node 2, the one with a loop. We do not consider node 2 
its own neighbor so it has only three neighbors, nodes 1, 4, and 5. In local comple- 
mentation we complement the edges between these neighbors. This means that the 
edges between pairs of nodes {1, 5}, {1, 4}, and {4, 5} are complemented. The result 
is that the edges between nodes 1 and 5 and between nodes 4 and 5 are removed, 
since these nodes are connected in the graph in (a), and an edge is added between 
nodes 1 and 4 since they are not connected in (a). Finally, in (c), we take the graph 
in (a) and perform local complementation along the edge connecting nodes 2 and 4. 
This is done by first performing local complementation on node 2, then performing 
local complementation on node 4, and then performing local complementation on 
node 2 again. If we had instead chosen to first perform local complementation on 
node 4, then the intermediate stages would have looked different but the final graph 
would have been the same. 
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and local complementation along an edge. Performing local complementation on a 
node complements the edges between all of the node's neighbors. Local complemen- 
tation along an edge is equivalent to a sequence of local complementations at the 
nodes defining the edge. The sequence is as follows: first perform local complemen- 
tation on one of the nodes, then local complement at the other node, and finally 
local complement at the first node again. Local complementation along an edge is 
symmetric in the two nodes defining the edge, so it does not matter at which node 
local complementation is first performed. 

2.3.2 Graph states 

Our scheme for drawing graphs corresponding to stabilizer states derives from a class 
of stabilizer states called graph states, which are privileged with a natural graphical 
representation [54, 55]. Consider a generator matrix G in standard form, Eq. (2.16). 
Let us call the rank of the left n x n matrix - this is the Gi matrix in Eq. (2.26) - 
the left rank of G. If the left rank of G is n, then the canonical form of G is 

B ) , (2.43) 



G = 

with the constraint that B^ = B. Since i? is a symmetric binary matrix, it is the 
adjacency matrix of some graph. Furthermore, if B has all O's on the diagonal the 
graph has no self loops. In the case where the generator matrix has full left rank, 
and where B has all O's on the diagonal, the generator matrix represents a stabilizer 
state called a graph state. We could just as well define a graph state to be any state 
whose generator matrix has full left rank, but for historical reasons graph states are 
defined not to have self loops. A graph state is completely determined by the graph 
with adjacency matrix T = B, and they are clearly a proper subset of all stabilizer 
states. 

From Eqs. (2.42) and (2.43), we can write the canonical stabilizer generators for 
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a graph state as 

g,=X, n Zk. (2.44) 

keAf(j) 

Recognizing that Z{X I) Z = X Z, we can see that graph states are the 
result of beginning with generators gj = Xj and then applying Z between qubits 
corresponding to neighbors in the graph. The generators Qj = Xj correspond to the 
state 1+)®", so a graph state's graph can be regarded as a blueprint to build the 
graph state from a simple initial state. 

2.3.3 Graphs from generator matrices 

We now have all the necessary tools needed to assign a graph to a general stabilizer 
state. Any stabilizer state with a generator matrix of canonical form can be converted 
by local Clifford operations to a state possessing a generator matrix with the form 
in Eq. (2.43). From Eq. (2.39), applying Hadamard operations to the last n — r 
qubits of the stabilizer state, where r is the initial left rank of the generator matrix, 
exchanges columns r + 1 through n in the generator matrix with columns ra + r + 1 
through 2n, so that a generator matrix as in Eq. (2.16) is transformed to 

/ / B a\ 
G=\ \. (2.45) 

\0 I A^ Q j 

The diagonal of B in this generator matrix can then be stripped of I's without oth- 
erwise changing the generator matrix by applying S to offending qubits and using 
Eq. (2.39). The resulting generator matrix has the form of a graph state and cor- 
responds to a stabilizer state that differs from that represented by Eq. (2.16) by at 
most a single H or S operation per qubit. 

The close relationship between graph states and stabilizer states suggests the 
possibility of a graph like representation of stabilizer states. Our approach to such 
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a representation is simply to transform the generator matrix of a stabilizer state 
into that of a graph state, draw the graph thereby obtained, and add decorations to 
each node indicating whether an if or S* was applied to the corresponding qubit in 
the process. Motivated by the standard graph convention that a 1 on the diagonal 
of an adjacency matrix denotes a self loop on that node, we choose to represent S 
operations by self loops. Indicating the application of an H is less clear and indeed 
there are several ways to do it [81, 86]. As there is no obvious choice, we employ 
the use of two different types of nodes. Normally we draw nodes as filled in dots, or 
as solid nodes, so we indicate qubits to which an H was applied as open circles, or 
hollow nodes. We will sometimes use the term fill state to reference the dichotomy 
in node types. 

To summarize, we can graphically represent any stabilizer state with the fol- 
lowing procedure. First, determine the canonical generators for the state, as ex- 
plained in Subsecs. 2.1.2 and 2.1.3, which yield the generator matrix in standard 
form, Eq. (2.16). Second, draw the graph, including self loops, with adjacency ma- 
trix given by 



r = I. , (2.46) 




making solid the nodes corresponding to the rows and columns of the submatrix B 
and hollow the nodes corresponding to the rows and columns of the submatrix 0. An 
example of a stabilizer state and an associated graph is given in Fig. 2.2. 

There are many comments that should be made before proceeding further. For 
one, notice that this procedure does not associate every combination of edges, self 
loops, hollow nodes, and solid nodes with a stabilizer state. Because the submatrix 
in Eq. (2.46) contains only O's, hollow nodes never have loops, and there are no edges 
between hollow nodes. We should also point out, here, that a stabilizer state might 
be representable by more than one graph. This follows simply because a stabilizer 
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Figure 2.2: (a) A generator matrix for a stabilizer state, (b) A canonical-form 
generator matrix obtained from (a) by row and qubit swapping. The dashed lines 
are to facilitate comparison to the standard form for generator matrices in Eq. (2.16). 
(c) The adjacency matrix indicated by (b). (d) The adjacency matrix of (c) with the 
qubit swaps undone, (e) The stabilizer graph associated with (a). In parts (a)-(d) 
the columns have been labeled by the corresponding qubit. In part (e) the nodes 
are labeled by the corresponding qubit. It is clear that the qubit swaps are not 
actually necessary, since we reverse them in the end. The generator matrix in (a) 
can be converted to graph form directly by exchanging columns 2 and 3 on the left 
with the matching columns on the right, an operation that corresponds to applying 
a Hadamard to qubits 2 and 3. In graph form, the adjacency matrix is just the right 
half of the generator matrix. Loops arise from I's on the diagonal of the adjacency 
matrix, and hollow nodes are used to indicate which columns were exchanged between 
the right and left halves of the generator matrix to get the adjacency matrix. Notice 
that there are no edges between hollow nodes, nor are there any loops on hollow 
nodes. 



state might be described by more than one set of canonical generators. We will 
learn to fully characterize all the stabilizer graphs representing the same state in 
Subsec. 2.3.6. 

It is also important to note that the nodes in our graphs are labeled in that each 
node is associated with a particular qubit. Swapping two qubits is a physical opera- 
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tion that generally produces a different quantum state. In our graphs a SWAP gate 
can be described either by relabeling the corresponding nodes or by exchanging the 
nodes and all their decorations and connections while leaving the labeling the same. 
Since the process of bringing the generator matrix into canonical form can involve 
swapping qubits, we must keep track during this process of the correspondence be- 
tween qubits and columns of the generator matrix and thus between these columns 
and the nodes of our graphs. 

Possibly the biggest caveat of which to be aware is that our graphs, so far, rep- 
resent stabilizer states in the same way that generator matrices represent stabilizer 
states; that is, they represent them up to an overall sign of the generators. This is a 
downside that will be remedied in the next section. 



2.3.4 Graphs from quantum circuits 

An alternative perspective of our graphs for stabilizer states is that, in analogy with 
graphs states, they provide a blueprint for building the stabilizer state. Given a 
graph that represents a stabilizer state, if we remove all the decorations, that is we 
fill in all the nodes and remove all the self loops, we obtain a graph state which can 
be prepared as described in Subsec. 2.3.2. Then we can apply an S to all qubits 
corresponding to nodes with a self loop and an H to all qubits corresponding to 
hollow nodes, thereby producing the stabilizer state represented by the given graph. 
Thinking of our graphs for stabilizer states as a recipe for creating the state allows 
us to enrich our graphical formalism through a connection with quantum circuits, a 
useful theoretical tool in the toolbox of quantum computation. An example quantum 
circuit is illustrated in Fig. 2.3.4(a), but for a more complete explanation see Ref. [70]. 
Figs. 2.3.4(b) and (c) list the quantum circuit version of some frequently used single- 
qubit and two-qubit Clifford operations, respectively. 
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Figure 2.3: (a) A simple example of a quantum circuit. Each horizontal line tracks 
the evolution of a single qubit, with time moving forward as the circuit progresses 
towards the right. The state to the left of each horizontal line indicates the initial 
state of the corresponding qubit. A line incident on a box that means that an 
operation, which in the context of circuits is called a gate, is applied to that qubit. 
A single-qubit gate U corresponds to a box with a U inside that only intersects one 
horizontal line. A multi-qubit gate V is also indicated by a box with a V inside, 
except that the box intersects those lines corresponding to qubits on which it acts. 
The meter on the second qubit at the end of the circuit indicates a measurement. The 
output of the measurement is c, which is a classical bit whose evolution is indicated 
by a double horizontal line. Common single-qubit Clifford gates are illustrated in 
(b), with the shorthand for common two-qubit Clifford gates listed in (c). The gates 



a 



in (c) are ^Z, ^X, and ^F, with the last gate being the SWAP gate. The SWAP gate 



swaps two qubits and is equivalent to three X gates as indicated 



Let us write the quantum circuit for the preparation procedure derived from a 
graph representing an ra-qubit stabilizer state. Such a quantum circuit consists of 
three layers of gates applied to n qubits, each initially in the state |0). In the first 
layer, the Hadamard gate, if, is applied to each qubit. This prepares the qubits in 
the state, |+)®", required to create graph states using the method of Subsec. 2.3.2. In 
the second layer, controlled-Z gates, *"Z, are applied between various pairs of qubits. 
Using only these first two layers, it possible to create any graph state. Finally, in 
the third layer, Hadamard gates and phase gates are applied to various subsets of 
qubits depending on whether the corresponding graph nodes are hollow or have a 
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self loop. As a mnemonic, one can think of H being applied to qubits corresponding 
to 'H'ollow nodes, and S to those corresponding to nodes with a 'S'elf loop. In this 
final layer, there are two restrictions, arising from the restrictions on what graphs can 
represent stabilizer states. The first restriction is that a phase gate and Hadamard 
gate never occur on the same qubit. This corresponds to hollow nodes never having 
self loops. The second restriction is that two qubits with iJ's in the third layer are 
never connected by a Z gate in the second layer. This is the restriction that hollow 
nodes are never connected to each other in a graph. Quantum circuits consisting of 
these three layers are deemed to be in graph form. 

We mentioned earlier that quantum circuits would enable us to represent the 
signs of stabilizer generators, and we now demonstrate how this is done. Let us 
begin with a graph state, so that we have the generators in Eq. (2.44). In this case, 
the application of Zk to the graph state changes the sign of generator g^ and leaves 
the other generators invariant. Using 6jk which is 1 ii j = k and otherwise, this 
can be written as 

ZkdjZk = Zk \Xj YY Zi \ Zk = ZkXjZk Y\^ Zi 

= (-1)^^'= Ix, H zA =i-l)'^-g,. 

Thus we can dictate the sign of graph state generators by the presence of Z's in the 
third layer of the quantum circuit; a Z on qubit j indicates that sign^gj) = —1. This 
procedure extends to all stabilizer states since a later application of if or 5 does not 
further change the sign of the generator, 

sign{HkgjHk) = signiX^/" Zj'' g^) = sign(^j), 

; ; (2.48) 

sign{SkgjSl) = signiX/'Y/'gj) = sign(^,). 

Just as the presence of an if or S* in a quantum circuit is represented graphically as 
a hollow node or self loop, for obvious reasons we indicate a Z in the quantum circuit 
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as a minus sign in the corresponding node of the graph. This allows us to represent 
generator signs, and therefore completes our graphical description of stabilizer states. 

Even though our graphs now provide a complete description of stabilizer states, 
a further generalization will prove beneficial in the context of Clifford operations 
and Pauli measurements. Our graphs for stabilizer states already suffer from the 
affliction that multiple graphs can represent the same state, so we generalize our 
graphical representation so that any graph consisting of solid and hollow nodes, with 
or without self loops and signs, connected by arbitrary edges represents some state. 
We call such a graph a stabilizer-state graph, or, more simply, a stabilizer graph, 
while graphs with the property that hollow nodes do not have self loops and are not 
connected to each other are called reduced graphs. Reduced graphs are capable of 
representing any stabilizer state and reduce to the usual notion of graphs for graph 
states. We will fully explore the connection between general stabilizer graphs and 
reduced graphs in Subsec. 2.3.6. 

Now we explain how quantum circuits permit us to associate a stabilizer state 
with any stabilizer graph. To begin with, allowing arbitrary edges in a graph, so 
that hollow nodes can be connected to each other, corresponds to allowing arbitrary 
Z gates in a graph-form quantum circuit. Now it only remains to explain what is 
meant by a hollow node with a loop. We choose this to mean that, in the third layer 
of a graph- form quantum circuit, an 5* and an H gate have been applied, in that 
order. The order is important since SH ^ HS. We have already chosen Z gates to 
precede H gates for good reason, so since Z and S commute it makes sense to have 
them both precede the H gate. 

To conclude this subsection, we concisely state the relationship between stabilizer 
graphs and graph-form quantum circuits. Stabilizer-state graphs and graph-form 
quantum circuits are in one to one correspondence. In the second layer of a graph- 
form quantum circuit, two qubits are linked by a Z gate if and only if they are 
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Figure 2.4: (a) A circuit in graph form, with the layers numbered and separated 
by dashed hues, (b) A stabihzer graph corresponding to the circuit in (a). Z 
gates between qubits are transformed into hnks between nodes, terminating Z gates 
become negative signs, terminating S gates result in loops, and terminating H gates 
are denoted by hollow nodes. 



connected by an edge in the corresponding stabilizer graph. In the third layer of 
the quantum circuit, Z, S, and H gates appear, in that order. A Z gate appears 
if and only if the corresponding node has a minus sign, an S gate appears if and 
only if it has a loop, and an H gate appears if and only if the node is hollow. Both 
stabilizer graphs and graph-form quantum circuits are in many to one correspondence 
with stabilizer states. Fig. 2.4 has an example of a stabilizer graph along with its 
equivalent graph-form quantum circuit. 



2.3.5 Graphical description of Clifford operations 

Having explained our graphical representation of stabilizer states in some detail, 
now we get to see their utility in action. We will examine the usefulness of our 
graphs more in subsections to come, but an appropriate introductory demonstration 
of their applicability is to describe how the graph representing a stabilizer state 
changes when Clifford operations act on the state. Both stabilizer graphs and reduced 
stabilizer graphs have their relative merits, so we include an explanation of how 
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Figure 2.5: Simple circuit identities, which we use to derive more complex identities. 
Identity (a) follows trivially from the steps shown, while identity (b) is easily verified 
in the standard basis. 



Clifford operations effect both types of graphs. This section generalizes the results 
in Ref. [33]. 

Aside from the terminology in Subsec. 2.3.1, there a few other terms we use to 
describe the action of Clifford operations. Two of these terms are flip and advance. 
Flip is used to describe the simple reversal of some binary property, such as the sign 
or the fill of a node. Advance refers specifically to an action on loops; advancing 
generates a loop on nodes where there was not previously one, and it removes the 
loop and flips the sign on nodes where there was a loop. Its action mirrors the 
application of the phase gate, since S^ = Z. The last bit of terminology is our 
sloppy use of language for the sake of brevity. That is, when we say that a gate acts 
on a node we mean that it acts on the qubit corresponding to that node. In fact, in 
what follows, qubit and node are used interchangeably. 

We begin with the action of local Clifford operations on stabilizer-state graphs. 
Since the local Clifford group is generated by H and S", we need only describe the 
action of H and S on stabilizer states. This is accomplished by the following trans- 
formation rules. 

Tl. Applying if to a node flips its fill. 

T2. Applying S* to a solid node advances its loop. 
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T3. Applying 5 to a hollow node without a loop performs local complementation 
on the node and advances the loops of its neighbors. 

If the node has a negative sign, flip the signs of its neighbors as well. 

T4. Applying 5* to a hollow node with a loop flips its fill, removes its loop, performs 
local complementation on it, and advances the loops of its neighbors. 

If the node does not have a negative sign, flip the signs of its neighbors as well. 

These transformation rules can be derived from the circuit identities in Fig. 2.6, 
which rely on the basic circuit identities given in Fig. 2.5. Given an understanding 
of the relationship between circuits and graphs, transformation rules Tl and T2 are 
trivial. Transformation rules T3 and T4 derive from Figs. 2.6(a) and (b), respectively. 

The transformation rules T1-T4 do not generally take reduced stabilizer graphs 
to reduced stabilizer graphs. From Sec. 2.3.3, however, we know that there exists a 
reduced stabilizer graph corresponding to each stabihzer state, so it is always possible 
to represent the effect of a local Clifford operation as a mapping between reduced 
stabilizer graphs. The appropriate transformation rules for reduced stabilizer graphs 
are listed below. 

T(i). Applying if to a solid node without a loop, which is only connected to other 
solid nodes, flips the fill of that node. 

T(ii). Applying if to a solid node with a loop, which is only connected to other solid 
nodes, performs local complementation on the node and advances the loops of 
its neighbors. 

Flip the node's sign, and if it now has a negative sign, flip the signs of its 
neighbors as well. 
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Figure 2.6: Circuit identities relevant to transforming a stabilizer graph under local 
Clifford operations. Both identities are illustrated for the case of four neighbors con- 
nected to a qubit of interest; equivalent expressions hold for other numbers of neigh- 
bors. The proof of identity (a) relies on decomposing '"Y using Y = iXZ, the equiv- 
alence of '^i to S on the control qubit, and the identity in Fig. 2.5 (b). The first step 
in identity (b) uses the fact that I = e-'^l\HSf, which gives SHS = e'^'^XHSH; 
the second and third steps follow from the identities in (a) and Fig. 2.5 (a). The 
identity in Fig. 2.5 (b) generates the thicket of controlled-Z gates in the lower left of 
the final circuits in (a) and (b), thereby giving rise to local complementation in the 
corresponding stabilizer graphs. In (b), the phase shift on the top qubit, e™/'^^ is a 
global phase shift and thus can be omitted. To complete the transformation rules, it 
is necessary to know what happens in both of these circuit identities when there is 
a Z gate on the top qubit immediately after the Z gates. Ignoring overall phases, 
the effect on the third layer of the final circuit is to include an additional Z on all 
qubits. 
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T(iii). Applying H to a solid node without a loop, which is connected to a hollow node, 
flips the fill of the hollow node and performs local complementation along the 
edge connecting the nodes. 

Flip the sign of nodes connected to both the solid and hollow nodes. If either of 
these two nodes has a negative sign, fiip it and the signs of its current neighbors. 

T(iv). Applying H to a solid node with a loop, which is connected to a hollow node, 
performs local complementation on the solid node and then on the hollow node. 
Then it removes the loop from the solid node, advances the loops of the solid 
node's current neighbors, and fiips the fill of the hollow node. 

Flip the signs of nodes that were originally connected to both the solid and 
hollow nodes. If the originally solid node initially had a negative sign, fiip it 
and the signs of its current neighbors, and if the originally hollow node initially 
had a negative sign, fiip the signs of its current neighbors. 

T(v). Applying H to a hollow node fiips its fill. 

T(vi). Applying S* to a solid node advances its loop. 

T(vii). Applying S* to a hollow node performs local complementation on that node and 
advances the loops of its neighbors. 

If the node has a negative sign, fiip the signs of its neighbors as well. 

Of these transformation rules, T(i), T(v), and T(vi) are trivial, and T(vii) is a 
rewrite of T3. To prove the others requires results from Subsec. 2.3.6, in particular, 
equivalence rules El and E2. Specifically, T(ii) is obtained by applying equivalence 
rule El, which gives an equivalent, but unreduced graph, and then applying the 
Hadamard, via rule Tl, which leaves a reduced graph. For T(iii), one first applies the 
Hadamard, via rule Tl, and then uses equivalence rule E2 to convert to a reduced 
graph. In the case of T(iv), one applies the Hadamard, using rule Tl, and then 
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applies equivalence rule El, first to the originally solid node and then to the hollow 
node. A key part of these transformations is the conversion of stabilizer graphs to 
reduced form. A section of Subsec. 2.3.6 explains this process in more detail. 

It is not hard to check that, in using rules T(iii) and T(iv), any other hollow 
nodes that are connected to the originally solid node do not become connected and 
do not acquire loops, in accordance with the need to end up with a reduced graph. 

To complete the Clifford group, we need to present a graphical description of the 
action of controlled-Z gates on stabilizer states. Since most of the interest in Clifford 
operations is in local Clifford operations [32, 55, 56, 57], and since the action of Z 
generates all stabilizer states from a given initial one, we relegate the results relevant 
to Z gates to Appendix A. 

2.3.6 Equivalent graphs 

In this subsection we explain how to test whether two graphs correspond to the 
same stabilizer state, and as a consequence we learn how to generate all graphs 
representing a given stabilizer state. Let us begin with an example of different 
graph-form circuits corresponding to the same stabilizer state. Such an example can 
be found in Fig. 2.6. Applying an additional S'^ gate to the top qubit in Fig. 2.6(b) 
makes the initial and final circuits both in graph form, but the two circuits lead to 
quite different graphs. These two graphs then represent the same stabilizer state 
and so are called equivalent graphs. Two graphs are equivalent exactly when their 
corresponding graph-form circuits are equal. 

We now present two equivalence rules for stabilizer graphs. Applying either of 
the following two rules to a stabilizer graph yields an new one which represents the 
same stabilizer state. We will then use these rules to derive a similar set for reduced 
graphs. The proof that both sets of equivalence rules generate all graphs equivalent 
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to the initial one is given in Appendix B. 

El. Flip the fill of a node with a loop. Perform local complementation on the node, 
and advance the loops of its neighbors. 

Flip the node's sign, and if the node now has a negative sign, fiip the signs of 
its neighbors as well. 

E2. Flip the fills of two connected nodes without loops, and local complement along 
the edge between them. 

Flip the signs of nodes connected to both of the two original nodes. If either 
of the two original nodes has a negative sign, fiip it and the signs of its current 
neighbors. 

The first of these equivalence rules is our previous example, and can be obtained 
by applying an additional S^ gate to the top qubit in the identity of Fig. 2.6(b). 
For the second rule, we need yet another circuit identity. Figure 2.7(a) shows how 
Hadamards can be removed from a pair of connected qubits without S gates. Fig- 
ure 2.7(b) extends this identity to a demonstration of rule E2. 

Equivalence rules El and E2 can be reworked to yield equivalence rules for re- 
duced stabilizer graphs. These rules take reduced stabilizer graphs to equivalent 
reduced stabilizer graphs. 

E(i). For a hollow node connected to a solid node with a loop, local complement on 
the solid node and then on the hollow node. Then remove the loop from the 
solid node, advance the loops of its current neighbors, and fiip the fills of both 
nodes. 

As for signs, follow the sequence in transformation rule T(iv). 
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Figure 2.7: (a) A circuit identity used in (b). The identity is shown with four 
neighbors connected to the top qubit; equivalent expressions hold for other numbers 
of neighbors. The crucial trick here is to decompose SWAP in terms of three X 
gates, (b) A circuit identity demonstrating transformation rule E2. The selection of 
neighbors of the top two qubits is chosen to display all possibilities. The first equality 
in (b) employs the identity from part (a), while the subsequent ones follow from basic 
circuit identities. It is straightforward to verify that edges are changed according 
to local complementation along the edge connecting the top two qubits. Circuits 
illustrating E2 for other initial fill states can be obtained by applying additional 
terminal Hadamards to the circuits in (b). The sign rule for E2 expresses the effect 
of a Z gate on one of the top two qubits before the terminal Hadamard gate. Such 
a Z gate can be pushed through the Hadamard, becoming an X gate; in the final 
circuit, this X gate can be processed using the identity in Fig. 2.5(a). 
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E(ii). For a hollow node connected to a solid node without a loop, local complement 
along the edge between them. Then flip the fills of both nodes. 

As for signs, follow the sequence in equivalence rule E2. 

There is a simple relationship between the two sets of equivalence rules. Equiv- 
alence rule E(ii) is identical to E2 for the case that the two connected nodes have 
opposite fill. Equivalence rule E(i) is simply rule El applied twice: first to the solid 
node with the loop and then once to the hollow node that has acquired a loop from 
the first application of El. The second application of El is needed because the 
resulting graph is not reduced after one employment of the rule. 

Both of these equivalence rules can also be derived by applying two Hadamards 
to a solid node, E(i) handling the case in which the solid node has a loop and E(ii) 
the case in which it does not. Thus equivalence rule E(i) is simply transformation 
rule T(iv) followed by use of rule T(i) to apply a second Hadamard to the originally 
solid node. Likewise, E(ii) is rule T(iii) followed by T(i) to apply a second Hadamard 
to the originally solid node. Notice that both rules preserve the number of hollow 
nodes. 

It is worth summarizing here the results in Appendix B, where the completeness 
of our equivalence rules is proved, because the ideas behind the proof are simple and 
can easily be missed when reading the details. In order to test whether two graphs 
represent the same state, the first step is to use equivalence rules El and E2 to put 
the graphs in reduced form. Then, if equivalence rules E(i) and E(ii) are used to 
put the hollow nodes in the same place, the two graphs represent the same state 
if and only if they are trivially identical. Given that we can test the equivalence 
of any two graphs using our equivalence rules, these rules must then generate all 
graphs equivalent to a given one. The most difficult part of the proof is showing that 
graphs represent the same state if and only if they are trivially identical, once the 
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equivalence rules have been used. Because of this, an example illustrating the proof 
of this step is given in Fig. B.l. Fig. B.2 is also worth looking at since it demonstrates 
how to use the equivalence rules to efficiently test whether two graphs represent the 
same stabilizer state. Also, as a consequence of the proof in Appendix B, we see 
that reduced graphs are special because they constitute all the graphs representing 
a state that use a minimal number of hollow nodes. 

2.3.7 Graphical description of Pauli-product measurements 



As another example of the usefulness of stabilizer-state graphs, we now discuss how to 
graphically describe the effect of Pauli-product measurements on stabilizer states [11]. 
A Pauli-product measurement is a measurement in the eigenbasis of an n-fold tensor 
product of the identity, /, and the Pauli matrices, X, Y, and Z. In other words, a 
Pauli-product measurement corresponds to a matrix M such that 

n 

M = (g)M,, (2.49) 

where Mj = I, X, Y, or Z. If Mj ^ I, we call node j a measured node; otherwise, if 
Mj = I, we say the node is not measured. 

Given such a measurement operator, M, our task is two fold: first, to find the 
probability that a measurement of M on a quantum system in the stabilizer state {ip) 
gives an outcome (—1)™, and, second, to determine the post-measurement quantum 
state of the system, i.e., a post-measurement stabilizer-state graph. This section de- 
scribes a general graphical rule, applicable to the graph that represents the stabilizer 
state, which accomplishes these tasks. 

The difficulty of formulating a Pauli-product measurement rule can be greatly 
reduced by using graph transformation and equivalence rules introduced in Sub- 
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sees. 2.3.5 and 2.3.6. The following three paragraphs describe a sequence of three 
simplifications that can be made to any measurement, thereby restricting its form to 
one more amenable to a measurement transformation rule. With this simplification 
carried out, we then present an applicable graphical measurement rule. 

The first simplification relies on the fact that a measurement where Mj = CjZjC- 
is equivalent to a measurement where Mj = Zj preceded by application of the local 
Clifford operation C- and followed by application of Cj to the post-measurement 
state. For an X in the measured Pauli product M, C = H, and for a F in M, 
C = SH. Thus, the first simplification is to transform the original graph, using 
rules T1-T4, so that on the new graph the measurement becomes a product of Zs 
on the measured nodes. This means that it suffices to determine the effect of Z-tjpe 
measurements, that is, measurements with the property that Mj = I or Z for all j. 
The post-measurement state must be transformed by application of the appropriate 
local unitaries to the measured nodes, i.e., Cj to measured node j; in terms of graphs, 
this post-measurement transformation is handled by rules T1-T4. 

The second simplification is to reduce the graph, as explained in Subsec. 2.3.6. 
For the purposes of our measurement analysis, we are not required to reduce the 
entire graph, but only the measured nodes. After this second simplification, there 
are no loops on hollow measured nodes and no edges between hollow measured nodes. 

The final simplification is to disconnect hollow measured nodes from unmeasured 
nodes using equivalence rules El and E2. Suppose that a measured hollow node 
is connected to an unmeasured node. In the case that the unmeasured node does 
not have a loop, applying equivalence rule E2 to the pair turns the measured node 
solid. In the case that the unmeasured node has a loop, an application of El to the 
unmeasured node gives the measured hollow node a loop. Now one can apply El 
to the measured hollow node to turn it solid. One can verify that in both cases, 
application of the equivalence rules leaves the remaining measured hollow nodes 
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loopless and unconnected to one another. Thus, this last simphfication terminates 
in a number of iterations no greater than the number of measured hollow nodes. 

The end product of these simplifications is a Z-type Pauli measurement on a 
graph in which measured hollow nodes are loopless and unconnected to one another 
and to unmeasured nodes. With the preceding simplifications carried out, we can now 
spell out the graphical description of the measurement. The proof of this description 
is given in Appendix C. 

Our graphical description of Pauli-product measurements is greatly facilitated 
by introducing special sets of nodes in the stabilizer-state graph: Ti is the set of 
hollow nodes, S is the set of solid nodes, Z is the set of nodes with a sign, Ai = 
{j I Mj 7^ 1} is the set of measured nodes, M.s = M.\H is the set of measured 
solid nodes, A^h = A^ fl 7i is the set of measured hollow nodes, and A^se = {j G 
Als I |A1h n A/'(j)| = (mod 2)} is the set of measured solid nodes that have an 
even number of connections to measured hollow nodes. Here A\B denotes the set of 
elements in A that are not m. B, A\^B denotes the intersection of A and B, and |^| 
denotes the number of elements in A. 

When a Pauli measurement is made on a stabilizer state, the outcome is either 
definite or random, with random outcomes giving equal probabilities of 1/2 for the 
two outcomes. Which case applies depends on A^se- The outcome in the determin- 
istic case is specified by 

6=|AljjnZ|, (2.50) 

the number of measured hollow nodes with a sign. 

The result of a Z-type Pauli measurement is as follows. 



1. If AlsE = 0, the measurement outcome is (—1) with certainty, and the state 
is unchanged by the measurement. 
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n>n^'^-'A 



Figure 2.8: The steps that perform a Pauh-product measurement M = I®Z®Z®Z 
on a stabihzer-state graph for the four qubit cluster state. The juxtaposition of 
Pauh operator and node indicates the presence of that operator in the intended 
measurement. We label the nodes in the graphs clockwise starting from 1 in the 
upper left corner. Thus, M = {2, 3, 4}, Ms = Mse = {2, 3, 4}, and Mr = 0. Since 
A^SE 7^ 0, the outcome is random; for the sake of illustration, we take it to be +1. 
Node 2 is taken to be the chosen node, so the edges between its neighbors, nodes 
{1,3}, and the unchosen nodes in A^se, nodes {3,4}, are complemented in step 1. 
The chosen node does not have a sign, and m =2 b ii we assume a +1 measurement 
outcome, so step 2 has the effect of giving a sign to the node in A^se that is also a 
neighbor of node 2, meaning node 3. Step 3 removes all edges involving node 2 and 
then connects it to nodes 3 and 4 while making node 2 hollow. Finally, step 4 has 
no effect since the chosen node does not have a loop. If either node 3 or 4 had been 
picked as the chosen node, the resulting graph could be transformed into this one 
using equivalence rule E2. 



2. If AisE 7^ 0, the measurement outcome, (—1)™, is random, and a graph for 
the post-measurement state can be obtained according to steps 1-4 below. 

To find the post-measurement state when A^se 7^ 0, it is necessary first to pick 
a node, which we call the chosen node, from AlsE- The post-measurement state is 
then obtained by the following four steps. 

1. For each neighbor of the chosen node, complement all of its edges to unchosen 
nodes in AIse- 

2. If the chosen node does not have a sign, flip the signs of all its neighbors that 
are also in AIse; otherwise, if the chosen node has a sign, remove that sign, 
and flip the signs of all other nodes in AIse that do not neighbor the chosen 
node. If m 7^2 b, flip the signs of the chosen node and all its neighbors. 
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3. Remove all edges involving the chosen node, and then connect the chosen node 
to all the other nodes in A^se- Make the chosen node hollow. 

4. If the chosen node has a loop, remove that loop, perform local complementation 
on the chosen node, advance the loops of its neighbors, and if m 7^2 b, flip the 
signs of the unchosen nodes in A^se- 

These steps constitute a complete graphical description for the effect of the measure- 
ment M on the state. Notice that, in step 1, an edge between two nodes that are in 
A^sE and are initially neighbors of the chosen node gets complemented twice, so it 
remains unchanged. Figure 2.8 illustrates the use of this measurement rule for the 
case of a three qubit measurement on a four qubit cluster state. 

2.3.8 Graphical description of single-qubit Pauli measure- 
ments 

Single-qubit measurements are a straightforward but important special case [ ], as 
illustrated, for example, by the use of such measurements in measurement-based 
quantum computation [ ]. While we could easily apply our measurement rules 
above to this special case, we instead formulate a simpler alternate description of 
single-qubit measurements. This description does not generalize to other Pauli- 
product measurements, but it is more simple than our general rule in this special 
case. 

Since there are three possible measurements of Pauli matrices that can be made 
on four possible signless node types, there are 12 possible measurement configurations 
to consider. These are illustrated in Fig. 2.9. We first consider each case in complete 
generality, and then we finish this subsection with a specific example of single-qubit 
Pauli measurements on a stabilizer state. 
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The simplest case to consider is that of a Z measurement on a sohd node without 
a loop, which is case 1 in Fig. 2.9(a). It is simple to determine that in this case 
the node disconnects from the graph [ ]. This simply means that all edges to that 
qubit are removed. Since the measured qubit is in an eigenstate of Z after the 
measurement, it is now represented by a hollow node with a sign determined by the 
measurement result. 

Let us see why this is true. If the measurement is on node j, then the claim is 
that all generators except for Qj have Z or / in the jth spot in the tensor product 
decomposition of the generator and Qj has an X in that spot. Thinking of the circuit 
that creates the stabilizer state, after the second layer we have a graph state and 
we know this is true from Eq. (2.44). Moreover, the only possible Clifford operator 
in the third layer on the jth qubit is a Z. Therefore it must also be true for the 
final stabilizer state as well. Now it is clear that only Qj anticommutes with the 
measurement and so must be replaced by g', = Zj after the measurement. We can 
remove the Z from all generators with a Z in their jth spot through multiplication 
by g'j. In this way we obtain the stabilizer generators for a stabilizer graph obtained 
from the original one by simply disconnecting qubit j. 

We will now show the remarkable fact that a graphical description of all other 
measurements can be determined from this simple case by using equivalence rules El 
and E2 and transformation rules Tl— T4. The plan is to use the equivalence rules to 
represent the stabilizer state by a graph which permits a simple graphical description 
of the given Pauli measurement. 

We now determine the effect of a Z measurement on a solid node with a self 
loop, case 2 in Fig. 2.9(a). The loop is a graphical representation of an S gate being 
applied to a solid qubit without a loop. Now this application of S followed by a Z 
measurement is equivalent to measuring in the SZS'' = Z basis followed by an S 
gate on the measured qubit. We know that a Z measurement on a solid node without 
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(a) (b) (c) 

l.Zon# l.Xon# l.Yon# 

2.Zon^ 2.Xon^ 2.Yon^ 

3.ZonO 3.XonO 3.YonO 

4.Zon(5^ 4.x on (5^ 4. Yon 



Figure 2.9: The different possible Pauli measurements on the different possible node 
types. The sign of the node is irrelevant for classification purposes and is therefore 
omitted. 



a loop disconnects that node, and therefore a Z measurement on a solid node with 
a loop also disconnects the node from the graph. The post-measurement qubit state 
is once again represented by a hollow node. This trick of using the node decorations 
to change the measurement basis is a useful one and is used repeatedly. 

We now finish the last two Z measurement possibilities. In order to graphically 
describe these measurement situations, we use equivalence rules El and E2. In 
the first situation, case 3 in Fig. 2.9(a), suppose a Z measurement is made on a 
hollow node without a loop. We want to flip the fill of the measurement node. This 
can be done via E2 if the node is connected to another node without a loop, local 
complementing along the edge connecting the node and its loopless neighbor. If all 
of the neighbors of the node have a loop, then we must first use El on a neighbor 
to produce a loop on the measurement node. The resulting task is to describe a Z 
measurement on a hollow node with a loop, which is case 4 in Fig. 2.9(a). This can 
be accomplished by applying El to the node to flip its fill and then disconnecting the 
node, since we are now measuring Z on a solid node. The post-measurement node 
is hollow. 

The next cases involve X measurements which are nearly analogous to Z mea- 
surements. The difference is that X measurements on hollow nodes have the effect 
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tu-i-a. •vH 



(b)l. ^Y 



l^-^fZ-^a 'IM 



(c)l. 



• ^^ 



',^'^^ oXX 



Figure 2.10: Examples of X, F, and Z measurements on various node types. See 
the text for a more detailed explanation. 



of disconnecting those nodes from the graph. This follows from recognizing that 
a hollow node, with or without a loop, represents a solid node, with or without a 
loop, followed by an H gate. This H gate changes the measurement basis from X 
to HXH = Z, and the result is equivalent to making a Z measurement on a solid 
node, which we know disconnects that node. In order to make an X measurement 
on a solid node, we need to first make the node hollow. This can be done via equiv- 
alence rules El and E2, analogously to cases 3 and 4 in Fig. 2.9(a). After the X 
measurement the disconnected node is solid. This covers cases 1 — 4 in Fig. 2.9(b). 

Now the only cases that are left are the Y measurements, cases 1 — 4 in Fig. 2.9(c). 
However, if we once again change the measurement basis by writing Y = SXS"^, then 
we can replace the Y measurement by an X measurement preceded by an S gate. 
The action of this S gate is described graphically by transformation rules T2, T3, 
and T4. Now we are left with measuring X on the transformed graph, which can 
be done as described for cases 1 — 4 in Fig. 2.9(b). The post-measurement state of 
the measured qubit is represented by a solid node with a loop. This completes our 
graphical representation of Pauli measurements on stabilizer states. 
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To illustrate the utility of our graphical representation of Pauli measurements on 
stabilizer states we use the examples in Fig. 2.10. We draw attention to the applica- 
bility of our results to measurement-based quantum computation [79] by analyzing 
X, Y, and Z measurements on a 2 x 3 cluster state, which is a graph state corre- 
sponding to a square lattice [19]. We purposefully follow an inefficient measurement 
order for pedagogical reasons. We also always assume the measurement outcome +1. 

We begin with the X measurement illustrated in Fig. 2.10 (a). Representing an 
X measurement on the indicated note can be accomplished by applying equivalence 
rule E2, thereby making the node, along with one neighbor of the node, hollow. Local 
complementation along the edge defined by the two nodes yields the graph on the 
right side of Fig. 2.10 (a)l. The two graphs in Fig. 2.10 (a)l are equivalent in that 
they both represent the 2x3 cluster state. The purpose of using the graph on the 
right, however, is that the post-measurement state can be obtained by disconnecting 
the measured node from the graph, as illustrated in Fig. 2.10 (a)2. 

Next we consider Fig. 2.10 (b) where a Y measurement is made on a solid node 
of the graph. This Y measurement has the same effect as an X measurement on 
the middle graph of Fig. 2.10 (6)1. This X measurement can be carried out by first 
applying equivalence rule El, which results in the far right graph in Fig. 2.10 (6)1, 
and then disconnecting the node, as illustrated in Fig. 2.10 (6)2. Since we changed 
the measurement basis from Y to X, we must apply an 5* gate to the measured node, 
giving it a loop. 

Finally, Fig. 2.10 (c) depicts a Z measurement on a hollow node with a loop. 
Equivalence rule El can be used in this case to turn the node solid, resulting in the 
graph on the right of Fig. 2.10 (c)l. Now that the node is solid, the Z measurement 
can be made by simply disconnecting the node. We present the graph for the post- 
measurement state in Fig. 2.10 (c)2. Note that the unmeasured nodes are in a graph 
state. 
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2.3.9 Graphical description of stabilizer codes 

The stabilizer formalism was originally developed for the purpose of using stabilizer 
codes in quantum error correction. Fortunately, it turns out that stabilizer codes, 
too, can be represented using our graphical formulation of stabilizer states. The 
method we use is an adaptation of the one in Ref. [86], extended in Refs. [84, 85], 
for codes using graph states. We note here that our description is easily adapted 
to subsystem stabilizer codes [58, 78], and even holds hope for being extended to 
codeword stabilized codes [23, 30]. 

A stabilizer code that encodes k logical qubits in n physical ones is defined to 
be a 2^ dimensional stabilizer subspace of n qubits. This means that all stabilizer 
codes are completely specified by listing n — k stabilizer generators on n qubits, 
Qi, ■ ■ ■ ,gn-k- In order to understand how stabilizer codes encode information, con- 
sider the following process. Pick any state in the stabilizer subspace that happens to 



also be a stabilizer state and call it |0 ■ ■ ■ 0), where the overbar means this is an n- 
qubit encoding of a /c qubit state. To define this state we need to list k more stabilizer 
generators, Zi, . . . , Zk, in addition to gi, . . . , gn-k which we already know stabilize 
the state. Given this, we can define a total of 2^ stabilizer states all of which are in the 
given stabilizer subspace. These states are written as |ci ■ ■ ■ c^), for Cj G {0, 1}, and 
are defined to be stabilizer states stabilized by f^i, . . . , gn-k, {—^Y^Zi, . . . , {—lY'^Zk. 
Hence we have defined our encoded basis states and our effective Pauli Z operators 
which put negative signs in the appropriate way on the basis states. It only remains 
to prove our encoded basis states are orthogonal. 



Cl ■■■Cfc 



-ipz, 



Cl ■■■Cfe 



(-iP(-iri(c;---4|ci---c,) (2.51) 



Cj = Cj or ( c'l ■ ■ ■ Ci 



ci---Cfc) = 0, 
and we have completely specified how the k qubits of information are encoded. 
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Since stabilizer graphs are constructed using canonical generators, we begin by 
using the given stabilizer generators to construct a canonical set of generators for 



|0 ■ ■ ■ 0). Doing this, we obtain canonical stabilizer generators, (^'(O ■ ■ ■ 0), which can 
be written in terms of the original generators as 

9',{0 ■ ■ ■ 0) = gT^^^ ■ ■ ■ gl-^^'^zf'^ ■ ■ ■ -ZY'\ (2.52) 

with the binary vector ( a{i) h{j) ) specifying how to obtain the jth canonical 
generator. This binary vector also permits the obtention of canonical generators for 
|ci ■ ■ ■ Cfc), since 

= (-l)K.)c-^;.(0---0). 

Thus all of the 2*^ stabilizer states in our stabilizer code are specified by the same 
canonical generators, up to an overall sign. In other words, all of the states |ci ■ ■ ■ c^) 
in our code are described by the same generator matrix. 

Now we are ready to represent any stabilizer code by a graph. First, using the 
canonical generators g'j (0 ■ ■ ■ 0) draw a reduced stabilizer graph on n nodes that 



represents the state |0 ■ ■ ■ 0). Then draw k more nodes, which we call input nodes, 
giving a total of n + A; nodes. These k additional input nodes correspond to the k 
qubit state being encoded. Next, we want connect the input nodes to the stabilizer 
graph by connecting the /th input node to the jth node of the stabilizer graph if 
bi{j) = 1, i.e., when Zi was used to obtain the jth canonical generator, g'AO---0), 



for |0 ■ • ■ 0). The resulting graph provides a complete description of the basis states 
and logical operators for the given stabilizer code. 

Using this graph we can reconstruct the stabilizer generators for |ci ■ ■ ■ c^) as 
follows. First of all, if q = 1 let us call the /th input node an active input node. 
To determine the jth stabilizer generator, simply begin with (7' (0 ■ ■ ■ 0) and multiply 
it by —1 for each edge connecting the jth node of the stabilizer graph to an active 
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input 




Figure 2.11: A graph that represents all features of the 5-qubit code. This code 
encodes 1 qubit using 5 qubits, and so our graph is made of 1 + 5 = 6 nodes. The 
round nodes, considered by themselves, represent the stabilizer state |0). The square 
input node connects to nodes 1, 4, and 5. This indicates that the canonical generators 
for |1) are obtained from those of |0) by switching the signs of generators 1, 4, and 
5. 



input node. Doing this we find that the generator is multiplied by —1 if and only if 
bi{j) = ci = 1. Hence we have, 

g'^ic ■ ■ ■ c) = (-l)St. M.>.^^.(o . . . 0) = (-l)^0>"^;(0 ... 0), (2.54) 

in agreement with Eq. (2.53). 

Describing in words how to construct graphs for stabilizer codes is difficult. 
Therefore, we choose to illustrate the procedure by constructing a graph for the 
well known 5-qubit code. This code encodes one qubit by using 5 qubits, and this is 
the smallest code that can correct any single-qubit error. According to the 5-qubit 
code, the stabilizer generators for |c), with c = 0, 1, are 

gi = X^Z^Z^X®I, g2 = I®X®Z®Z®X, 

g^ = X(g)I(g)X(g)Z®Z, g4 = Z(g)X®I(g)X(g)Z, (2.55) 

and {-lyz = {-iyZ®Z®Z®Z(^Z. 

By writing out the generator matrix for these generators and row reducing, we 



71 



Chapter 2. The stabihzer formahsm 

determine a set of canonical generators to be 

g[ = {-ly+^X (g)I(g)Z(g)I(g)X, g'^ = I®X®Z®Z®X, 

g'^ = Z®Z®X®I®X, g'^ = {-ly+H ® Z ® I ® X ® X, (2.56) 

and g'^ = {-iyZ®Z®Z®Z®Z. 

As illustrated in Fig. 2.11, we obtain a graph for this stabilizer code by first drawing 
the stabilizer graph for |0) and then drawing an edge from the input node to all 
those graph nodes whose corresponding canonical stabilizer generators contains a 
factor (— 1)"^. 

We often have occasion to think of stabilizer-state graphs as a prescription for how 
to create the state. It turns out this is true for stabilizer codes as well, in the sense 
that the graph gives a set of instructions for how to encode a k qubit state. Fig. 2.12 
shows a circuit corresponding to a graph that represents a stabilizer code, where 
edges to input nodes are also represented by Z gates, followed by a measurement 
in the Z basis on the k input nodes. This circuit is capable of encoding any basis 
state, but fails to encode arbitrary unknown states. This failure of the circuit to 
encode arbitrary states fundamentally arises because quantum circuits are linear in 
the input, while our method of encoding is not linear in the input. However, a slight 
modification of the circuit in Fig. 2.12 does give rise to a circuit that is capable of 
encoding arbitrary unknown quantum states. 

To encode an arbitrary state, lip) = J^a c =o ^ci-cjci ■ ■ ■ Ck), we simply switch 
the measurement basis to the X basis. From Fig. 2.12 and linearity, we know that 
before the measurement the full n + k qubit state is 

1 

y^ aci-cJci---Cfc)|ci---Cfc). (2.57) 

ci---Cfc=0 

Now measure in the X basis by applying a Hadamard and measuring in the Z basis. 
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|o) H 


H 


|o) H 


H 




|Cl) 



Cfc) 



Z, 5, and H gates 



Cl 



Cfc 



Figure 2.12: A circuit that encodes the basis state |ci ■ ■ ■ c^). This circuit is directly 
obtained from the graph that describes the stabihzer code. The Z gates, as well 
as the Z, S, and H gates, are determined by the stabilizer graph part of the code's 
graph. The final k qubits correspond to input nodes, and their controlled-Z gates 
act on those qubits connected to that input node in the graph. If ci, . . . c^ = 0, then 
the controlled-Z gates do not act, and the top n qubits are in state characterized by 
the stabilizer graph, i.e., the |0 ■ ■ ■ 0) state. If some Cj are 1, then the Z gates from 
those qubits deposit a Z onto qubits connected to them. These Z gates flip generator 
signs so that the resulting state encodes |ci ■ ■ ■ c^). Before the measurement, the final 
n + k qubit state is |ci ■ ■ ■ Cfe)|ci ■ ■ ■ Cfc). 



Hadamards on the k input qubits transform the state to 



^ aci...cJci---Cfc) 

ci---Cfc=0 

1 1 



io) + (-irni) 



2fc/2 



ci---Cfc=0 di---dj,=0 



|0) + (-l)^'=|l) 



v^ 7 V v^ 

kn ^ Y^ (-l)''^^«ci...cJci---Cfc)|rfi---4). 



(2.58) 



Now we make the Z measurements and obtain the results Xi, . . . , x^ G {0, 1}, making 
the post-measurement state of the first n qubits 



Y^ (-I)'="^aei...cjci---Cfc). 

ci---Cfc=0 



(2.59) 



This can be made the correct encoded state by applying the logical {ZjY^ operator 
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for all j, since 

k 1 

j=l ci---Cfc=0 

1 A: 

ci--Cfe=0 j=l 

1 1 

= Yl (-l)''''^aci...Cfe(-l)''=^|ci---Cfc) = E aei...cJci---Cfc). 

ci---Cfe=0 ci---Cfc=0 

Thus we see, in analogy with stabilizer states, graphs for stabilizer codes gives a 
circuit that is capable of encoding an arbitrary state. 



2.4 Concluding remarks 

This chapter explored the stabilizer formalism for qubits. Since stabilizer states are 
eigenstates of products of Pauli operators, representing Pauli operators as binary 
vectors gives a powerful description of stabilizer states. We described how to ma- 
nipulate these binary vectors to accomplish various tasks, such as multiplying Pauli 
operators together, applying Clifford operations, and achieving a canonical form for 
the stabilizer generators of a state. 

Inspired by this binary representation of the stabilizer formalism, we introduce 
stabilizer-state graphs. These are simple graphs consisting of solid and hollow nodes 
possibly with loops and signs attached to them. Aided by the quantum circuit 
formalism, we prove that Clifford operations and Pauli measurements are simply 
represented using stabilizer graphs. Moreover, stabilizer graphs are even capable of 
representing stabilizer codes in quantum error correction. Because of the relationship 
to quantum circuits, graphs for stabilizer states and codes have the interpretation 
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that they represent a set of instructions for how to construct or encode the state 
given a small set of elementary gates. 

Applying our graphical representation of stabilizer codes to problems in quantum 
error correction is an appealing direction for future research. In particular, subsets of 
nodes in our graph correspond to stabilizer elements, or syndrome measurements in 
the language of quantum error correction. Knowing which syndrome measurements 
facilitate a simple and effective error decoding strategy is an important open problem 
in quantum error correction. Results in this area may be forthcoming using our 
graphical approach. It also seems that our graphical representation of stabilizer 
codes could be adapted to represent subsystem stabihzer codes and possibly even 
codeword stabilized quantum codes. These advances are on the frontier of known 
error correcting codes and a unifying graphical description that proves useful would 
be desirable. 

Another interesting open problem is the question of which stabilizer states are 
local unitarily equivalent to a given stabilizer state. In the case of local-Clifford 
equivalence, we already know the answer; two stabilizer states are local-Clifford 
equivalent if and only if their graph edges are related by local complementation. 
However, recent progress suggests that local-Clifford equivalence is different than 
general local unitary equivalence [57]. Thus, it may be worthwhile to apply our 
circuit identities to a non-Clifford local unitary operation, thereby extending the 
results here to general local unitary equivalence of stabilizer states. 
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Local realism 



Much of the reason we learned about stabihzer states in such detail was so that we 
could examine their properties from the point of view of local realism. Local realism 
was introduced in Ref. [ i T-] to expose an apparent inadequacy of quantum theory. 
The basic elements of that argument exposed a conflict between quantum mechanics 
and physical intuition, namely that physical reality conforms to locality and realism. 
This conflict remained untested in the laboratory until Bell's theorem [13] provided 
a physical experiment that could be performed. When the experimental results 
were published [5, 82, 98], the results strongly supported quantum mechanics and 
refuted local realism. This result puts quantum mechanics fundamentally at odds 
with concepts of classical physics, and truly sets it apart from all other accepted 
physical theories. 

Simply put, realism is the assumption that all possible measurements are assigned 
outcomes. For instance, quantum mechanics states that one cannot learn both the 
position and momentum of a particle, due to the Heisenberg uncertainty principle. 
Realism states, however, that even though we can only learn the value of position 
exclusive of learning momentum, both position and momentum have well defined 
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Figure 3.1: Realism can be thought of as asserting that there exist tables attached 
to each system which determine the outcomes of all possible measurements. Here 
we imagine that two particles, drawn as filled in circles, interact and then separate, 
bringing with them a table that determines their measurement outcomes. The left 
column of the table is the measurement Aj or B^, and the right column is the 
outcome, if one would decide to make that measurement. In this case outcomes are 
either ±1 in some appropriate units. Each time the experiment is repeated a new 
pair of particles and tables is generated, with the randomness of quantum mechanics 
implying that the tables change in some probabilistic fashion. 



values prior to the measurement giving us one of them. A nice conceptual tool is to 
think of a list of measurements attached to each system, as illustrated in Fig. 3.1. 
This list of measurements has a value associated with each measurement which is the 
outcome that will be given if the measurement is made. Throughout this chapter, 
we always assume that there are two measurement outcomes +1 and —1. 

The idea of locality is just that operations and measurements on one system 
cannot instantaneously affect a spatially separated second system. For our purposes, 
this implies that the tables in Fig. 3.1 cannot change once the particles become 
spatially separated. 

So the task at hand is first to learn about the power inherent in local realistic 
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models to capture correlations in measurement outcomes. When we say correlation 
we mean the average of the product of measurement outcomes. We can then effec- 
tively compare local realistic models with quantum mechanics to find that quantum 
mechanics produces stronger correlations than local realism allows. Specifically, af- 
ter reviewing some important historic work, we find a general approach to testing 
whether a set of experimental results could be reproduced by a local realistic theory. 
This method, however, is limited since the general case of this problem is believed to 
be very hard [ ]. So, we derive next a large class of efficiently computable necessary 
conditions on the experimental results in order for them to be describable by a local 
realistic theory. Finally, we review a mathematical framework for learning about 
all probabilistic theories, which include local realistic models as well as quantum 
mechanics, to better understand the relationship between the two. 



3.1 LHV tables 



The list of measurements and outcomes associated with each separated system in 
Fig. 3.1 is called a local-hidden-variable table, or LHV table for short. The term 
local hidden variable comes from the supposed hidden variables that determine the 
specific values in the table, so that randomness need not be invoked. If we assume 
local realism, then we assume these tables can be constructed and that observed mea- 
surement results are consistently described by such tables. In this section we learn 
how to mathematically construct LHV tables, of which there are two main methods, 
and then we will find that such tables are in confiict with quantum mechanics. 
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3.1.1 Probability distributions 

The most common, and historically relevant, method of LHV table construction 
is by using probability distributions. By looking at many instances of an LHV 
table, an instance meaning a particular set of values for the table entries, we can 
generate statistics for measurement results. Then we can speak of the probability 
that a particular set of outcomes is realized for a given set of measurements. These 
probabilities can be written as P{Ai — > 01,^2 — > 02, . . . ; i?i — > &i, -B2 ^ ^2, • • •) = 
P{{Aj —>■ ttj}; {Bk -^ &fc}) for two systems A and B and where Aj and B^ are the 
measurements and aj,bk G { — 1,1} are the outcomes. These probabilities must be 
positive and normalized, 

P{{Aj -. a,}; {Bk -^ bk}) > 0, 

^ ^ (3 1) 

and J2 E P{{A,^a,};{B,-.b,}) = l. 

{a,}=-l{bfc}=-l 

Given this probability distribution we know that the LHV table assigning values of 
aj to Aj and bk to Bk occurs with probability P{{Aj ^ aj}; {Bk ^ bk})- 

3.1.2 Random variables 

An alternative formulation of LHV tables that we give a bit more attention is their 
construction through random variables. Much work has been done analyzing local 
realism through probability distributions [10, 13], while this alternative perspective 
using random variables has been largely ignored. However, this is the perspective 
for much of our work. 

Random variables are used to give measurement outcomes, and since we only 
consider measurement outcomes ±1, we only consider two valued random variables 
that produce values of ±1. 
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Definition 2. A random variable, R{a), is such that 

+1 with probability (1 + a)/2 
— 1 with probability (1 — a)/2 




Note that the parameter a is the expectation value of the random variable since 
{R{a)) = ((l+a)/2)(+l) + ((l-a)/2)(-l) = a. To get a feeling for this definition of a 
random variable, let us consider a few cases. One special case is the random variable 
i?(±l) which takes the value ±1 with unit probability. We take the product of two 
random variables to mean a random variable whose outcome is the product of the 
outcomes for the two random variables being multiplied. With this understanding, 
we have i?(±l)i?(a) = ±R{a) = R{±a). Also another case of interest is when a = 0. 
In this case -R(O) is a 50-50 random variable, like an unbiased coin fiip. 

The way LHV tables are constructed using random variables is to use an in- 
stance of R{{Aj)) to determine the value of aj, and likewise for Bk- This guar- 
antees that the average value of aj agrees with the measurement statistics since 
(R{{Aj))\ = (Aj). What is left is to make sure the individual random variables 
are correlated so as to reproduce all the correct measurement correlations, meaning 
that (R({Aj))R({Bk))) = (AjBk). If there are more than two systems, then we 
must find correlated random variables that also reproduce correlations involving the 
other systems. Once the proper random variables have been found, instances of LHV 
tables are produced through instances of the random variables. 

It may or may not be obvious that random variables and probability distributions 
are equivalent methods of generating LHV tables. If this is not obvious, simply 
note that not only can probability distributions and random variables generate LHV 
tables, but also a probability distribution or a set of random variables could be 
constructed, in the asymptotic limit, given access to instances of an LHV table. 
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Thus, these two approaches must be equivalent. 

3.1.3 Allowed correlations for random variables 

Now we back away from LHV tables and consider random variables in full abstract- 
ness. Our task is to figure out the possible correlations that can be found among 
random variables. We do this in two steps. First, we consider only two random vari- 
ables and learn what correlations can exist between them. Then we can generalize 
this result to n random variables inductively. In the next section, with this task 
completed, we use random variables to construct LHV tables. 

A simple calculation yields the product of two uncorrelated random variables. 
R{a) and R{b) multiply to +1 when R{a) and R{b) are both +1 or both —1. The 
probability that this happens is 

l + a\ fl + b\ / l-a \ f l-b \ _ 1 + ab 



2 J \ 2 J \ 2 J \ 2 ^ 
which means that R{a)R{b) = R{ab). 

For correlated random variables, R{a)R{b) = R{p) and we wish to find the pos- 
sible values of p = (R{a)R{b)), the average value of the product of the random 
variables. To do this, consider Fig. 3.2(a). The bars in the figure have length 1 and 
if you pick a point in the top bar, you will pick it on the left part with probability 
(1 — a)/2 and on the right part with probability (1 + a)/2. Therefore randomly 
picking a point on the line is equivalent to an instance of the random variable with 
the left and right sides corresponding to —1 and +1 outcomes, respectively. 

First let's see how correlated our variables can be. We want to arrange things so 
that we pick both random variables to have the same value as often as we can. This 
is illustrated in Fig. 3.2(b). Clearly every time R{a) = —1, that is (1 — a)/2 of the 
time, we can also pick R{b) = — 1 and every time R{b) = +1, which is (1 + 6)/2 of 
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l-a 



1 + a 



\ fl(o) = -l 



R(a) = +1 



1-6 



1 + 6 



R{b) = -1 



R(b) = +1 



(b) 



(c) 



l-a 



1 + 6 



I R(ajR(b) = +1 I R{n)R(b) = -1 



R{a)R{b) = +1 



1-6 



l-a 



a + b 



R{a)R(b) = -1 



R{ajR(b) = -1 I R{a)R{b) = +1 1 



Figure 3.2: Visual aids useful for finding the possible values of two correlated random 
variables. The figure in (b) represents the case of maximally correlated random 
variables, while (c) represents the case of maximal anticorrelation. 



the time, we can simultaneously pick R{a) = +1. After this, though, we are stuck 
having R{a) and R{b) with opposite signs at least (a — b)/2 of the time. In this case 
the probability that the product of R{a) and R{b) is +1 is 



1 



+ 



1 + 6 



l + (l-(a-6)) 



(3.3) 



2 y V 2 

which means that R{a)R{b) = R{1 — [a — b)) = R{1 — \a — b\). Note that this 
argument only requires that a < b. However, since reversing a and b gives the same 
formula, it must be valid for all a and b as well. Random variables R{a) and R{b) 
satisfying R{a)R{b) = R{1 — \a — b\) are called maximally correlated. 

Now let's see how anticorrelated our random variables can be. In the figure 
(1 — a)/2 < (1 + b)/2, so we can pick R{a) to be negative every time R{b) is positive. 
Likewise we can pick R{b) to be negative every time R{a) is positive. As is illustrated 
in Fig. 3.2(c), this means the probability of getting a negative result is 

b\ l-((a + 6)-l) 



l-a\ /I 
+ 



(3.4) 
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which means that R{a)R{b) = R{{a + 6) — 1) = R{\a + b\ — 1). The assumption here 
was that a > —b. But if a < —b, then —a > —{—b) and the above result, being 
insensitive to replacing a by —a and b by —6, gives us the maximal anticorrelation for 
R{—a)R{—b) which is the same as for R{a)R{b). Random variables R{a) and R{b) 
satisfying R{a)R{b) = R{\a + b\ — 1) are called maximally anticorrelated. Combining 
these two results and noting that these maximal correlations can always be achieved 
gives us the following result^. 

Lemma 1. Let a,b & [—1, 1]. Then R{a)R{b) = R{p) is possible if and only if 

\a + b\-l <p<l~\a-b\. 

Let's test out this result on a few special cases. Suppose 6 = ±1. Then Lemma 1 
gives us that p = ±a as expected. Also, we see that p can take any value in [—1, 1] if 
and only if a = 6 = 0, which is also not surprising. A special case that will permit a 
generalization to n variables is when the random variables are more likely to be +1 
than —1. For this case we get the following corollary to Lemma 1. 

Corollary 1. Let a,b E [0, 1]. Then R{a)R{b) = R{p) is possible if and only if 

a + b-l<p<l-\a-b\. 

This corollary is the one we will be able to generalize to n random variables. 
When we do so, we get the following theorem. 

Theorem 1. Let Oi > 02 > ■ ■ ■ > a„ G [0,1]. Then i?(ai) ■ ■ ■ i?(a„) = R{p) is 
possible if and only if p E [—1, 1] and 

(ai H ^an) - {n-1) <p < {n -1) - {ai-\ h a„_i - a„). 



^In Ref. [ ] it was shown that R{a)R{b) = R{p) ^ |a + 6| - 1 < p < 1 - |a - 6| for 
some nonlocal-hidden-variable theories. 
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The proof of Theorem 1 is a bit lengthy, and it is not terribly enlightening, so 
it is given in Appendix D. We can now use this result to prove the corresponding 
generalization of Lemma 1. 

Corollary 2. Let \ai\ > \a2\ > ■■■ > |a„|, with ai,a2,--- ,an G [—1,1]. Then 
R{ai) ■ ■ ■ R{an) = R{p) is possible if and only if p & [—1, 1] and 

(|ai| H h \an\) - (n-l) < Si---SnP <{n-l) - (|ai| H h |an-i| - knl), 

where Sj = aj/\aj\ is the sign of Qj. 

Proof. R{ai) ■ ■ ■ R{an) = R{p) -^ -R(|ai|) ■ ■ ■ -R(|an|) = -R(si ■ ■ ■ s„p) and the result 
follows from Theorem 1. D 

What we do now is consider for a moment the content in Theorem 1. First, notice 
that if any aj is then the interval is symmetric; this is why we do not need to define 
the sign of aj = 0. Also, recall that for two random variables, both need to have 
expectation value in order for the average of their product to be anything in the full 
interval, [—1,1]. So, suppose that all random variables have the same expectation 
value a > 0. In this case, the upper bound in Theorem 1 is always greater than or 
equal to 1. For the lower bound to be —1 we need a < {n — 2)/n. For n = 2, a of 
course needs to be 0, but we now see that a ^ 1 in the limit that n -^ oo. That is, 
for very large n, it is possible to find random variables all with expectation values 
close to 1 and yet with arbitrary correlations between them. With a bit of thought, 
though, it is not a surprising result that increasing the number of random variables 
also generally increases the possible correlations among them. 

Also note that proving this theorem for a particular value of n automatically 
proves it for all smaller n. For, inserting ai = 1 into the result of Theorem 1 yields 
the same formula for one less random variable. 
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3.2 Bell inequalities 

We return now to considering a local realistic description for measurement results 
on quantum systems. It turns out that we can develop an experimentally testable 
quantity capable of revealing nature's nonconformity to local realism. The idea is 
to construct a Hermitian observable, B, called a Bell operator, whose expectation 
value under LHV tables is bounded by some value, (i3)LR, the maximal local realistic 
value, yet quantum mechanics predicts a possible expectation value exceeding this 
bound [17]. Thus given a Bell operator, we have the Bell inequality 

(B) < {B)lr. (3.5) 

A predicted violation by quantum mechanics of any Bell inequality would prove that 
quantum mechanics is not both local and realistic. 

Developing experimentally useful Bell inequalities, or even Bell inequalities at 
all, has been the subject of an immense body of research, not a fraction of which we 
could hope to cover here. As just a taste of what has been done see Refs. [2, 16, 22, 
26, 47, 53, 59, 63, 66, 72, 76, 99, 100]. Instead, we will review the most important 
early work on Bell inequalities before discussing some novel results and approaches. 

3.2.1 CHSH inequality 

The most well known, useful, and historically relevant Bell inequality was invented by 
Clauser, Home, Shimony, and Holt and is therefore called the CHSH inequality [26]. 
This Bell inequality emerges from the following thought experiment. Suppose that 
two qubits are allowed to interact for some time and then they are sent in opposite 
directions, one to Alice and the other to Bob. Alice performs one of two projective 
measurements, Ai or A2, while Bob measures either Bi or i?2- This hypothetical 
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experiment is performed many times so that statistics can be gathered on the mea- 
surement results, but with the chosen measurements being varied each time. With 
these statistics, one can calculate the average value of the CHSH Bell operator, 

-BcHSH = Ai0Bi + Ai0B2 + A20Bi-A20 B2. (3.6) 

If the results of the experiment are describable by a local realistic model, then 
we can bound the value of (-Bchsh)- Each time this experiment is run, the table 
expressing the measurement results has 0^,6^ = ±1, and in this run we can write the 
value of i3cHSH as aihi + 0162 + 02^1 — ^262 = ai(6i + 62) + a2(&i — ^2) which has the 
value ±2ai or ±2a2 depending on whether hi = 62 or 61 = —62- Thus when averaging 
over all the trials we have that {Bcusu)lr < 2, giving the CHSH inequality, 

(-Bchsh) = (^1 (S) Bi + Ai (g) B2 + A2 (S) Bi - A2 (S) B2) < 2. (3.7) 

We now demonstrate how quantum mechanics defies local realism by finding that 
this inequality can be violated. To see this, write the projective measurements, with 
assumed eigenvalues ±1, of Alice and Bob as Aj = Sj ■ a and Bk = bk ■ <?, where Sj 
and bk are unit vectors. This gives us a quantum mechanical upper bound. 



{Ai ®Bi + Ai®B2 + A2®Bi-A2® B2) 

= (fli ■ a ® {bi + b2) ■ a + a2 ■ (y ® (bi - 62) ■ f?) 

= (ai ■ (T (g) 6+ • a)\bi + 62I + (^2 • <? ® &- ■ <?)|&i - h] 

< 1^1 + ^2! + 1^1-^21 <2V2, 



(3.8) 



where b± is a unit vector along 61 ±62- The second to last inequality comes from the 
fact that 6-t- ■ (T is an observable with outcomes ±1. The last inequality comes from 
considering the angle 6 between 61 and 62, 

\bi + hi + \bi-b2\ = V2 ("Vl + cos^ + Vl-cos^) , (3.9) 

and maximizing over 6. 
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This quantum mechanical upper bound can furthermore be achieved by setting 
Ai =X, A2 = Z, Bi = (1/V2)(X + Z), and B2 = (1/V2)(X - Z), and considering 
the state \ip) = (l/v^)(|00) + |11)). Since 2v^ > 2, we then have that quantum 
mechanics makes predictions beyond the capabihties of any local realistic theory. 



3.2.2 Linear programming approach 

Eq. (3.7) is a necessary condition for local realism to hold, meaning that its violation 
by a theory disqualifies it from being local and realistic, but that does not mean 
it is sufficient, meaning that its satisfaction guarantees that the theory violates no 
Bell inequality. The issue of whether the CHSH inequality is sufficient has been 
considered by others, with current progress suggesting that it is not [29, 45, 46, 94]. 
This leads to the question of whether one could find a set of Bell inequalities that, 
when taken together, are sufficient for a local realistic description of measurement 
outcomes to exist. This prospect has been studied by Peres [ ] and he has indeed 
found such a complete set. He also found that determining all the Bell inequalities 
is an NP hard problem. 

Peres' method for constructing a sufficient set of Bell inequalities answers the fun- 
damental question of whether such a finite set exists and allows one to generate Bell 
inequalities for a given number of measurement settings. In a typical Bell inequality 
experiment, one of these Bell inequalities is chosen and then the experimental setup 
is designed to demonstrate a violation of that Bell inequality. In this subsection, we 
outline a much simpler approach to generating Bell inequalities than the one taken 
by Peres. Our approach is also particularly relevant to experimental violations of 
local realism because, given a robust experimental setup, our approach generates 
Bell inequalities that should be violated by the experimental results. 

Suppose, for simplicity, that we have two separated systems on which measure- 
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ments Aj, 1 < j < rriA, and Bk, 1 < k < rriB, are made. The generalization to 
more than two systems is straightforward. If we suppose a local realistic descrip- 
tion of the measurement outcomes exists, then there must be, as per Subsec. 3.1.1, 
a column vector P, which if there are two outcomes per measurement has dimen- 
sion 2™-*+™^, whose entries are the probabilities that a certain table is realized, 
P{{Aj —>■ aj}; {Bj -^ bj}). For instance, in the case as before of two measurements 
we have, 



(3.10) 



Note that this permits a generalization to measurements with any number of out- 
comes if desirable. 

— * 

Any correlation can be written as a linear combination of entries in P. Let us 
write the correlations that can be observed in some experiment, of which there will 
be {rriA + l){mB + 1) if we include the / measurement, as entries in a vector, C. Then 
this implies that a local realistic description of these correlations exists if and only 
if, for the appropriate matrix M, there exists a vector P with non-negative entries 
such that MP = C . In the case that measurements have two outcomes ±1, M will 
be an (ttia + 1)(^^b + 1) x 2^^^^^ matrix whose entries are ±1. In general, M is 
completely determined by the number of measurements and the number of outcomes, 
as well as their allowed values. Thus the problem is, given M and C, to solve, 

MP = C subject to P > 0. (3.11) 

Notice we do not include the normalization constraint separately because it is con- 
tained in the equation MP = C. In fact, we assume that the first row of M has all I's 
and that the first entry in C is 1, as this takes care of normalization. We also adopt 
the notation that A > B when every entry of A is bigger than the corresponding 
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entry of B. Problems like these have been studied extensively and are called linear 
programs . 

The most general form of a linear program is, given L, M, and C, to solve 

minimize L'^P subject to MP >C and P > 0. (3.12) 

Just about any linear optimization problem involving linear constraints can be put 
into this form. We will see later how our problem can be put into this form, but first 
there is an important feature of linear programming needing discussion. For every 
linear program, Eq. (3.12), there is a related dual linear program. 

maximize C'^Q subject to M^Q < L and Q>0. (3.13) 

Notice that P appears nowhere in the dual program. This program seeks a completely 
different quantity Q that even generally has a different dimension than P. Yet there 
is an important relationship between these problems. The strong duality theorem for 
linear programming states that the minimum value of L^P in Eq. (3.12) is exactly 
equal to the maximum value of C^Q in Eq. (3.13) [21]. We will see that duality 
in linear programming reveals a duality between probability distributions and Bell 
inequalities. 

Returning to Eq. (3.11), we try to rewrite our problem as a linear program, 
minimize 1^{MP-C) subject to MP>C and P > 0, (3.14) 

where 1 is a vector of all I's of appropriate dimension. A solution to Eq. (3.11) exists 
if and only if this minimum is 0. The reason is that MP > C assures that every 
entry of MP — C is positive. Thus the sum of all the entries in MP — C is if and 
only if each entry of MP — C is zero, i.e MP = C. Now we easily turn this into a 
linear program. 

minimize (M^l )^P subject to MP>C and P > 0. (3.15) 
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If the minimum value is l^C, then there exists a probabihty distribution, P, pro- 
ducing the given measurement correlations, C. Otherwise, no such probability dis- 
tribution exists. 

As a practical matter, this linear program gives us a method for determining 
whether a set of correlations violates local realism. In terms of theory, it gives a 
deep connection between probability distributions and Bell inequalities as we show 
below. The dual program of Eq. (3.15), which eventually produces Bell inequalities, 
is 

maximize C^Q subject to M^Q < M^l and Q > 0. (3.16) 

From the strong duality theorem we know that there exists a local realistic description 
of the correlations in C if and only if this maximum value is I'^C. Since Q = 1 
satisfies the constraints and achieves this value, we know that generally the maximum 
value must exceed I'^C. 

From this we deduce that there exists a local realistic description of the correla- 
tions if and only if 

maximize C^{Q-1) subject to M^{Q -1)<0 and (Q -!)>-! (3.17) 

is 0. This gives us the set of Bell inequalities C^{Q — 1) < for all vectors Q — 1 
satisfying the constraints. This is a complete set of Bell inequalities in that they are 
both necessary and sufficient for a local realistic description to exist. They all have 
the form of a linear combination of correlations being less than 0. 

We can simplify these Bell inequalities in a few ways. The first thing is to 
recognize that the constraint {Q — \) > —1 is not needed. For, suppose that we have 
a vector Q — 1 satisfying M'^{Q — 1) < 0, but such that {Q — 1) > —ql for some 
g > 1. Then the vector [Q — l)/g satisfies all constraints and thus we get the Bell 
inequality C^{Q — l)/^' < -v^ C^{Q — 1) < 0. So in fact we get a Bell inequality 
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for all vectors Q — 1 satisfying M'^{Q — 1) < 0, regardless of whether they satisfy 
the second constraint. 

A second simplification can also be used to deal with the constraint M'^{Q — 1) < 
0. Since we assume the first row of M is all I's, due to P needing to be normalized, we 
can break upQ — lasQ — 1 = ( —go Q ) > where g has one less dimension than Q. 
If we also break up (7 = ( 1 c ) , then the Bell inequalities are —qo+c^q < for all 
— gol + M'^q < 0, where M is M without the first row. To have a complete set of Bell 
inequalities it suffices to pick the smallest go and so we have go = (M-^g)max, where 
max denotes the maximum entry of the vector, giving us all the Bell inequalities, 

c^g-'<(M^g-)^ax for all q. (3.18) 

In some sense this is an obvious result. It states that all the Bell inequalities are 
obtained by considering an arbitrary linear combination of correlations, and then 
saying that this must be less than the maximum possible value of those correlations 
under extremal LHV tables. However, this result is very important because it says 
we need only consider linear combinations of correlations, a fact that is far from 
obvious to begin with. It is also important because if a local realistic description of 
a set of correlations does not exist, it provides an algorithm that generates a Bell 
inequality which the correlations violate. 

What we have shown here is a computational method that determines whether 
a set of correlations can be described within a local realistic framework. If the set 
is so describable, then the linear program gives a probability distribution resulting 
in those correlations. If the set is not describable by a local realistic model, then 
the linear program is capable of generating a Bell inequality that the correlations 
violate. Practically, this is a useful thing to be able to do. In terms of theory, 
linear programming gives a way to generate all the Bell inequalities and through the 
principle of duality in linear programming we see that probability distributions and 
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ABC 

• • • 

Ai Bi Ci 

M B2 C2 



■Am A Bms Cmc 

Figure 3.3: To derive Bell inequalities from random variables we assume the following 
setting. We allow an arbitrary number of parties, listed here as party A, B, C, etc. 
Each party controls a subsystem of the total system, denoted as filled in circles. Each 
party can furthermore make an arbitrary number of measurements, although these 
measurements must give only two outcomes, ±1. This last restriction is the only 
lack of generality in the setup. 



Bell inequalities are also dual to each other. 



3.2.3 Random-variable Bell inequalities 

It is typically computationally difficult to solve the linear program in Subsec. 3.2.2 
because the number of columns in the matrix M grows exponentially with the number 
of measurements. Therefore, it is useful to use a different approach to derive explicit 
Bell inequalities, although such a set may not be complete. In order to come up 
with Bell inequalities for a wide range of configurations, we take a novel approach 
involving random variables. We derive Bell inequalities using random variables for 
the situation illustrated in Fig. 3.3, which is indeed quite general. 

Before considering the most general case, let us consider the specific example of 
two parties with two measurements each. The types of correlations that exist are then 
of the form {AjB^), giving us a list of four correlations that must be simultaneously 
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satisfied in order for a local realistic description of the correlations to exist. 

i?((A))i?((5i)) = i?((A5i)) 

R{{A^))R{{B,)) = R{{A,B,)) ^^ ^^^ 

R{{A2))R{{B^)) = R{{A2B^)) 

Ri{A2))Ri{B,)) = RiiA^B,)) 

A necessary condition for these four equations to have a solution is that their product 
have a solution. Thus multiplying these four conditions together gives us the consis- 
tency condition, 1 = R{{AiBi))R{{AiB2))R{{A2Bi))R{{A2B2)). This condition can 
be turned into an algebraic quantity by using Corollary 2 in Subsec. 3.1.3, giving us 
the Bell inequalities 

\{AiBi)\ + \{A,B2)\ + \{A2Bi)\ + 1(^252)1 - 3 < SUS12S21S22 

(3.20) 

and siiSi2S2iS22 < 3 - \{A,B,)\ - K^i^a)! - \{A2Bi)\ + K^a^a)], 

where Sjk is the sign of {AjB^). If S11S12S21S22 = 1 then the first inequality is always 
satisfied, while if S11S12S21S22 = —1 the second inequality is always satisfied. Because 
of this, we can combine both above inequalities into a single Bell inequality, 

\{AiBi)\ + \{AiB2)\ + 1(^2^1)1 - snSuS2iS22\{A2B2)\ < 2. (3.21) 

After a little investigation, it becomes clear that this is equivalent to the familiar 
CHSH inequality. 

The general case is a straightforward generalization of the previous example. 
The types of correlations that we must capture have the form {A'jB^Ci ■ ■ ■), where 
a,b,c,... G {0, 1} determine whether a measurement is included in the correlation. 
This means the random variables that produce the measurement outcomes must be 
correlated so that /[R{{Aj))]^[R{{Bk))]''[R{{Ci))Y---\ = {A'^ B^Cj; ■ ■ ■) . This gives 
us a list of correlations that must be simultaneously satisfied in order for a local 
realistic description to exist. Multiplying together subsets of these equations and 
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using Corollary 2 gives us that 

n n n— 1 

Y, m,)\ -in-l)<l[s,<{n~l)- 5^ |(M,)| + |(M„)|, (3.22) 

where the Mj, of which there are n, are measurements in the subset that do not cancel 
and Sj is the sign of (Mj). One thing we should point out is that each individual 
system measurement must appear an even number of times throughout all the Mj, 
which implies that all the Mj must multiply to ±/. This is a constraint that is 
somewhat difficult to understand unless one writes down a few examples along the 
lines of our random variable derivation of the CHSH inequality. Now, as before, one 
of the inequalities in Eq. (3.22) is always satisfied. This allows us to write the most 
general form of a random- variable Bell inequality as 

71—1 / " \ 

Y,\{M,)\-il[sA\{M^)\<n-2. (3.23) 

3=1 \i=i / 

According to Theorem 1 we should choose \{Mj)\ < |(M„)| for all j, but in fact 
Eq. (3.23) needs to hold for any ordering in order for a local realistic model to exist, 
so we define random- variable Bell inequalities to include any choice for M„. We 
can easily write down this Bell inequality after picking a set of correlations, (Mj), 
although this set of inequalities is typically not complete. For an example of a Bell 
inequality inequivalent to any in this set, see Ref. [59]. 

3.2.4 Two qubits 

In this subsection we analyze random-variable Bell inequalities, Eq. (3.23), in the 
special case of two parties, but an arbitrary number of measurements. In this case 
we prove that all random- variable Bell inequalities are equivalent to the CHSH in- 
equality. That is, if the CHSH inequality holds for all choices of measurement then 
so do all random-variable Bell inequalities on two qubits. Conversely, since the 
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CHSH inequality is a special case of a random-variable Bell inequality, satisfaction 
of Eq. (3.23) guarantees that the CHSH inequality holds. 

If we consider the value of n in Eq. (3.23), which is the number of joint mea- 
surements, we can phrase the problem as follows. For n = 1, Eq. (3.23) says that 
- 1 (/(g)/) I < -1 which always holds. For n = 2, we have that \{A®B)\-\{A®B)\ < 
which also always holds. Continuing with n = 3, we get that \{A ® I)\ + \{I ® B)\ — 
\{A^ B)\ < 1, and sign variations thereof. Because of how random- variable Bell in- 
equalities are constructed, we know that / must appear as indicated so that the three 
operators can multiply to the identity. This inequality must also be satisfied because 
it is equivalent to R{{A (8 /))/?((/ (8 B)) = R{{A (8 B)). One can find such random 
variables as follows. Prepare the quantum system in a lab and make measurements 
oi A ^ B. The outcomes of A and B determine R{{A ® /)) and /?((/ ® B)), and 
these multiply to R{{A^B)). The first instance of Eq. (3.23) that is not necessarily 
satisfied by quantum states is n = 4, which gives rise to CHSH type inequalities. 
Our claim is that the n = 4 case being satisfied is sufficient for all values of n to be 
satisfied. 

Lemma 2. For two qubits, all random-variable Bell inequalities are satisfied if and 
only if the special case of n = 4 is satisfied. 

Proof. The forward direction is trivial. To prove the converse, suppose that 

\{Ai ® fii)| + \{Ai ® ^2)1 + 1(^2 ® Bi)\ - SuSl2S2lS22\{A2 ® ^2)1 < 2. (3.24) 

To show that Eq. (3.23) must now be satisfied we consider the possible forms of M„. 

The first case we consider is M„ = / ® Bi, with Bi ^ I since we can always 
reorder the Mj to have M„ ^ I ® I provided a non-identity measurement exists. 
Since Bi must occur in some other Mj, we choose it to appear in M„_i = Ai® Bi. 
Now we add and subtract (M) = {Ai ® I) from Eq. (3.23) and what we want to 
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prove becomes 



n-2 



V I (M,-) I - (Si ■ ■ ■ Sn-2S) I (Ai ® /) I + (Si ■ ■ ■ Sn-2S) \ {A^ ® /) | 

j=i (3.25) 

+ |(Ai ® 5i)| - (si ■ ■ ■ s„_2s)(s„_is„s)|(/ ® fii)| < n - 2, 

where s is the sign of (M) . Now if we can prove that 

n-2 

J] I (M,-) I - (si ■ ■ ■ s„_2s) I (Ai ® /) I < n - 3, (3.26) 

then combining this with the fact that (si ■ ■ ■ s„_2s)|(yli ® /)| + |(Ai ® i?i)| — 
(si ■ ■ • s„_2'S)(s„-is„s)|(/ ® -Bi)| < 1, which follows after some thought from the 
n = 3 case, we will have proven Eq. (3.23). Notice that Eq. (3.26) is an instance of a 
random- variable Bell inequality with one fewer measurement. This follows from the 
fact that we eliminated an even number of individual measurements, thereby leaving 
the equation with an even number of individual measurements. Thus we have used 
the n = 3 case to reduce the value of n by one. 

This strategy of reducing the value of n works with the other possible forms of 
Mn as well. If M„ = Ai ^ I, then we can use the same strategy as before to reduce 
the value of n. The only remaining case is if M„ = Ai ^ Bi, with Ai, Bi ^ I. In 
this case we choose (M„_i) = {A2® Bi), which we know must be a term in the sum 
because Bi must appear again. We then choose M„_2 = Ai (^ B2, which also must 
exist because Ai must appear again but cannot appear in M„_i. Now we add and 
subtract M = A2 ® ^2 from Eq. (3.23) to get 



71—3 

J2 I (Mj) I - (Si ■ ■ ■ S„„3S) I {A2 ® ^2) I + (Si ■ ■ • S„„3S) I {A2 ® B2) I 

+ \{A,^B2)\ + \{A2(S)B,)\ 

- {Si- ■ ■ Sn-3s){Sn-2Sn-lSnS)\{Ai ® 5i)| < U - 2. 
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As before, this equation follows provided we can prove that 



n— 3 



J2\{Mj)\~{si--- s„_3s) \{A2®B2)\<n- 4, (3.28) 



and 



(Si ■ ■ ■ Sn-3S)\{A2 ® ^2)1 + 1(^1® ^a)] + 1(^2® Bi)\ 
-(Si ■ ■ ■ Sn-3s){Sn-2Sn-lSnS)\{Ai ^ Bi)\ < 2. 



(3.29) 



This last equation, again with some thought, follows from the n = 4 case. The 
first equation, Eq. (3.28), is simply a random- variable Bell inequality with two fewer 
measurements because, as before, we removed an even number of individual mea- 
surements. 

Thus we can always reduce the problem to proving a random-variable Bell in- 
equality involving either one or two fewer measurements. In this way we reduce the 
problem, by step sizes of one or two, to a random-variable Bell inequality involv- 
ing either three or four measurements. We conclude that all random-variable Bell 
inequalities follow from the n = 4 case. D 

For more than two qubits a similar result likely holds, but the details for more 
parties has not yet been worked out. The conjecture would be that for m qubits 
one only need consider 2™ joint measurements or fewer, with all larger number of 
measurements following as a consequence of these cases. 



3.3 Generalized probabilistic theories 

So far we have discussed the fact that measurement correlations in quantum me- 
chanics cannot, in general, be described in local realistic terms. In order to better 
understand this violation of local realism, it helps to look at quantum mechanics in 
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terms of a broader framework. By comparing quantum mechanics to other theories, 
many of which may also make predictions contrary to local realism, we can better 
understand the nature of quantum correlations and pose questions such as why do 
quantum correlations not violate local realism more. 

Just to motivate why one should consider quantum mechanics as part of a larger 
framework, let us consider, once again, the CHSH Bell operator Bqusu = ^i-Bi + 
A1B2 + A2B1 — A2B2- Since we have that —1 < {AjB^) < 1, the naive bound one 
would place on (Schsh) would be (Schsh) < 4. This is the fundamental bound on 
the expectation value of the CHSH Bell operator. We found that (;Bchsh)lr = 2 
and (i3cHSH)QM = 2-\/2, where {B)qm is the maximum value of {B) over all quantum 
states [25, 97]. So neither local realism nor quantum mechanics achieves the CHSH 
fundamental limit. 

Popescu and Rohrlich asked whether the fact that quantum mechanics does not 
achieve the CHSH fundamental limit is a consequence of a principle called no signal- 
ing [ " ] . The no-signaling principle states that a party, Alice, cannot do anything to 
her system that would affect the distribution of measurement outcomes on a spatially 
separated system controlled by Bob. Note that this is different than locality. Local- 
ity assumes that Alice cannot affect Bob's system in any way, while no signaling only 
states that Bob cannot measure the effects produced by Alice. In a realistic theory, 
the two concepts coincide. In the end, Popescu and Rohrlich found that there could 
indeed exist a no-signaling theory in which the fundamental limit of 4 for the CHSH 
Bell operator is achieved. To prove this, they constructed a fictitious device called a 
PR box, named after its inventors, which accomplishes this feat. 

Thus the question arises of why quantum mechanics does not violate the CHSH 
inequality more than it does [15, 20, 21, 61, 93]. While there has been some success 
arguing that no physical theory should be able to achieve (iScHsn) = 4, no one has 
yet been able explain why quantum mechanics stops at 2-\/2. We will review the idea 
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Alice 



Bob 



Aj 



Gj ^ 




(-l)(^-i)^^-i)=a,6 



Figure 3.4: A visualization of a hypothetical PR box. The device has two ends, one 
at the location of Alice and the other at Bob. At her end, Alice can input what 
measurement she wants to make, and likewise for Bob. After Alice's measurement of 
Aj, with j = 1, 2, is input into the box, she immediately gets an output aj G { — 1,1}. 
This measurement result, from the perspective of Alice, is completely random. Bob 
gets an output as well, and when Alice and Bob meet later to compare results, they 
find that their results are perfectly anticorrelated if j = fc = 2 and perfectly correlated 
if either j = 1 or A; = 1. The box did this by giving one party a random result, and 
then giving the other party either the same result or its negative depending on Alice 
and Bob's measurement choices. 



of a PR box, and also consider a framework for considering all probabilistic theories. 
Doing this will enable a comparison between the correlations in quantum mechanics 
and those in other no-signaling theories. Perhaps, in the future, this will lead to a 
fundamental understanding of quantum correlations. 

3.3.1 PR box 

In this subsection we give an example of a theory, other than quantum mechanics, 
that makes predictions violating local realism. In this theory there exists a device 
called a PR box [12, 77] which cannot be used to signal instantaneously between 
Alice and Bob, but which gives Alice and Bob access to correlations that achieve 
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{BcHsu) = 4. 

The way a PR box works is explained in Fig. 3.4. The outputs of a PR box 
are perfectly correlated for AiBi, A1B2, and A2B1, and perfectly anticorrelated for 
A2B2. This gives (-Bchsh) = 1 + 1 + 1 — (—1) = 4. This device furthermore cannot be 
used to signal between Alice and Bob, because regardless of one party's measurement 
choice the other party simply gets random measurement results. Thus this device 
preserves no signaling yet achieves the CHSH fundamental limit. Obviously no such 
device, then, could be constructed in the real world. 



3.3.2 General framework for a single system 

We now discuss a general framework in which local realistic theories, quantum me- 
chanics, and PR boxes can all be discussed and compared [7, 8, 10, 62]. A basic 
property of these theories is that they all accept the existence of systems on which 
operations can be performed and measurements can be made. Measurements then 
have outcomes that occur with some probability, with determinism being recovered 
when probabilities are 1 and 0. Theories in this framework are therefore called gen- 
eralized probabilistic theories. Our first order of business is to discuss single systems 
in generalized probabilistic theories. Once we do this we can then understand how 
to deal with multiple systems, which are needed in order to discuss correlations and 
local realism. 

If we have many independent copies of a system, then by making a measurement 
on those copies we can determine the probabilities of the measurement outcomes. 
In general we need many copies because a measurement might disturb the system. 
We can repeat this process until we have done all measurements, or just until we 
have done enough measurements to be able to predict measurement probabilities for 
all measurements. For example, in quantum mechanics, we know the probability 
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for all measurement results given that we know {X), (Y), and (Z), so we can stop 
making measurements once we know the probabilities for the Pauli measurements. 
Then for any measurement A = a ■ a, (A) = {a ■ a) = Xlfci 0'j{(^j)- Thus for any 
generalized probabilistic theory we assume that there exists some number of fiducial 
measurements, Fj, whose outcome probabilities determine the outcome probabilities 
for all measurements. 

Once we know the set of fiducial measurements for a system, of which we assume 
there are M, we can then completely characterize the system by the fiducial measure- 
ment probabilities. The state of a system in a generalized probabilistic theory then is 
described by a vector, P, whose entries are the probability of a fiducial measurement 
outcome, 



P 



( P{Fi - 
P{F,- 

P{Fm 
\ P{Fm 



1) \ 
-1) 



(3.30) 



-1)/ 



Note that this vector is different than the vector P used to generate LHV tables. In 
particular, we now have that the entries of P sum to M rather than 1. We assumed 
in Eq. (3.30) that there are two measurement outcomes, ±1. The state vector P 
is easily generalized to more than two outcomes, but the two outcome scenario is 
sufficient for our purposes. In the case of a qubit where M = 3 and the state is 
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characterized by the Bloch vector a, we have that 
f l + aA 



P 



1 — ai 

1-02 
1 + 03 



(3.31) 



Notice that 3 parameters are sufficient to describe the state of a qubit, while P 
has 6 entries. This redundancy conies from the fact that we have not yet enforced 
that the state be normahzed. Thus we require that a vahd state, Eq. (3.30), be 
normahzed in that the probabihties for each fiducial measurement sum to 1. This 
requirement can be written as 



P(F,^l) + P(F,-.-l) 



/ 



0---0 1 1 0---0 



\ 



P(F, 



-1) 



V 



njP, 



(3.32) 



for j = 1, 



M and where nj = ej 



1 1 



with Cj being an M dimensional 



vector with a 1 in the jth spot and zeroes elsewhere. Our normalization constraint 
is now conveniently written as n^ P = 1, and so we call the vectors fij normalization 
vectors. Normalization vectors are orthogonal with 



rij rik = Cj Ck ■ 



1 1 




26. 



'jk- 



(3.33) 



Since P is 2M dimensional, and since we have M orthogonal normalization vec- 
tors, we only need M more vectors and we will have a convenient basis for all 
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states. This role is filled by the mean vectors, fhj. These vectors are defined by 
rnj = e^ ® ( 1 — 1 ) , and their inner product with P gives the mean values of 
the fiducial measurements in that state, 

/ \ 



fhJP 



0---0 1 -1 0---0 



P(F, -. 1) 
P(F, -. -1) 



V 



(3.34) 



/ 



= P(F,^1)-P(F,^-1) = (F,). 
The mean vectors are orthogonal to each other and also to the normalization vectors, 



rhjrhk = efck 



rrij rik = e/ e^ 



-1 

1 
1 



2Sjk, and 



(3.35) 



■ 5ik = 0. 



So now we have 2M orthogonal vectors which serve as a basis for all states. This 
gives us that 



^ 1 *^ 



(3.36) 



i=i 



for some CLj,bj, and with an overall factor of (1/2) put in for convenience. There 
are still 6 parameters here, but by using the normalization condition, Eq. (3.32), we 
eliminate 3 of them, 



niP = l =^ 



1 ^ 

- "^{ajfilnj + hfilrhj) = 0^ = 1, 



(3.37) 



i=i 



where we used Eqs. (3.33) and (3.35) to evaluate the inner products. So now we 
actually have. 



^ 1 ^' 



(3.38) 



j=i 
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Finally, since the numbers bj represent the average value of the fiducial measure- 
ments, we have that —1 < bj < 1. The vector b whose entries are bj is analogous 
to the Bloch vector for qubits, since the entries of the Bloch vector are the average 
values of Pauli measurement outcomes. 

Any generalized probabilistic theory must have systems whose states, P, conform 
to Eq. (3.38) since the states are normalized and assign positive probabilities to 
fiducial measurement outcomes. This does not mean that any vector P in Eq. (3.38) 
with — 1 < 6j < 1 represents a state in a given theory. To specify a particular theory, 
one must specify the set of allowed states. States in this set must be represented by 
normalized vectors, P with positive entries, but the set need not contain all allowed 
states. The theory which does allow all states is called GNST, or the generalized 
no-signaling theory [']. Quantum mechanics, for instance, allows only vectors b 
whose norm is no greater than 1. This permits a representation of qubit states by 
vectors b lying inside a sphere. GNST, on the other hand, permits any vector whose 
entries are less than 1 in magnitude. These vectors lie not inside a sphere, but 
inside a hypercube. Just as quantum bits are called qubits, we call states in GNST 
generalized bits, or gbits. Gbits and qubits are contrasted in Fig. 3.5. 



3.3.3 General framework for multiple systems 

Now that we have a complete description of single systems in generalized probabilistic 
theories, we move to tackle the case of multiple systems. For obvious reasons we 
insist that no generalized probabilistic theory allow instantaneous signaling, and this 
is the fundamental principle on which we base our study of multiple systems. The no- 
signaling principle is implied by stating that for two systems, A which is controlled 
by Alice and B which is controlled by Bob, measurements and operations acting 
separately on A and B must commute. The reasoning is thus. Suppose all operators 
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Gbit 



Mi^2K 





Figure 3.5: Comparison of gbits, for M = 3, and qubits. On the left is a 2 x 2 x 2 
cube and on the right is a sphere of radius 1. Vectors inside the cube give rise to 
states in GNST, or gbits, and vectors inside sphere give rise to states of qubits. The 
sphere, if moved on top of the cube, would fit perfectly inside it. This means that 
states of qubits are also states of gbits, but the converse not true. The corners of 
the cube, for instance, do not correspond to any quantum state. 



acting separately on A and B commute. Then results of measurements on A are the 
same whether Bob does his operations before or after Alice. Since operations Bob 
does after the measurement by Alice cannot affect its result, neither can operations 
Bob performs before Alice's measurement have any effect. In what follows we assume 
only two systems for simplicity. The extension to more than two systems is fairly 
obvious. 

If we assume that operations acting separately on A and B commute, then the 
probability P{A ^ a; B ^ b) is unambiguously defined. Now suppose for the 
moment that systems A and B are completely uncorrelated with each other. In this 
case, P{A -^ a\B ^ h) = P{A — > a)P{B -^ b). Now we want to represent the 
state of these two systems with one state vector, using reasoning analogous to our 
motivation of the tensor product in Subsec. 1.3.4. Since the state of system A is 
determined by the fiducial measurements F^, and likewise for B, we conclude that 
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the state of the entire system is the tensor product of the individual system states, 



/ P{F^^1) \ 



P 



AB 



P^®P^ 



y p{F^, ^-i ) 



^ P{F^ -^ 1) ^ 
^ PM ^ -1 ^ 



^ P{F^ -^ l)PiF^ ^1) ^ 



y P(FA ^ -l)P(i^B 



M 



■n) 



(3.39) 



P(Pi^^l;Pf 



yp(F^^^-l;Fi,^-l) J 

where we assumed M fiducial measurements for both systems. Extending this to 
correlated systems leads us to the assumption that the state of multiple systems is 
a vector in the tensor product space of the individual system vector spaces. The 
overall state vector looks like the final equality in Eq. (3.39). 

Since we have a nice basis for a single system, involving normalization and mean 
vectors, we can just take a tensor product of these basis vectors to get a basis for 
multiple systems. Thus a state of two systems can be written as, 

, M M 

y^ X^(^i ® ^k + ajk'mj ®nk + bj^Uj m^ + Cjkrhj (g) m^), (3.40) 
j=i k=i 



^AB ^ ^^^^^ 



where we know, as before, that the coefficient oiuj^nk is 1 because of normalization. 
We will see later that aj^ and bj^ determine the average values of fiducial measure- 
ment results at A and B, and that the numbers Cjk express correlations between 
fiducial measurement outcomes. We find below the restrictions on ajk, bjk, and Cjk- 

To learn more about the states of multiple systems, let us consider the reduced 
states of Eq. (3.40). If we want to recover the state of Alice's system alone then 
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we need to know the probabilities of her fiducial measurement outcomes. We can 
determine P{F^ -^ aj) from P^^ by assuming that Bob makes a measurement of 
Ff and summing over the outcomes. Mathematically, 



3 
k "" -"--"-^ 

A , „ N D/rA , „ . pB , 1 \ I D^ rA , „ . rB 



P{F^ ^ a,) = P(F/ ^ a,; F^ -. 1) + P(i^^ ^ a,; F^ ^ -1). (3.41) 

This is accomplished via normalization vectors by the action (/ ® fi1)P^^, which 
needs to obtain the same result regardless of k. Thus we obtain the reduced state 
for Alice's system, 

P^ = {I®nl)P'^'', (3.42) 

and similarly for Bob. Using this result on Eq. (3.40) gives us that 



M 
5A 



P'' = \Y.^n,+a,km,), (3.43) 



and similarly for Bob. Since this result needs to be independent of fc, we get that 
ajk needs to be independent of /c, i.e., ajk ^ o,j. A similar argument also gives that 
bjk -^ bk- This gives us the state of two systems as, 

M M 

4 



P^^ = - ^ ^(^i ® ^A: + a.j^i ® ^fc + &fc% ® mfc + Cj^rhj ® rfik). (3.44) 



i=i fe=i 

Now we see that a^ = {F^) and bk = {F^) determine the individual system results, 
and we have the restriction that — 1 < aj,bk < 1 . 

Now we seek to find the allowed values of Cjk, given values of aj and 6^. The 
only constraint we have left is that the probabilities of outcomes need to be positive. 
Taking this view will lead us to the allowed values of Cjk- However, there is a simpler 
method. If we use that aj = (Fj^) and bk = {F^), then the outcomes at A and B 
can be thought of as instances of random variables R{aj) and R{bk)- The average 
value (FjFk) = {fhj ®mk)^ P^^ = Cjk must be equal to {R{aj)R{bk)), giving us that 

l^j + bk\ — 1 < Cjk < 1 — \aj — bk\ (3.45) 
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by Lemma 1. So we find that two system states in GNST have the form in Eq. (3.44) 
with —1 < aj,bk,Cjk < 1 satisfying Eq. (3.45). States of multiple systems in other 
theories are a subset of these states. 

With this formalism of generalized probabilistic theories, we can now compare 
quantum mechanics to other theories which allow a greater violation of CHSH. We 
write, for instance, the state of a PR box as 

^ / Af M 

^^=''=' (3.46) 

+ rn-i ® mi + mi ^2 + ^2 ® ^i — ^2 ® ^2 

which satisfies all the constraints on coefficients to be in GNST. We can compare 
this to a quantum state, which we call the Bell state, that maximally violates the 
CHSH inequality, 

^^=''=' ^ (3.47) 

1,. ^ -. -. -. -. -. ^ \ 
+ —={mi ®mi + mi®m2 + m2®mi—m2® m2) I , 

which also satisfies the needed constraints. Now it is absolutely clear that the cor- 
relations in the quantum state are simply weaker than the PR box correlations by 
exactly a factor of 1/a/2. 

There is considerably more to this framework than what is presented here. We 
have illustrated those aspects of generalized probabilistic theories which can be used 
to answer questions about local realism, which is the topic of interest for this work. 
Refs. [7, 8, 10, 60] contain a lot of results particularly regarding quantum information 
processing tasks in generalized probabilistic theories. 

Using this general framework might lead us to a better understanding of quan- 
tum correlations, since now we can compare quantum mechanics to other theories 
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that allow stronger correlations. For instance, an interesting question is to consider 
the theory which allows only locally quantum states and measurements but allows 
arbitrary correlations between systems. How does this theory differ from quantum 
mechanics? This formalism allows one to ask many interesting questions regarding 
the power of quantum correlations. 

3.4 Concluding remarks 

This chapter discusses aspects of local realism and the quantum mechanical violation 
of it. We introduce the random variable and probability distribution representations 
of local realism, and use both to demonstrate different results. Using probability 
distributions we come up with a linear program that solves the problem of whether 
a set of correlations, usually thought of as being generated by measurements on a 
quantum state, can be reproduced in a local and realistic way. If they can be, the 
program generates an appropriate probability distribution. Otherwise, the program 
generates a Bell inequality that is violated; a Bell inequality being a correlation 
inequality that all local and realistic theories must satisfy. Since this program scales 
exponentially in the number of measurements, we then use random variables to 
efficiently generate a set of Bell inequalities that apply to any number of parties and 
any number of dichotomic measurement settings. 

We finish the chapter by introducing a general probabilistic model in which local 
realistic theories and quantum mechanics can be discussed on equal footing. We 
find that this model also permits other non-signaling theories which violate Bell 
inequalities more than quantum mechanics. The question thus arises of why nature 
does not permit such super strong correlations. 

Current research in Bell inequalities suggests that nonlinear Bell inequalities may 
be more powerful than linear ones. Usually, a single nonlinear Bell inequality is either 
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capable of replacing many linear Bell inequalities [100], or it is capable of demon- 
strating a violation of local realism using far fewer measurement settings [91]. It is 
interesting to note that random-variable Bell inequalities are nonlinear. Therefore, 
future research that investigates nonlinear Bell inequalities seems appropriate. On a 
related subject, it would be exciting to know if one could supplement random- variable 
Bell inequalities with either one additional Bell inequality, or possibly a small set of 
Bell inequalities, so that the new set of Bell inequalities is sufficient. Doing this even 
for the case of two qubits would be remarkable. 

There are many open problems in generalized probabilistic theories. Perhaps 
the most interesting is to consider that theory which insists that states be locally 
quantum mechanical, but allows arbitrary correlations between local systems. By 
this we mean that the reduced state is that of a qubit, and furthermore that positive 
probabilities are assigned to all joint-measurement outcomes. States in this theory 
can be represented by density matrices that, instead of having positive expectation 
values on all states, have positive expectation values on all product states. If this 
theory does not violate any Tsirelson bound [ ' ], then quantum mechanics would 
essentially be as nonlocal as possible; if this theory does violate some Tsirelson 
bound, then quantum mechanics is not as nonlocal as it could be. Preliminary 
results suggesting the former statement is true for two-level systems can be found 
in Ref. ['*]. Either way, an answer to this question would provide deep insight into 
quantum nonlocality. 
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Stabilizer states and local realism 



The last two chapters on stabihzer states and local realism have both been building 
up to this chapter, whose namesake is the title of this work. In this chapter we 
will study the correlations present in outcomes of Pauli measurements on stabilizer 
states. Pauli measurements are a good class of measurements to start with because 
their outcomes exhibit correlations that violate local realism [53, 67, 83], yet they are 
well understood within the stabilizer formalism [ ] as per Subsecs. 2.2.3 and 2.3.7. 

The general problem we consider is as follows. Suppose that an ra-qubit stabilizer 
state, which we will assume corresponds to a connected graph, is distributed between 
n parties so that each qubit is spatially separated from the rest. On each qubit we 
assume that either a measurement of /, X, y , or Z is made, with / corresponding to 
no measurement at all. Once this is done, we imagine an independent investigator 
collecting a subset of the results and multiplying them together. This subset of 
results may be correlated in a particular way, and it is these correlations we wish to 
study. 

The way we study the correlations present in Pauli measurements breaks up into 
two main ingredients that should be clearly distinguished. The first is that we begin 
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with a local realistic description of all measurement results. This means that we 
attach an LHV table to each qubit which states what each measurement result, X, 
Y, and Z, will be. It turns out that an essentially complete description of all possible 
LHV tables is possible. However, these LHV tables cannot correctly reproduce all 
measurement correlations because no realistic description of all Pauli measurements 
generally exists. Thus we need to augment our LHV tables somehow in order to 
faithfully reproduce all quantum correlations. 

We do this with our second ingredient, classical communication [74, 80, 91, 92]. 
Specifically, after the measurements are made we allow some of the measurement 
choices to be communicated to other parties. With this new information parties 
are allowed to change their measurement outcomes so as to reproduce the correct 
correlations for the measurements that were made. By examining how many mea- 
surement outcomes need to be sent and to how many parties they need to be sent, we 
learn how difficult it is to reproduce quantum correlations. If LHV tables performed 
well, then we would only need a little communication. On the other hand if LHV 
tables necessarily get many correlations wrong, then it is reasonable to assume that 
a complicated communication strategy will be needed to fix them. 

A model of measurement correlations that uses a local-hidden-variable table and 
supplements it with classical communication is called a classical- communication- 
assisted local-hidden-variable model, or more simply just a communication model. 
Classical-communication-assisted LHV models were pioneered by Toner and Ba- 
con [ ], and Tessier et al. [ ]. 



4.1 LHV tables for stabilizer states 

In this section we learn about constructing LHV tables for Pauli measurements on 
stabilizer states. Of course no table will reproduce all the predictions of quantum 



112 



Chapter 4. Stabilizer states and local realism 

mechanics, as can be seen with the simple example in Ref. [07]. Suppose we have 
an LHV table that provides results of X and Y measurements on all the qubits of a 
GHZ state, 

|GHZ) = i=(|000) + |lll)). (4.1) 

Among the stabilizer elements of this state are —X^Y^Y, —Y^X^Y, —Y^X^Y, 
and X ® X ^ X, which only involve the measurements under consideration. As 
in Subsec. 3.2.3, this gives four conditions on the random variables that generate 
instances of our LHV table, 

R{{X,))R{{Y,))R{{Ys)) = R{{X,Y2Ys)) = R{-1) = -1 
Ri{Y,))Ri{X,))Ri{Y,)) = Ri{Y,X,Y,)) = R{-1) = -1 
R{{Y,))R{{Y2))Ri{X3)) = R{{Y,Y2X,)) = R{-1) = -1 
R{{X^))R{{X2))R{{X;)) = R{{X,X2X;)) = R{1) = 1. 

Multiplying these four equations together yields 1 = — 1, and so a simultaneous 
solution to these equations does not exist. Thus any LHV table must get at least 
one of these wrong. 

Even though we cannot hope to get all measurement results correct, LHV tables 
can capture much about the quantum correlations. This is because the correlations, 
even though stronger than local realism allows, are still rather simple. As discussed 
in Subsec. 2.2.3, a tensor product of Pauli matrices, g, can be classified as one of 
two types for a given stabilizer state. If ±g is in the stabilizer group, meaning 
that ±g = gl"" ■ ■ ■ g'^ where gi, . . . ,gn are canonical stabilizer generators, then the 
result of the measurement is +1 or —1 with unit probability. We will call these 
definite measurements. In all other cases the measurement gives ±1 with equal 
probability; hence we call these random measurements. For example, if the stabilizer 
state corresponds to a connected graph, then no individual Pauli measurement or its 
opposite is in the stabilizer group, and therefore all individual Pauli measurement 
outcomes are random. 
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It is useful here to reiterate the results of Chapter 2 regarding products of Pauli 
group elements. The product of a set of Pauli group elements can always be written 
as in Eq. (2.20), where Eq. (2.21) is the sign that sits in front of the Pauli product. 
For the measurement of a stabilizer element g, the result of the measurement is 
+ 1 with unit probability. For such stabilizer elements, we generally consider the 
definite measurement siga{g)g, because this is a Pauli product with a plus sign in 
front. The outcome of a measurement of sign{g)g is sign((yf) = (—1)'*" with certainty. 
Throughout this chapter we will also be interested in the LHV prediction for a 
definite measurement sign{g)g, called the value oi g, which we write as (— l)''^. We are 
careful to distinguish between LHV results and quantum results, giving the outcomes 
different symbols and different names. If Sg =2 Vg then the LHV table and quantum 
mechanics agree on the result of that measurement, while if Sg 7^2 Vg then they 
disagree. 



Our goal is to find LHV tables that, since they cannot capture all measurement 
correlations, assign values of ±1 to all definite measurements and assign a uniformly 
random variable -R(O) to all random measurements. This will ensure that all random 
measurements are correctly assigned random values. The way such a table could 
make a prediction contrary to quantum mechanics is either to give a value of +1 to 
definite measurements needing a —1 outcome, or to assign a —1 value to those needing 
to be +1. We say that tables with this property are probability preserving in the 
sense that they assign outcomes with the right probability, either certain or random, 
but they may assign the wrong outcome with that probability. In what follows we 
explain a standard method for constructing probability-preserving LHV tables for 
Pauli measurements on stabilizer states, and then we show that any probability- 
preserving LHV table can be constructed in this way. 
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4.1.1 Standard LHV table construction 

This method for constructing LHV tables is originally due to my advisor, Carlton 
M. Caves, and the material presented in this subsection is more or less a review of 
his document on the subject. Once we review how this construction works, we can 
then prove some novel results about the tables. 

As usual we begin with random variables corresponding to individual Pauli mea- 
surement results, and then we derive how to correlate them appropriately. The way 
we write the individual Pauli results is 

RiiX,)) = (-l)-S RiiZ,)) = i-ir^ R{{X,))R{{Y,))R{{Z,)) = {-IfK (4.3) 

The notation above writes the outcomes as —1 to a power, with the power being 
50-50 random or definite, depending on the stabilizer state. Essentially, tables are 
created by assigning values to X and Z, and then, instead of assigning a value to Y, 
by stating the fixed correlation of X, Y, and Z. 

Using these values for the individual measurements, we can calculate the values 
assigned to definite measurements under the table. The first step is to find the 
values assigned to the definite measurements corresponding to our canonical stabilizer 
generators. This can be done by writing the part of the stabilizer generators that is 
a tensor product of Pauli matrices in binary notation, 

N 
k=l 

If we use the individual Pauli results in Eq. (4.3) then the value assigned to this 
tensor product of Pauli matrices, which is labeled as (— 1)*'^, is 

N 
k=l 

At this point, if we wanted to get right the outcome of sign{gj)gj, we would set this 
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equal to Sj modulo 2. However, we want to give LHV tables the freedom to assign 
any value whatsoever to sign{g j)gj. 

The exponent of Eq. (4.5) gives the value of Vj. 

^i =2 Yl i(^k[ri{gj)]k[r2{gj)]k + Xk[ri{gj)]k + 2fch(t/j)]fc) 

\ (4.6) 

K k / 

where x and z are binary vectors whose entries are x^ and z^- Now we simplify this 
by introducing a binary vector, V, whose jth entry is Vj + X]fc'^fc[^i(fi'i)]fc[^2(fi'j)]fe- 
Note that V depends on the generators as well as the values of Ck- What we get is 
that 



V=,^x z )G\ (4.7) 

with G being the generator matrix for the stabilizer state. 

What we have done is to relate V, which involves the values assigned to the 
signless stabilizer generators as well as the correlation between X, Y, and Z, to 
[ X 2 ) by using the generator matrix G. Now, if we could invert this equation 
to solve for i x z ), then by picking the values we want to assign to the signless 
generators as well as the values of Ck, we could actually solve for the required values 
of X and z. These values of x and z, along with the chosen c^, would then determine 
a table that assigns the desired values Vj to the signless generators. Unfortunately, 
solving for f a; ^ j in Eq. (4.7) is problematic because G does not have a unique 
inverse. 

We solve this problem by using the dual generator matrix, H, introduced in 
Subsec. 2.1.4. Specifically, since the dual generators, hj, give rise to random mea- 
surements, we assign them random values (—1)^^. This gives us a binary vector R 
whose jth entry is Vj + Y2k '^k[f^i{hj)]k[f^2{hj)]k, and we have the analogous version of 
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Eq. (4.7) 

R=,[x z)h^. (4.8) 

Combining this equation with Eq. (4.7) gives us that 

V r)=2[x ^ ) ( G^ H^)- (4.9) 

Now we can use the inverse of ( G^ H^ ) in Eq. (2.32) to invert Eq. (4.9), 



V R)i G"" H^ ) =2i V r)a\ a. (4.10) 




X z 



This, of course, works the same for any set of generators since dual generators exist 
for any set of given stabihzer generators. However, we only consider, here, the case 
of canonical stabilizer generators. 

With this result we can construct an LHV table once we are given three quantities. 
First, we need to know what values are to be assigned to tensor products of Pauli 
matrices corresponding to a set of canonical stabilizer generators, Qj^ for the state 
under consideration. For instance, we can choose to get all the generators correct 
by setting (—1)^^ = sign(5fj), but note that this does not guarantee (—1)^" = sign(5f) 
for all stabilizer elements. Then we need to know the random values assigned to 
the dual generators hy Since these generators are not in the stabilizer group for the 
state, they must be assigned random outcomes. Finally, we need to know the value 
of Cj in Eq. (4.3) for all j. Once we know these values, we can construct an LHV 
table that respects them. 

Standard LHV tables for Pauli measurements on stabilizer states are those tables 
constructed by picking values for Vj and Cj and using Eq. (4.10) to determine all 
individual results. We now prove the following result concerning such tables. 

Lemma 3. Standard LHV tables for Pauli measurements on stabilizer states are 
probability preserving. 
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Proof. When we say that standard LHV tables are probabihty preserving, we mean 
that they assign random values to all random measurements and fixed values to all 
definite measurements. 

First we observe that the 2n rows of are linearly independent, where n is 

the number of qubits in the stabilizer state for which we want to construct the table. 
As a consequence, any n-iold tensor product of Pauli matrices can be written as a 
product of generators Qj and hj, possibly with a sign. Since only stabilizer elements 
correspond to definite measurements, if the decomposition of the measurement g 
contains any hj, quantum mechanics will predict its result to be random. Thus 
what we need to prove is that a measurement of a tensor product of Pauli matrices 
g = ±gi^ . . . gn"hi ■ ■ ■ h^n is random if and only if some hj is nonzero. 

We can write down the value of Vg for a measurement of sign{g)g in the same 
way we wrote down the value of Vj in Eq. (4.6), 



^5=2 



^Ck[ri{g)\k[r2{g)\k + ^Xk[ri{g)\k + ^Zk[r2{g)\k- (4.11) 

fc=i fe=i fc=i 

We want to write this in terms of the values of the generators, Vj and r^, because 
those are the quantities we know. So by writing the binary vector r[g) in terms of 
r{gj) and r{hj) we get that 



^9=2 



^Ck[ri{g)]k[r2{,g)]k + ^Xk \^aj[ri{gj)]k + hj[ri{hj)]t 

k=l fc=l \j=l 

n / " 

+ 5Z ^M 5Z "iN(^i)]fc + bj[r2{hj)]k 

k=l \j=l 



n n / n 

"2 



(4.12) 



^Ck[ri{g)]k[r2{g)]k + ^aj i ^Xk[ri{gj)]k + ^Zk[r2{gj)]k 

k=l j=l \fc=l fc=l 



+ X1^J I ^^k[ri{hj)]k + ^Zk[r2{hj)]k 

j=l \k=l k=l 
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Now we can use Eq. (4.6) and its counterpart for rj to write Vg in terms of the known 
quantities Vj, Vj, and Ck, 

n 

Vg =2 ^Ck[ri{g)]k[r2{9)]k + ^aj i Vj + ^Ck[ri{gj)]k[r2{gj)]k 

k=l j=l \ k=l 

n / n 

+ Zl^M "^J + ^Ck\ri{hj)]k[r2{hj)]k j (4.13) 

j=l \ k=l 

Terms that do 



not involve r^s / j=i 



'j'^ji 



Now it is easy to see that the result is random if and only if there exists a hj 7^ 0. D 

To summarize, we can now construct LHV tables, which assign whatever values 
we want to generators, that produce definite results for all definite measurements and 
random results for all random measurements. What we will discover now is that our 
method for constructing LHV tables allows us to construct any probability-preserving 
LHV table in existence. 

4.1.2 All probability-preserving LHV tables 

Suppose we are given a probability-preserving LHV table that assigns values to 
all individual Pauli measurements on some stabilizer state. In this subsection we 
will prove that a standard LHV table exists that reproduces all correlations of the 
given table. Thus standard LHV tables are capable of representing all probability- 
preserving LHV tables. 

Lemma 4. For any probability-preserving LHV table describing Pauli measurements 
on a stabilizer state, there exists a standard table with the same correlations. 

Proof. Since the given table assigns values to all individual measurements, we can 
consider one particular instance of the table and write, for that instance, the values 
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of X, Y, and Z measurements, 

R{{X,)) = (-1)4, R({Z,)) = (-1)4, R{{X,))R{{Y,))R{{Z,)) = (-1)^'.. (4.14) 

where x'j,z',,c'j G {0,1}. Although all three of these variables might be random in 
the LHV table, primed values are fixed because they correspond to one particular 
instance of the given LHV table. The reason I consider one instance is that, in order 
to construct a standard table, all I need to know is the values of the outcomes for 
definite measurements corresponding to stabilizer generators, and these values are 
fixed for all instances of the table. Therefore the value of Vj for this instance must 
be the value of Vj for all instances, and so we get that 

^i =2 Yl (4[^i(^j)]fch(^i)]fc + x'k[ri{gj)]k + 4h(^j)]fe) • (4.15) 

k 

The standard LHV table we want to construct uses this value for Vj, and has c^ = c'^. 
From this it follows that, for all j, 

Yxk[ri{gj)]k + Zk[r2{gj)]k =2 J]4[n(t/i)]fc + 4h(t/i)]fe, (4.16) 

k k 

It only remains to prove that this standard table assigns the same values to all 
definite measurements as the given one, and we have proven what we wanted. 

So consider the value v' of a definite measurement under the given table, once 
again decomposing its binary vector into a linear combination of binary vectors for 
stabilizer generators. 



n 
V 



YciXri{g)]k[r2{g)]k + Y{XfJ(ri{g)]k + 4h(^)]* 



9 

k=l k=l 



5^4[n(^)]feh(^)]fc + ^4 I 5^aj[ri(^j-)]* 

fc=i fc=i \j=i 

(4.17) 

fc=i V j 

n n / n 

fc=l j=l \k=l 
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Now we can use Eq. (4.16) along with Ck = c'^. to write Vg in terms of standard LHV 
table quantities, 



^3 



^Ck[ri{g)]k[r2ig)]k + ^aj i^Xk[riigj)]k + Zk[r2igj)]k 
fc=i j=i \fc=i 

n n / n 

fc=l fc=l \j=l 



(4.18) 



fc=i \j=i 

n n 

=2 ^Ck[ri{g)]k[r2ig)]k + ^ixk[ri{g)]k + Zk[r2{g)]k) =2 Vg- 
fc=i fc=i 

This says that v' =2 Vg, so the standard table assigns the same values to all definite 
measurements as does the given table. Since the standard table automatically gives 
random results to random measurements, this completes our proof. D 

The essence of the above proof is that all definite measurement values are deter- 
mined by the signless generator values, along with the values of c^. Since these are 
variables in creating a standard table, standard tables are capable of representing 
any probability-preserving LHV table. 



4.1.3 One additional measurement 

For universal measurement-based quantum computation using stabilizer states, it is 
crucial that one additional measurement, other than a Pauli measurement, be in- 
cluded. Since our investigation in this chapter is based on probability-preserving 
LHV tables, the existence of probability-preserving tables for non-Pauli measure- 
ments deserves some attention. Measurements involving non-Pauli operators are no 
longer either definite or uniformly random, so we first make explicit what is meant 
by probability preserving in this case. Probability-preserving tables for Pauli mea- 
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surements are defined so that (M)lhv = ±(^)qm, and so we take this as the general 
definition of a probabihty-preserving LHV table for non-Pauli measurements as well. 
Unfortunately, we prove in this subsection that such tables do not exist, in gen- 
eral, when a non-Pauli measurement is allowed. Therefore, extending our results 
to include a non-Pauli measurement seems likely to require more complicated prob- 
abilistic communication strategies. Whether probability-preserving tables exist for 
Pauli measurements on more general states [31] is still open. 

We now work through an example where the construction of probability-preserving 
LHV tables fails when any non-Pauli measurement is allowed. In this vein, consider 
the stabilizer state 

|V') = ^(|00) + |11)) (4.19) 

and an arbitrary single-qubit measurement T{6, (p) = u{6, (p) ■ a, for some unit vector 
u{6, 0). The angles 6 and used to describe the vector are not necessarily the usual 
spherical polar angles. To define 6 and (p, we need to first define new axes x', y', and 
z' as follows. First write u in terms of the unit vectors ±x, ±y, and ±z, such that 
the coefficients of u are not negative. For example, ii u = —x, then u = +1(— x). 
Once this is done, we choose z' as the axis with the smallest positive coefficient and 
x' as the axis with the largest positive coefficient. Notice that a primed axis, a' , 
need not be ±a and that the resulting axes need not form a right-handed coordinate 
system. In the example that u = +1(— a;), we must choose x' = —x, but we can 
choose z' to be either ±y or ±z. If we choose z' = y and y' = —z, then the resulting 
coordinate system is not a right-handed one. 

The angles 6 and are illustrated in Fig. 4.1 and are simply spherical polar 
angles with respect to the primed axes. Specifically, 6 is the polar angle from the 
z' axis and is the azimuthal angle about the z' axis, measured in a right-handed 
way for a right-handed coordinate system and in a left-handed way for a left-handed 
coordinate system. The primed axes are defined so that u lies in the tetrahedron in 
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Figure 4.1: Definitions of 6 and (p once tlie axes x', y', and z' liave been identified. 
The primed axes shown result in a right-handed coordinate system, but this need 
not be the case in general. Since u is furthest from z' and closest to x', it must lie 
in the tetrahedron outlined in the figure. 



the figure. The fact that the vector u lies in the tetrahedron constrains the angles 
so that 



0<(j)<- ande^<e <-, 

-Y - ^ 0- -2' 



(4.20) 



where 6*,^ satisfies cos 9^ = sin 9^ sin 0. The upper face of the tetrahedron lies along 
the plane z' = y' ^ cos 9 = sin 6' sin 0, hence 9^ is the smallest allowed value for 9. 

Now consider measurements on the state {ip) of X' (8> X', Y' ® X', Y' (g) T{9, (p), 
and X' (g) T{9,(f)), where X' and Y' correspond to measurements along x' and y', 
respectively. We calculate that 



(X' ® X') = ±1 
(r'®T(^,0)) = ±sin^sin( 



(Y' 0X') = O 
(X'(g)T(^,0)) = ±sin^cos. 



(4.21) 



We want to show that it is impossible to preserve these quantum probabilities with 
an LHV table unless T{9, (p) is proportional to a Pauli matrix. 

Following Subsec. 3.2.3 we want, for a probability-preserving LHV table, random 
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variables so that 

R{{X[))R{{X',)) = ±1 

R{{Y;))Ri{X',)) = R{0) 

(4.22) 

R{{Y;))R{{T{e))) = ±i?(sin^sin0) 
R{{X[))R{{T{e))) = ±R{smecos(f)). 

Multiplying these four equations together, and absorbing all signs into -R(O), gives 
us the following consistency condition. 

i?(sin 9 cos 0)i?(sin 6 sin 0) = R{0) (4.23) 

Using Theorem 1, what we want to know is the values of 6 and satisfying the 
inequality, 

I sin6'cos{/) + sm6sin(j)\ — 1 < < 1 — | sin0cos(/) — sm6sm(j)\. (4.24) 

For the allowed values of 6 and (p given in Eq. (4.20), our claim is that | sin6'cos{/) + 
sin 6' sin 01 — 1 < is true if and only if 6' = = which is the case that T{6, 0) is 
proportional to a Pauli matrix. 

First, we rewrite the left hand side of Eq. (4.24) as |cos0 + sin0| < 1/sin^, 
where we use that sin^ > supposing that 6^0. Now from Eq. (4.20), we get that 
1/sin 9 < 1/sin 9^ = y 1 + sin^ 0. So let us compare | cos + sin 0| and y 1 + sin^ 0. 
Squaring them both gives | cos0 + sin0p = 1 + (2cos0) sin0 and (y 1 + sin^0)^ = 
1 + (sin0) sin0. Supposing that 7^ and given that < < 7r/4, it is clear that 



2 cos > sin =^ I cos + sin 0| > -w 1 + sin > 1/sin 9. (4-25) 

Hence we have a contradiction with Eq. (4.24), and we conclude that | sin ^ cos + 
sin 9 sin 0| — 1 < holds true only for = 0^^ = 0. Therefore we have shown that 
the probabilities associated with the correlations in Eq. (4.21) can be preserved only 
if T{9,(f)) is proportional to a Pauli matrix. 



124 



Chapter 4. Stabilizer states and local realism 

4.1 A Comparing LHV tables and quantum mechanics 



It is an interesting question to ask how well our standard probability-preserving tables 
for Pauli measurements perform, a task that can be accomplished, for example, by 
calculating the probability of failure. However, in order to do this, we need to know 
the exact prediction of all definite measurements under the table. We can then 
compare these predictions to the quantum mechanical predictions and see when they 
agree or disagree. 

We learned in the beginning of this section that quantum mechanics predicts the 
result of a definite measurement, Mg = sign{g)g, to be sign(f7) = (— 1)"*®. So, we 
actually already know the quantum mechanical prediction of definite measurements 
from Subsec. 2.1.3. We are interested in the sign of products of canonical generators 
9i^ ■ ■ ■ Qn" given in Eq. (2.36) of Result 5. What we will now find is the corresponding 
formula for standard LHV tables so that we can compare measurement outcomes. 

To determine the result, (—1)^®, assigned to a definite measurement under a 
standard LHV table, insert the generator values in Eq. (4.6) into the value of Vg in 
Eq. (4.18), thereby eliminating Xk and z^. Doing this gives us that 



n n 

i=\ j=i 



(4.26) 



T 

av 



+ ^ci[ri{g)]i[r2{g)]i + ^ ajCi[r^{gj)]i[r2{gj)]i. 
1=1 j,i=i 



Now we want to write the binary vector for g a.s a. linear combination of the binary 
vectors for gj. When we do this, the quantity Cjk = (X]r=i Ci[^i(fl'j)]«['"2(fl'A:)]«) appears 
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several times. 



Vg =2 av 



n In 

.T 



1=1 \j=l J \k=l 

n 

+ ^a,Q[ri ((?,•)]« hfe)]/ (4.27) 

n n 

= av + 2_. '^j'^kCjk + /_.'^fijj =2 Ci'tJ + / ^O'jO'kCjk- 
j,k=l j=l jj^k 

What we have here is that the value under the standard LHV table of a definite 
measurement of g is the sum of two quantities. The first, rather expectedly, is 
the product of the values for definite measurements corresponding to the generators 
needed to obtain g, (—1)"'' . The other quantity involves looking at the generators 
involved in obtaining g and carefully taking into account the value of Ck for those 
generators. 

To compare this with the quantum mechanical result, we should write the LHV 
measurement outcome in terms of r{g) by replacing the binary vector a with r{g) 
via Eq. (2.34). The result is most simply written in terms of C, which is C with its 
diagonal entries set to zero. In that case, 



Vg =2 av 



+ aCa^ = r{g){vHAf + r{g) {AH^) C (HA) r^{g). (4.28) 



To carry out the calculation of this quantity, let us write C = G1CG2, where G 
[ Gi G2 ) and C is a diagonal matrix whose entries are the values c^, i.e Cjk 
CkSjk- Using the canonical generator matrix, we have 



C = GiCGl 



I a\ C('-) 

y I (:7("-^) 
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where r is the left rank of the generator matrix G, and we spht up the matrix C into 
its first r and last n — r rows and columns. Now we get C by making the diagonal 
entries zero, 

C=\ , (4.30) 

\ ) 

where B is obtained from B by setting its diagonal entries to zero. Using this matrix, 
we finally get the outcome of an LHV table. 

Result 6. Let Mg = sign{g)g be a tensor product of Pauli matrices corresponding 
to a stabilizer element g = g^" ■ ■ ■ g^"- A standard LHV table, defined by a binary 
vector V and a matrix C with Cjk = CkSjk, assigns to this measurement an outcome 

{-iy\ with 

Vg =2 av'^ + a^'^) ( C^^)B C^"-^ A + AC^'^'"^ ) a^ 

^ ' (4.31) 

= r{g)[vHKf ^r^{g)yC^-)B C^A + AC(^-') j r^(^). 

The superscript r or n — r indicates whether to take the first r or last n — r entries, 
and we have Bjk = Bjk{l — Sjk)- 

Now, in order for an LHV table to predict all the same results as quantum 
mechanics, one needs to have that Vg =2 Sg for all g. For every g for which this is 
not true, that is a measurement the LHV table gets wrong. 

It would be nice to know the best LHV tables for Pauli measurements on stabilizer 
states. However, the answer to this question is not known for the general case. 
There are some preliminary results in this regard. For instance, permutations of C 
corresponding to a symmetry of the graph yields a table that gets the same number 
of outcomes correct. It is also known that given a table, another table that assigns 
opposite results to all measurements does not exist. Results like these lead towards 
understanding how many correct results a table can predict, but since we do not 
have a general answer we omit the proof of these claims. 
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4.1.5 Equivalent standard LHV tables 

Every probability-preserving stabilizer LHV table has a standard counterpart, how- 
ever it can happen that two apparently different standard LHV tables actually yield 
identical correlations. We now focus on when this happens using the tools we just 
developed to calculate LHV measurement outcomes. 

Clearly for each of the 2"' ways to assign canonical generator values we get a 
unique table. For a given choice of the Vj, there are also 2" choices for the Cj, 
however it is possible that two different choices of Cj yield equivalent tables. If we 
just consider the special case that the stabilizer state is a graph state, it turns out 
that flipping all the Cj yields an equivalent table, and we can prove that this leads 
to exactly a factor of 2 redundancy. 

Lemma 5. For an n-quhit graph state corresponding to a connected graph, there are 
2n2'n.-i _ 22"-i ufiiqiiQ probability-preserving LHV tables. 

Proof. First of all, specializing the value of Vg in Eq. (4.31) to graphs states gives us 

Vg =2 av'^ + aCBa^ = av^ + aBCa^, (4.32) 

since in this case r = n and GiCG^ = ICB^ = CB already has all zeroes on the 
diagonal. 

So now we begin with two tables, one table that assigns values to the stabilizer 
generators according to v and values to the product of X, F, and Z measurements 
according to c and the other table that assigns values according to w and d. The 
idea is to show that the two tables assign the same values to all stabilizer elements if 
and only iiw = v and either d = c oi d =2 c+1, where c + 1 adds 1 to each element 
of c. 
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First, suppose w = v and d = c+ 1. Then using Eq. (4.32) with Djk = dk6jk, we 
see that 

aw^+aBDoF = av^ +aB{C+I)a^ = av^+aBCa'^+aBa'^ =^ av^+aBCa^, (4.33) 

where Bj^ = B^j was used in determining that aBa^ is even, and therefore equal to 
modulo 2. So these two tables assign the same outcomes to all measurements. 

Now suppose that the two tables assign identical values to all stabilizer elements. 
If we write w =2 v + u and D =2 C + P, then this condition becomes that, for all a, 

au^ + aBPj =2 0. (4.34) 

This is the condition that the table specified by u and P assign +1 to all stabilizer 
elements. If we let a = Cj, i.e., the measurement is of stabilizer generator Qj, then 
this becomes 

n 

=2 Uj + ^ SjkBkiPaSip5pj = Uj + BjjPjj = Uj. (4.35) 

k,l,p=l 

So we get that u = ^ w = v, which is the first part of what we needed to show. 
What is left is to find all possible P such that aBPa^ =2 for all a. The claim is 
that either P = or P = I. To prove this, suppose that there were two nodes such 
that Pjj = and Pkk = 1- Now we can assume that Bjk = 1 because the graph is 
connected, and so there must be two connected nodes that have different diagonal 
values of P. Otherwise all nodes would have either Pjj = or Pjj = 1, contradicting 
our assumption. So we find that for a = Cj + Ck, 

n 

=2 aBPa^ = ^ {5ji + 5ki)BipPpp5pq{5gj + 5qk) 

l,p,q=l (4.36) 

~ BjjPjj + BjkPkk + BkjPjj + BkkPkk = 1- 

Hence the contradiction, and we conclude that P = OotP = I^D = C or 
D =2 C + I, as desired. D 
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We also have the following corollary. 

Corollary 3. For an n-quhit stabilizer state corresponding to a connected graph, 
there are 2^""-*^ unique probability-preserving LHV tables. 



Proof. First assume that two tables assign the same values to all definite measure- 
ments. From Eq. (4.31), and following previous notation, we get that 

(4.37) 
^ ua^ + a^'-'^P^'-^Bia^'-'^Y + a^'^P^"-^ A + AP(^-'-))(a(^-'-))^ =, 0. 

As before, setting a = Cj and a = Cj + e^ for some neighbors j,k < r we get u = 
and P*-''^ = or P^"^^ = L The resulting constraint is now 

r n 

=, a^^\P<^^'>A + AP(^-))(a(^-'-))^ = S S ci,Ajka,{P,, + P^k) (4.38) 

j=l k=r+l 

for all a. But by connectivity, there is some Ajk = 1 so for that j and k we get 
Pjj = Pkk- Now since all the Pjj are equal, we conclude that all the Pkk are also 
equal, which gives us that P = or P = /. 

Similar to before we can check that such choices for P generate equivalent tables, 
and so we get the same redundancy factor of 2 that we had for graph states. D 



We also note that the local-Clifford equivalence of stabilizer states to graph states 
also gives us Corollary 3. 

This concludes a very thorough examination of LHV tables for Pauli measure- 
ments on stabilizer states. We now move on to learning how classical communication 
can correct these tables when they make incorrect predictions. 
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4.2 Simple classical- communication models 

In a communication model we assume that Pauli measurements are made on a sta- 
bilizer state and that the results are obtained from an LHV table. We furthermore 
choose this table to be the probability-preserving tables of the previous section so 
that the decision to change LHV results need not be probabilistic. Measurement 
choices are then communicated between nodes in the graph. The party at each node 
can then, based on this new information, choose to either output their LHV result 
or output its negative. The decision of which choice to make depends on the details 
of the particular model. 

We begin by discussing the simplest kinds of communication models. These 
models share two properties that make them simple. The first is that communication 
take place only between neighbors in the underlying stabilizer graph. This restriction 
to communication only with neighbors in the graph makes intuitive sense if we think 
of a graph as a recipe for constructing the corresponding stabilizer state. In that 
case, nodes that are connected have interacted in the past and therefore occupy a 
privileged position with regard to exchange of information. We call a model that 
restricts communication to neighbors a nearest-neighbor communication model. 

Another simplification we make is to consider site-invariant models. Site invari- 
ance is the property that nodes in symmetric situations perform the same action. 
This property is natural because identical particles that have identical histories 
should behave identically, site-invariant nearest-neighbor communication models, 
i.e., models that obey both of these simplifications, are explored in Refs. [11, 90]. 

Finally, in what follows, instead of considering stabilizer states in general we only 
consider Pauli measurements on graph states whose stabilizer generators all have a 
sign of -|-1. Fortunately, we do not lose any generality by doing this because any 
communication strategy for stabilizer states can be translated into a communication 
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model for a graph state. This follows from the local-Clifford equivalence of stabilizer 
states to graph states and the local-Clifford manipulation of LHV tables described 
in Refs. [89, 90]. In addition, the underlying graph of a stabilizer state describes the 
same interaction pattern as for the local unitarily equivalent graph state, so commu- 
nication in the graph model will still be confined to be along edges between neighbors 
in the graph. However, we do point out that in general this adds the complication 
of having to communicate node type as well, so that each party knows which local 
Clifford operations were applied to the parties with which they are communicating. 

It turns out that these two properties that make models so simple also assure 
their failure. In this section and the next, we consider the effects of each of these 
properties in some detail and show that if a model has either property it necessarily 
fails to capture some quantum correlations. However, it is still helpful to consider 
these models because they give us minimal conditions for other models to work and 
because they can still capture many quantum correlations as explained below. 



4.2.1 Definite measurements and definite submeasurements 

Before getting to a simple communication model, we first review quantum mechanical 
results concerning Pauli measurements on graph states. When these measurements 
are made, the product of all measurement outcomes is the same as the value of mea- 
surement of the corresponding Pauli product. This tensor product of Pauli matrices, 
M, is called a global measurement because it results from considering all the mea- 
surement outcomes that were made. A suhmeasurement of a global Pauli product M 
is a Pauli product M such that the non-identity elements of M all appear in M, i.e., 
Mj = Mj or Mj = I for all j. LHV models implicitly predict the result of measuring 
such a subset of the Pauli operators of a global measurement, the measurement of 
an identity operator being simply the omission of the corresponding LHV entry. A 



132 



Chapter 4. Stabilizer states and local realism 

proper communication-assisted LHV model for graph states should not only repro- 
duce the predictions of quantum mechanics for global measurements but also for all 
possible submeasurements. LHV tables by themselves are capable of reproducing 
neither type of measurement. 

If we consider a global measurement M, which may be random or definite, of a 
tensor product of Pauli matrices then a submeasurement of M can be specified by 
a binary vector e such that Cj = 1 if Mj is included in the submeasurement and 
Cj = if it is not. Given the vector e that specifies the submeasurement, we can use 
[ri{M)]jej and [r2{M)]jej to write the binary vector for the submeasurement. In this 
way we force all measurements with Cj = to be J, and leave the rest unchanged. 

Our models will automatically predict correctly the quantum result of random 
global measurements and random submeasurements because the underlying LHV ta- 
ble assigns random results to those measurements. Thus in all of what follows, we 
concern ourselves only with definite measurements and definite submeasurements. 
Recall that any definite measurement or submeasurement Mg corresponds to a sta- 
bilizer element g, with Mg = sign{g)g. 

The condition to be a definite measurement is simple to calculate. If Mg = 
sign{g)g, then it must commute with all the generators since Mg = sign{g)g1^ ■ ■ ■ g'^ =^ 
[Mg,gj] = 0. Whether Mg commutes with the generators gj can be determined by 
the symplectic product in Eq. (2.9), so this tells us that 



r{Mg)AG^ =. ^ ( n{Mg) r^iMg) ) I ° 








(4.39) 



^ r^iMg) =, r,{Mg)T 



where F is the adjacency matrix that defines the graph state. We can use this result 
to determine whether a submeasurement of M, corresponding to binary vector e, is 
random or definite. We can then recover the above formula in the case that e = 1. 
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So, if we take the binary vector for a submeasurement and plug it in to the above 
condition, we get the condition that the submeasurement be definite. This is an 
important result. 

Result 7. The condition for a submeasurement of M corresponding to e to be 
definite is that 

n 

[r2{M)],e, =, J2 ^jk[ri{M)]kek, (4.40) 

fe=i 

for all J = 1, . . . ,n. 

Now we turn to the quantum mechanical prediction of a definite measurement or 
submeasurement. The quantum mechanical prediction for the outcome of a definite 
submeasurement of a global measurement M which is specified by the binary vector 
e can be determined by plugging in the binary vector of the submeasurement into 
Eq. (2.36) of Result 5. We find this result to be i to the power, 

n n 

J2 [ri(M)],e,r,fc[ri(M)],e, - ^[ri(M)],e,h(M)],e, (4.41) 

j,k=l j=l 

Now if we use Eq. (4.40), denoting quantities that should be taken modulo 2 with 
an overbar, we get the quantum mechanical result of a definite submeasurement. 

Result 8. The quantum prediction for a definite submeasurement of M correspond- 
ing to e is i to the power 



^[n(M)],e, 5^r,fc[n(M)]fcefc - ^r,4ri(M)],e, . (4.42) 

j=l \k=l fc=l / 

4.2.2 A nearest-neighbor model 

In this subsection we present a site-invariant nearest-neighbor communication model 
which predicts correctly the global outcome of the measurement of any Pauli product. 
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This communication model will supplement one particular choice of LHV table, the 
one with v = and c = 0. The outcome of any definite measurement under this 
table can be calculated using Eq. (4.31), which looks rather complicated until it 
is realized that v is the zero vector and C is the zero matrix. Thus we have that 
Vg = for all definite measurements Mg] this table assigns the value +1 to all 
definite measurements. Obviously, then, this table correctly predicts all definite 
measurements that are +1 according to quantum mechanics and incorrectly predicts 
those assigned a —1 outcome. 

This table is one we will use throughout this chapter, and so it deserves a moder- 
ately involved discussion. Since we only consider graph states, so that G = i I T ) 
with r being the adjacency matrix of the graph, we find that the dual generator ma- 
trix H is i I ]. For these generator matrices, and setting f = c = 0, Eq. (4.10) 
reduces to 



r 


^ 


=^ z = r and x =2 rT 



Or) \ ^ z = r and x =2 rT ■ (4.43) 



X z 



Hence we find that the values of z are and 1 with equal probability, and for this 
reason we drop r in favor of z and write x =2 zT. 

In summary, our LHV table uses n binary random variables, zi, . . . ,Zn, each tak- 
ing on values and 1 with equal probability. The random variable (—1)^^ is thought 
of as value for Zj under the LHV table, i.e., the value of the z spin components of 
the jth qubit. For the corresponding values of the x and y spin components, we have 
that 

(-1)^^ = n (-l)'^ 

'^^'^^ (4.44) 

(-1)% = (-1)-. Yl (-1)"^ 

k&Af{j) 
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The values Xj assure that +1 values are associated with the generators gj, i.e., 

(_1)-. Yl {-ly^ = +1 , (4.45) 

k€AfiJ) 

in analogy to Eq. (2.44). The values 

(-1)^^(-1)W(-1)^^ = 1 (4.46) 

give instructions for how to multiply LHV entries, as does XYZ = i for Pauli 
matrices. 

We assume now that each party is given a measurement Mj to perform, chosen 
from /, meaning no measurement, and X, Y, and Z. After the measurement, there is 
a round of communication between neighboring sites, and then each party outputs a 
value +1 or —1 as the result of the measurement. When no measurement is performed 
at a site, the output can be regarded as +1. 

During the round of communication, site j sends a bit bj to each site k G A/'(j), 
where bj = if Mj = I,Z and bj = 1 if Mj = X,Y. The value rrij, for the 
measurement of Mj, output at site j is determined by the LHV table result for the 
measurement at that site and by the quantity 

n 

tj= Y. h = Y,^ikbk, (4.47) 

fceATQ-) k=l 

which is computed from the bits sent to site j from neighboring sites in a graph with 
adjacency matrix F and which is equal to the number of neighboring sites that make 
an X or y measurement. The output rrij is determined by rules that decide whether 
to flip the sign of the LHV entry associated with the measurement at site j: 



1. If Mj = /, ruj = 1. 



2. li Mj = Z, rrij = (-1)^^. 
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(-1)-:' if t, =,0,1, 



3. If Mj = X, rrij 



4. If Mj = Y, mj 



-Ip if tj =4 2, 3. 



-i)y^ if tj =,0,1, 
-l)y^ if tj=,2,3. 



If we take the product of the outputs from all the sites, the result, written as (—1)"'^, 
corresponds to a prediction for the measurement of the operator M = (S)?=i Mj. The 
number of bits communicated in this model is twice the number of edges in the graph. 

4.2.3 Proof that the model works for global measurements 

The communication model introduced in the last subsection reproduces the quantum 
predictions for any global Pauli measurement on graph states. The proof of this fact 
proceeds as follows. First we will determine the result (—1)^*-' of the measurement 
M according to the model. Once this is done we will determine the result of the mea- 
surement according to quantum mechanics, denoted as (—1)^^^. Finally, we compare 
the two quantities and see that sm =2 vm for definite measurements. We will also 
see that quantum mechanics and the model agree for random measurements. 

Let us first determine the result of a measurement according to our model. If we 
write the result of the measurement according to the LHV table as (—1)'''^^, then the 
model says we fiip this value under certain circumstances. According to the rules of 
the model only measurement results of X and Y can change the value of (—1)^'^, 
and they do so only when tj =4 2, 3. This can be summarized as 

where, once again, we have the notation that quantities under an overbar are to 
be taken modulo 2. Notice that tj — tj is even, so we are assured that the result 



137 



Chapter 4. Stabilizer states and local realism 

will either be +1 or —1. For random measurements, (— 1)^^-^ is random according to 
Lemma 3 and therefore so is (—1)^". For definite measurements, (—1)^" = +1 and 
so our model gives the result of a definite measurement to be 

(-1)""= Y[ i'^-^K (4.49) 

{i|bi=i} 

These values for definite measurements are the ones we wish to compare to quantum 

mechanics. 

We want to show that the quantum mechanical result of a global measurement 
M is equal to the result predicted by our model. The key to relating the quantum 
mechanical result to our model is by recognizing the relationship between r{M^ and 
values of 6j and ty First, recall that [ri(M)]j = 1 when Mj = X or F, which tells 
us that [ri(M)]j = 6j, the quantity communicated in the model. Second, if M is a 
definite measurement then Eq. (4.39) must be satisfied, which is to say 

r2(M) =2 ri(M)r 

(4.50) 



^ [r2(M)], =, 5^r,fc6, = t, 



fc=i 



So now we can write sm from Eq. (4.42) in terms of quantities used in the com- 
munication model. Doing so gives the value of Ism to be 



J][ri(M)], 5^r,.,[ri(M)],-5^r,,[ri(M)], =,^6,(t,-t,). (4.51) 

i=i \fc=l k=\ / j=l 

So we find that (—1)'^*^, the quantum mechanical outcome for a definite measurement 
M, can be written as 

{i|bj=l} 

which is exactly the result of definite measurements predicted by our model according 
to Eq. (4.49). 
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Even though we have now completed our proof, we make a few comments here. 
According to Eq. (4.50), definite measurements have the property that tj is even 
when Mj = X, and tj is odd when Mj = Y. This means that we can omit the 
case when tj is odd and Mj = X or tj is even and Mj = Y, because in these cases 
the measurement result is random. The only necessary part of our model is that X 
measurements change their LHV result when tj =4 2 and Y measurements change 
their LHV result when tj =4 3. This means that there are variants of our model 
that also predict correctly the result of all global measurements. One special such 
variant [ ] is that X change its LHV result when tj =4 1, 2 and Y change its result 
when tj =4 0, 3. Of course, this change means nothing for global measurements, 
but when considering submeasurements of a global measurement such variations can 
make a difference. 



4.2.4 Site-invariant communication models 

Both the numbering and the arrangement of nodes in a graph are arbitrary, so it 
seems reasonable to suppose that a communication model should be insensitive to 
these things. We refer to this property as site invariance and define it formally as 
follows. Given a graph, each of whose nodes has been assigned a measurement, 
a permutation that leaves the graph invariant is one that interchanges nodes and 
their measurements, letting edges move with the nodes, such that the new graph 
is identical to the old graph in the sense that they could be placed on top of each 
other with all nodes, measurements, and edges overlapping. A site-invariant model 
is one for which nodes in identical situations, as defined by permutations that leave 
the graph invariant, make the same sign flipping decision. Surprisingly, we find this 
trait to be at odds with the modeling of submeasurements. 

We demonstrate the limitations imposed by site invariance using the example of 
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12 3 



6 5 4 

Figure 4.2: Example demonstrating that any communication model based on the 
LHV table from Eq. (4.44) which is site invariant fails to reproduce some submea- 
surements. The global measurement M = YiY2Y^Y4Y^Yq has a random outcome, but 
contains a submeasurement M = Y1Y2YSI4Y5IQ such that —M is an element of the 
stabilizer group. This means that an overall sign flip is required to correct the +1 
prediction of the LHV table for a measurement of M. The two qubits measuring Y at 
nodes 1 and 3 are in symmetric situations, as are the qubits at nodes 2 and 5. Thus, 
under a site-invariant model, 1 and 3 must make the same sign flipping decision, as 
must 2 and 5. For each pair, the sign flipping decisions cancel one another, producing 
no overall sign flip and thus giving an incorrect result of +1 for the measurement of 
M. 



a 2 X 3 cluster state, which is depicted in Fig. 4.2. The two relevant measurements 
for this example are M = Y1Y2Y3Y4Y5YQ, which has a random outcome, and M = 
YiY2Y^I4Y^Iq, which has the certain outcome —1. When either of these is considered 
as a global measurement, our model yields the correct prediction, as we have already 
shown in general, but when the second is considered as a submeasurement of the first, 
the model fails. In this second case, rules 1-4 say that the two qubits measuring Y at 
nodes 2 and 5 should introduce a sign flip, but the two qubits measuring Y at nodes 
1 and 3 should not. The result is no overall sign flip and an outcome +1, showing 
that the model gets the submeasurement outcome wrong. In contrast, when M is 
considered as a global measurement, rules 1-4 dictate a sign flip for qubit 2, but no 
other qubit, thus giving the correct, certain outcome —1. The same measurement M 
can lead to different sign-flipping decisions in the two situations because the nearest- 
neighbor environments of the qubits differ depending on whether a submeasurement 
or a global measurement is under consideration. As is shown in Figure 4.2, the 
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counterexample is not limited to the communication model used in this paper. In 
fact, any site- invariant model based on our LHV table yields an incorrect result for 
the submeasurement M. 

This example can easily be generalized by adding p rows and q columns to opposite 
sides of the 2x3 cluster state. Doing this results in a class of {2 + 2p) x (3 + 2g) cluster 
states for which LHV models based on the LHV table of Eq. (4.44) and assisted by 
a site-invariant communication model fail for some submeasurement s. 

It may be possible to find a site-invariant model for a different LHV table that 
could reproduce all submeasurement correlations. In this case the symmetry would 
not be broken by the communication, but rather the table. The point is, nodes in 
symmetric situation still must behave differently. It is interesting to point out that in 
this example, the amount of communication and the set of qubits that communicate 
are both irrelevant. Regardless of how much communication we allow, a site-invariant 
model cannot reproduce all submeasurements. 



4.2.5 Graphs for which the model works 

We have shown that our communication model agrees with quantum mechanics for 
correlations involving global measurements. We have also demonstrated that because 
it is site invariant, it will generally fail on submeasurements. In this section we 
demonstrate that there are some classes of graphs for which the model correctly 
predicts not only all global measurement correlations, but all correlations present in 
submeasurements of those global measurements. 

The submeasurement results predicted by our model come from a generalization 
of Eq. (4.49). That equation gave us the result of a global measurement to be i with 
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an exponent 



j=i \k=i fc=i / 

Generalizing this to a sub measurement, described by a binary vector e, is not as 
simple as replacing all b values with b times e; if it were so simple our model would 
predict all submeasurements correctly. In the above equation, bj indicates which 
measurements were present in the measured Pauli product, and b^ indicates which 
nodes sent their measurement information. Therefore, we need to replace bj by bjCj, 
since we are only multiplying together those results with Cj = 1, but we leave bk 
alone because everyone still sends a bit 6^ if they made an X or F measurement 
regardless of whether their measurement is chosen as part of the submeasurement. 
Thus we get that the prediction of a definite submeasurement by our communication 
model is i to the power 



E ^^^^ E ^i^^^ - E r^-^^^ • (4-54) 

i=i \k=i k=l / 

Since we know our model correctly predicts random submeasurements, our task 
is to find graphs F such that Eq. (4.42) and Eq. (4.54) agree for all binary vectors 
e and definite submeasurements satisfying Eq. (4.40). We find two such classes of 
graphs below. 

Lemma 6. Our nearest-neighbor communication model correctly predicts all submea- 
surements for complete bipartite graphs. 

Proof. We know our model predicts correctly random submeasurements, so we only 
need consider definite submeasurements. 

The nodes in a complete bipartite graph can be partitioned into two sets, A and 
B, such that each node in A is connected to every node in B, yet no node in A is 



142 



Chapter 4. Stabilizer states and local realism 

connected to another node in A nor is any node in B connected to another node in 
B. The adjacency matrix for such a graph looks hke 

r.(°;j. (4.55) 

where 1 denotes a matrix with all entries equal to 1. The nodes in the sets A and B 
are clumped together to obtain this nice form. 

The quantum mechanical result of any submeasurement of a global measurement 
M is, from Eq. (4.42), 

y^ bjCj I ^ bkCk - ^ bkCk I + ^ bjCj I ^ bkCk - ^ bkCk | • (4.56) 

jeA \keB keB J jeB VfceA fceA / 

The result according to our model of any submeasurement is, from Eq. (4.54), 

Y^b,eJj2^k~Y.^k]+J2b,eJj2bk-Y,bk]. (4.57) 

Notice that for both the quantum mechanical result and the communication model 
result, the quantities in parentheses are even. Finally, the constraint to be a definite 
submeasurement, Eq. (4.40), becomes for this graph. 



[r2(M)]j6AejeA =2 ^ ^feCfc 

keB 

[r2{M)]j^BejeB =2 ^bkCk- 



'■^"^ (4.58) 



fceA 

Now we are ready to prove that our model correctly predicts all definite submea- 
surements. First, since we already know the model predicts the correct result for the 
global measurement, we can assume that some Cj = 0. Without loss of generality, let 
us suppose this occurs for some j G A. Then Eq. (4.58) says that J2keB bk^k is even. 
Since the product of two even numbers is zero modulo 4, this eliminates the second 
term from Eqs. (4.56) and (4.57). Now there are two cases to consider, depending on 
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whether Cj = for some j E B or Cj = 1 for all j G B. In the first case, Eq. (4.58) 
implies that J2keA ^k^k is even and both quantum mechanics and our model predict 
the outcome to be +1. In the second case, setting e^ = 1 for all A; G -B in Eq. (4.56) 
gives us Eq. (4.57), and so our model agrees with quantum mechanics in this case as 
well. 

D 

Complete bipartite graphs are interesting because they include the star graphs of 
GHZ states [..,]. There is also another class of graph for which our model works. 

Lemma 7. Our nearest-neighbor communication model correctly predicts all submea- 
surements for the symmetric difference of two complete graphs. 

Proof. Once again, we know our model correctly predicts random submeasurements 
so we only consider the definite submeasurements. 

The symmetric difference of two complete graphs [ ] is constructed by beginning 
with a complete graph, meaning a graph in which every node is connected to every 
other node, and removing from it a smaller complete graph on a subset B of nodes. 
The adjacency matrix for this graph looks like 



/ 



V 1 



1 

1 

1 



1 \ 



1 





(4.59) 



/ 



The nodes in the graph can be partitioned into sets A and B, where A includes all 
of the remaining nodes not in B. 
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The quantum mechanical result of any submeasurement of a global measurement 
M is, from Eq. (4.42), 



y^ bjCj I ^ bkCk - 1 - X^ hek ^ 1 ) + X^ ^jCj I ^ 6^6^ - ^ fefcCfc | • (4.60) 

jeA \k=l k=l J j(^B VfeeA fceA / 

The result according to our model of any submeasurement is, from Eq. (4.54), 



J^b.e, [J2^k-l-J2^k-l] +J2^je, 5^6. - 5^6. . (4.61) 

j&A \k=l k=l ) j&B \keA k&A / 

The constraint to be a definite submeasurement, Eq. (4.40), becomes for this graph, 

n 

[r2(M)]j6AejeA =2 ^hek 

^=1 (4.62) 

[r2{M)]j(.BejeB =2 ^hek- 

Now, as before, there must be some Cj = and there are two possibilities, j & A 
and j G B. If Cj = for some j G B, then Eq. (4.62) says that Ylik^A^k^k is even. 
This makes the first term in Eqs. (4.60) and (4.61) zero modulo 4, since they are the 
product of two even numbers. Now we look at the second term and consider whether 
there is some j & A such that Cj = 0. If there is, then Eq. (4.62) implies that 
J2t=i ^fc^fc ^^ even, and since we know that J2keA ^k^k is even that means X^fces ^k^k 
must also be even. This makes the last term zero and makes Eqs. (4.60) and (4.61) 
equal. If Cj = 1 for all j G A, then Ck can be set to 1 in Eq. (4.60) and the two 
results are still the same. 

The last case is when Cj = 1 for all j G B, but there exists a. j & A such that 
Cj = 0. In this case, Eq. (4.62) tells us that J^k^i bk^k is even. There are two ways 
this can happen. Either both Ylk&A bk(^k and XlfeeB bk^k are even, in which case both 
Eq. (4.60) and Eq. (4.61) are always zero, or they are both odd. Assuming Ylik&A ^k^k 
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and J2k(^B^kGk = J2keB^k are both odd, Eqs. (4.60) and (4.61) reduce to 



y^ hek - 1 - X] bk^k - 1 + X] hek - ^ bkCk 

k=l k=l k&A k&A ,. „„\ 

(4.63) 

n n 

and y^hk - I -y^hk - I + y^hk -y^hk, 

k=l k=l fceA feeA 

because the product of an even and an odd number modulo 4 is just that even number 
modulo 4. Now we use that X]fc=i ^fc^fc i^ even and that J2keB ^k^k and Ylik<^A ^k^k 
are both odd to reduce the quantum result further, 



y^ hkek - 1 - X] ^fcCfc - 1 + X] ^fcCfc - ^ fefcCfc 

fc=l A:=l fceA feeA 

n 

= J]fcfcefc-l-l + ^6fcefc-l (4.64) 

fc=i fceA 

= 5]] ^fcCfc + 2 ^ fofcCfc - 3 =4 ^ fofcCfc - 1 = ^ 6fc - 1, 
feG-B fceA fceB fceB 

where once again, in the second to last step, we had the product of an even and an 
odd number giving just the even number. As for the result of our model, we use that 
YlikeB ^k is odd to reduce it to 



fc=l k=l k<^A fceA 



J]^fc + 5Z^^~^^5Z^^ + 5Z^^-~5Z^^ (4.65) 

fceB fceA fceA k<^A keA 

keB VfceA fceA / fceB 



and so our model and quantum mechanics agree in all cases. 



D 



The symmetric difference of complete graphs includes complete graphs them- 
selves. The symmetric difference of two complete graphs and complete bipartite 
graphs actually coincide in the case of star graphs. 
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In all the classes of graphs above, it is possible to get from any node to any other 
node by traversing no more than two edges. So, in a sense a nearest-neighbor model 
communicates through half the graph. This may be why our model works, and so it 
would be interesting to find graphs without this property for which our site-invariant 
nearest-neighbor communication model works. 

4.3 Other communication models 

We have looked at one of the simplest kinds of communication models imaginable, 
this being a site-invariant model that limits communication to nearest neighbors in 
the graph. In this section we will go beyond this. We will consider site-invariant 
models with unlimited communication, models whose communication distance is 
limited but whose decision algorithm can vary from site to site, and even a model 
with no restrictions at all. 

Site-invariant models will always use the LHV table in the last section. The rea- 
son is that allowing the LHV table to break the symmetry of the graph defeats the 
purpose of considering a site-invariant model. We proved earlier that a site-invariant 
model is incapable of reproducing all submeasurement results on a two-dimensional 
cluster state. Motivated by results from measurement-based quantum computa- 
tion [ ■ :], we compare this result to one-dimensional cluster states, i.e., linear chains. 
We find that a site-invariant model is indeed capable of reproducing all submeasure- 
ment results on a linear chain, although more than nearest-neighbor communication 
is needed. 

We also consider models that are not site invariant, and for this we must account 
for an arbitrary underlying LHV table. We analyze such models based on communi- 
cation distance, which we define to be the number of successive edges through which 
information can be sent. For the one-dimensional cluster state, we find that no com- 
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munication distance 1, meaning nearest-neighbor, model is capable of reproducing 
all sub measurement results even if it is not site invariant. In general, we prove that 
any model that reproduces all submeasurement results for all graph states must have 
a communication distance that scales linearly with the size of the graph. 

4.3.1 Site- invariant model for ID cluster states 

In Subsec. 4.2.4 it was shown that for certain graph states, no model based on the 
LHV table of Eq. (4.44) and assisted by a site-invariant communication model is 
capable of reproducing all submeasurement correlations. This result holds even if 
the model allows unlimited distance communication. 

This subsection further explores the symmetry of site invariance by considering it 
in the context of one-dimensional cluster states, showing that linear chains do indeed 
permit successful site-invariant models. The model we describe involves communi- 
cation over a distance equal to the number of edges in the one-dimensional cluster 
state, i.e unlimited communication. At present it is unknown whether the submea- 
surement correlations of a linear chain could be reproduced by a model with limited 
distance communication, although this is unlikely since in the next subsection we 
show that nearest-neighbor communication is insufficient. 

The fact that unlimited communication is allowed is in the same spirit as our 
counterexample of Fig. 4.2, where communication spans the entire graph and the only 
restriction is site invariance. The key simplification in the case of one-dimensional 
cluster states is that all qubits, except those at the ends of the chain, have exactly two 
neighbors. As a consequence, the form of stabilizer elements whose LHV result from 
Eq. (4.44) requires correction is constrained so that the correction can be effected by 
a site-invariant model. 

Before explaining our site-invariant model for linear chains we should learn more 
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about them, introducing terminology as appropriate. For an n-qubit chain, the n 
stabihzer generators are given by gi = X1Z2, Qj = Zj^iXjZj^i for j = 2, . . . , n — 1, 
and Qn = Zn_iXn. Any stabihzer element is a product of generators. An arbitrary 
product of generators can be decomposed into a product of terms each of which 
is a product of successive generators. We call these terms primitive stabilizers or 
just primitives . The primitive stabilizers are separated by the omission of one or 
more generators in the product of generators. An example of a stabilizer element for 
n = 10 qubits is 919293959699 = — ^i-^2^3-^4^5^6^7^8-^9^io- The primitives in this 
example are 919293, 9596, and 5(9. 

Associated with each primitive is a Pauli product, with the sign omitted, for the 
qubits corresponding to the generators in the primitive. We call these Pauli products 
words. For the 10-qubit example above, the words are YiX2Y^, Y^Yq, and Xg. At 
each end of a word, there is an / if one generator is omitted and a Z if two or more 
generators are omitted. We can make these word boundaries apply even at the end 
of the linear chain by embedding our cluster state in an infinite linear chain. The 
generators for the qubits to the left of j = 1 and to the right oi j = n are always 
omitted, and we redefine 91 = Z0X1Z2 and 9n = Zn-iXnZn+i- We can list the entire 
set of words by considering all possible primitives: 

1. X for a primitive with one Pauli operator. 

2. y (g) y for a primitive with two Pauli operators. 

?). Y ® X®*^-'^^^ ® Y for a primitive with j > 3 Pauli operators. 

If a word is bounded by an /, there must be another word immediately on the 
other side of the /. A sentence is a Pauli product consisting of a set of words separated 
by singleton J's and bracketed by Z's at both ends. Words are not stabilizer elements, 
but sentences are. The example above contains two sentences, ^0^1X2^3/4^5^6 ^7; 
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including the zeroth qubit, and ZgXgZiQ. The Z bookends on a sentence separate 
it from other, nonoverlapping sentences in the same overall stabilizer element. Be- 
tween the Z's in successive sentences, there can be an arbitrary number of J's. Any 
stabilizer element is a product of nonoverlapping sentences, and for these elements 
J's occur only between sentences or as singletons between words, X's and F's occur 
only in words, and Z's occur only as the boundaries of sentences. 

Recall that the goal of the communication model is to introduce a sign flip into 
the product of LHV table entries for those Pauli products that are the negative of a 
stabilizer element. The only words that introduce a minus sign into the corresponding 
product of generators are those of the form Y X®*^-^~^^ Y with j odd. Thus a 
candidate for a site-invariant communication model is the following. 

1. Each site at which an X or a Z is measured broadcasts the measurement 
performed upon it. 

2. Each site that measures X determines if it is the middle, implying an odd 
number of X's, qubit in a word of the form Y0X'^^^~'^^0Y in a submeasurement 
sentence, and if so, flips its LHV entry, i.e., changes (—1)^ to — (— l)"*. 

This clearly gets any stabilizer right and thus all global measurement correlations 
right. 

The only question remaining is whether this model works for submeasurements. 
We answer this question by showing the following: two sentences, 5*1 and 5*2, that are 
submeasurements of the same global measurement, generally not a stabilizer element, 
must be identical on the region where they overlap, except possibly at bracketing 
Z's. This property implies that 5*1 and 5*2 have exactly the same words in the region 
of overlap. Thus, for any pair of submeasurements of the global measurement, a 
sign flip arising from a word of the form Y X"^^^^"^^ 0Y in the overlap region is 
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common to both submeasurements. Since both the word and the sign flip occur 
in both submeasurements, our model correctly predicts both outcomes. Therefore, 
proving this property will complete the proof that our site-invariant model correctly 
predicts all submeasurements. 

To prove this property, notice flrst that if 5*1 and 5*2 overlap, there are two cases: 
the region of overlap coincides with one of the sentences, or it does not. In the former 
case, we choose 5*2 to be the sentence that coincides with the region of overlap, and 
in the latter case, we choose 5*1 to be the sentence on the left and 5*2 to be the one on 
the right. With these conventions, the left boundary of the overlap region coincides 
with the Z that bounds the left end of 5*2, and the right boundary of the overlap 
region coincides in the former case with the Z that bounds the right end of 5*2 and 
in the latter case it coincides with the Z that bounds the right end of 5*1. 

To be submeasurements of the same global measurement, the two sentences must 
satisfy the following basic rule: in the overlap region, sites within a word of one 
sentence must be occupied in the other sentence by the same Pauli operator or by 
an /. Since Z's do not occur in words, this rule implies that the Z's that bound the 
overlap region at either end in one of the two sentences cannot occupy a site within 
a word in the other sentence and thus must be a bounding Z or a singleton I in 
the other sentence. The submeasurement requirement, by itself, implies that in the 
overlap region, the site of a singleton I in one sentence can be occupied by anything 
in the other sentence, but the available words impose a much stronger constraint, as 
we now show. 

Consider the left boundary of the overlap region, which is occupied by the leftmost 
Z in 5*2 and by a Z or a singleton / in 5*1. Immediately to the right in both 5*1 and 5*2 
is a word. When one of these words is shorter than the other, the basic rule implies 
that the shorter word must be a preflx of the longer one. A glance at the allowed 
words 1-3 shows, however, that none is a preflx of another. Thus 5*1 and 5*2 must 
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have the same word in this first overlap position, which is followed by a singleton / in 
both sentences. Applying the same logic to this and subsequent singleton J's shows, 
as promised, that 5*1 and 5*2 are identical in the overlap region, except possibly at 
the boundaries. 



4.3.2 Nearest-neighbor models for ID cluster states 

Two properties that make for a simple communication model are site invariance 
and nearest-neighbor communication. We showed that a site-invariant model can 
be made to work for linear chains if unlimited communication is allowed. In this 
subsection we prove that nearest-neighbor communication models are incapable of 
reproducing all submeasurement correlations for a linear chain. 

In this subsection, for the first time, we need to consider communication models 
that are not site invariant. This means that we also need to incorporate arbitrary 
LHV tables on which the model is based. We do this by considering the outcome, 
according to a nearest-neighbor communication model, of a measurement of a Pauli 
matrix a on qubit j. The output of this measurement in such a model depends on 
both the qubit in question and the measurements made on neighboring qubits. We 
write this in the form a," where j is the qubit being measured, s is the lower case 
version of the Pauli operator measured upon it, and a consists of the Pauli matrices 
being measured on the neighboring qubits. For instance, xj^ is the output for the 
model of an X measured on the jth qubit when the left neighbor measures Y and 
the right neighbor measures Z, and z^ is the outcome of a Z measurement on a 
qubit j who is at the end of the chain and whose neighbor chooses to measure X. 

We construct our counterexample on a chain of 11 qubits. By considering ap- 
propriate measurements, we arrive at constraints on the communication model that 
yield a contradiction. It turns out that we can get a contradiction after considering 
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Measurement Outcome of the Model 

ZiX2YsX4^X^,XqZ7 Z^ X2 X^ Xq Z-j =1 

Z11X2I •;iX^Y -^J\q/jj Z^ X2 3^4 Xq Z^ ^1 

7 Y Y Y Y Y 7 ^IX ^ZX^XX^XZ ^XI _ ^ 

Z2A3A1A5A6A7Z8 Z2 X3 X5 X-j Zg 

/ _ 

^8 — ~ 



7VYVYVYV7 ^lY „,ZX ^YY „,XX .,XX ^YY „,XZ ^Y I 

Z2A3A4r5A6l7^8 2:2 0:3 J/g a;g ?/7 Zg 



7 Y Y Y Y Y 7 ^IX „ZX ^XY JYZ ^XI _ 1 

Z3A4A5A6J rAgZg 2:3 0:4 Xg Xg Zg — i 

7 Y Y Y 7 ^IX^ZX^XZ^XI _ 1 

Z3A4A5AgZ7 Z^ X^ X^ Z-j — i 

7 Y ^' Y ^' Y 7 ^IX„ZY„YY„YZ^XI _ 1 

Z3A4J -jAgi vAgZg Z^ X^ Xg Xg 2:9 — i 

^3X4^ ',X,Z, zi^x^x^z^' = 1 

z,x,x,y,x,Y,z,, f/iyx^^z^'h^''''' ^ "^ 

Z/^Xc^Ji(iXjZg z^ x^ Xj Zg ^1 

7 Y ^' Y ^' Y 7 ^IX^ZY^YY^YZ^X _ 1 

^sAeJ vAgJ 9A 10^11 2:5 Xg Xg x^Q z^-^^ — i 

7 Y V Y 7 ^IX„ZY„YZ^XI _ 1 

^sAeiTAsZg 2:5 Xg Xg Zg — i 

7 IT \" IT \" "T 7 ^IX„ZY^YY^YZ^X _ 1 

^sAeivAsioA 10^11 2:5 Xq Xq a;;^Q 2;]^;^ — i 

7 "T V "T 7 ^/x™zy„yz^x/ _ 1 

^SAeJ^TAgZg 2:5 Xg Xg 2;g — i 

Table 4.1: In the left column is listed Pauli measurements on a linear chain with 
11 qubits. For each global measurement, a definite submeasurement is shown in 
black. These definite submeasurements lead to the contraints on the outcome of 
a nearest-neighbor communication model listed in the right column. If these 16 
equations in the right column are carefully multiplied together, all the outcomes will 
cancel and what will remain is the equation 1 = — 1. Since this is impossible, no 
nearest-neighbor communication model can reproduce all 16 of these correlations. 
Note that this argument is independent of the underlying LHV table and assumes 
nothing about the details of the communication strategy other than communication 
being only between nearest neighbors. 



the 16 measurements listed in Table 4.1. This proves that no nearest-neighbor com- 
munication model can reproduce all submeasurement correlations arising from Pauli 
measurements on a linear chain. This result is independent of the underlying LHV 
table and does not assume a site-invariant model. 

We can most likely generalize this result to say that no limited distance commu- 
nication model can reproduce all submeasurements on a linear chain. The reasoning 
would be similar to the strategy of the next section. However, we have not yet done 
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this and so the question of whether a hmited distance communication model can 
reproduce Pauh measurement correlations on a ID cluster state still remains open. 

4.3.3 Non-nearest-neighbor communication models 

We have shown that site-invariant communication models are incapable, for at least 
some measurements, of reproducing all submeasurement correlations. Therefore any 
successful communication model must reject this property. In this section we show 
that any model whose communication distance is constant or scales less than linearly 
with the number of qubits fails as well to predict some submeasurements correctly. 
As a result we have that any model reproducing all quantum correlations must have 
communication that spans nearly the entire graph, as well as not being site invariant. 

In a communication model with communication distance d, nodes j and k can 
signal to each other if there exists within the graph a path from j to k that traverses d 
or fewer edges [■:>(>]. Put another way, this is the statement that information can only 
be transmitted along edges and that the number of successive edges through which 
some piece of information can be sent is at most d. In this section we prove, via con- 
tradiction, that no communication-assisted LHV model for which the communication 
distance satisfies 



d<4 



n 1 
24 ~ 2 



+ 1 (4.66) 



correctly reproduces the predictions of quantum mechanics for all submeasurements 
on all graph states of n qubits. 

The proof relies on an infinite class of graph states for which a set of five global 
measurements can be chosen that are not locally distinguishable. Each of these global 
measurements includes a submeasurement that can be written in terms of stabilizer 
elements and is thus definite. The output of each qubit, however, must be such 
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that the correct values are obtained for all submeasurements that are consistent 
with its observable surroundings. This requirement, for the particular states and 
measurements chosen, yields a contradiction. 

To begin, consider the graph state corresponding to an n node ring where n = 12 f 
and / is an odd positive integer. Let the qubits be numbered sequentially, starting 
with 1 at an arbitrary point on the ring and moving clockwise along it. Additionally, 
define the following subsets of the n labels: 

V = {4/,8/,12/}, 

M = {2f,6f,10f}, 

3^ = {j|j = lmod2}, 

(4.67) 
£ = {j\j ^V,M and j = 2 mod 4} , 

n = {j\j ^V,M and J = mod 4} , 

Sk = {j\2f{k-l)<j<2fk}. 

For our purposes, it is useful to think of the ring as arranged in an equilateral triangle 
with vertices specified by the subset V as shown in figure 4.3. The midpoints of the 
legs of the triangle are then given by the subset Ai, and the segments between 
adjacent vertices and midpoints are given by the S/s. We use the notation Sj^k as 
shorthand for Sj U Sk, and A\B is used to denote the set consisting of the elements 
of A that are not in B. 

Now consider global measurements of the form ^jJ^ Mj, with 



M, 



X or r if j G V, 

Y iijey, (4.68) 

X otherwise, 



such that the number of vertices measuring Y is not one. These global measurements 
include the following submeasurements, for which quantum mechanics predicts an 
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outcome with certainty. 



For M4f = X, Msf = X, M12/ = X, 

6/ 

X2fX4^fXefXsfXiofXi2f II Xj = ||5f2j, (4.69) 

jecun j=i 

implying a measurement outcome of +1. 



For M4/ = Y, Msf = Y, M^y = Y, 



(4.70) 



(4.71) 



X2/X6/X10/ 1 1 ^i 1 1 -^fc 
jey kec 

3/ 
= [[[-94j-3g4.j-2g4.j-lJ , 

implying a measurement outcome of (— 1)^'^ = —1. 



For M4/ = Y, Msf = X, M^f = Y, 

Y4fX6fXsfXiofYi2f II ^i II Xk 
jeyr\Si,2 keTZuic\Si-2) 

/-I 6/ 

= giYl[-g4j-ig4jg4j+ijg4f-i JJ g2k , 

j=l k=2f 

implying a measurement outcome of (— 1)-^^^ = +1. 



Cyclic permutation of this last measurement yields two more with +1 outcomes. 



For M4/ = Y, Mgf = Y, M^j = X, we have 

X2fY4fYsfXiQfXi2f II Yj II Xk , (4.72) 

j&ynS3,4 kenuic\S3,4) 
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Figure 4.3: Example demonstrating that any communication-assisted LHV model 
with communication distance d < 1 fails to reproduce some submeasurements for 
the ring with n = 12 nodes. Five of the global measurements shown on the ring, 

YiX2Yo^X4Y^XQYTXsYc)Xif)YiiXi2 , 

YiX2YsYiY^XeY7Ys,YgXiQYiiYi2 , 

YiX2Y3Y4Y^XeYTXsY;jXiQYiiYi2 , 

YiX2YpX4Y5X(iYjYgYc)XioYiiXi2 , 

submeasurements shown in black that are useful for ex- 

a contradiction. These submeasurements imply the fol- 

on a nearest-neighbor communication model, 

ry <-y» rif /y /y /y \ 

yXx2ylylx^y^ylx^QyX^ = -1 , 

Y Y 1 

Y Y 1 

•^21/41/5 1/7 Vs-^io-^u ~ -*- ) 

Y Y 1 

•^2-^4-^6^8^9 ^11^12 ~ -*- 5 



contain 
hibiting 
lowing 



constraints 



which when multiplied together yield a contradiction. 
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and for M4/ = X, Mgj = Y, Muf = F, we have 



X,fX,fX,jYsfY^y n ^^ n ^^- (4.73) 

jeyn55,6 k€TZu{c\S5,(i) 



Now assume there exists a distance d = 2f — 1 communication-assisted LHV 
model that correctly replicates the predictions of quantum mechanics for all Pauli 
measurements on n qubits. As in the previous subsection, the output of such a model 
can be fully described in terms of single-qubit LHV entries whose value depends both 
on the qubit in question and on the measurements made by other qubits within its 
communication range. We write these hidden variables in the form Sj", where in 
this case a indicates the measurements made on all qubits within its communication 
range. The global measurements utilized for Eqs. (4.69) through (4.73) have the 
virtue that the communication range of each qubit includes at most one other qubit 
whose measurement is changeable, and that is the qubit at the nearest vertex. Thus, 
in comparisons between them, the measurement performed on, at most, a single 
qubit need be included in a. Moreover, the qubits at the center of each side of the 
triangle cannot see the changes at the vertices. Consequently, the constraints implied 
by Eqs. (4.69) through (4.73) on a hidden variable model with communication range 
d can be expressed as follows: 
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1 —X2fX4fXQfXsfXiQfXi2f J_J_ Xj , 
-1 =X2fX6fXiof YlyJYlxl , 

1 =y'ifX6fXsfXiofyi2f Yl yJ 
jeyr\Si,2 

n -? n 

ke{cuTZ)nS4,5 i&{cnS3,6)u{ii\S4,5) 

1 =x2fyify8fXiofXi2f Yl yJ 

i6yn53,4 

n -? n 

fcG(£u7^)n56,l ^e(£n55,2)u(7^\56,l) 






(4.74) 



X 



Y 

I ; 



1 =X2fXAfX(ifysfyi2f JJ ^. 



y 

ieyn55,6 



n -? n 

A:e(£u7^)n52,3 ie(£n5i,4)u(7e\52,3) 



xr. 



Using the identity ^ = ^ n (5i U ^2 U ^3 U ^4 U .Ss U Sq) for A = y, C, otTZ 
and the fact that all variables square to 1, it can be shown that the right hand 
side of the first equation in Eq. (4.74) is equal to the product of the right hand 
sides of the other four equations. Thus, we have the contradiction 1 = — 1, showing 
that no distance d communication-assisted LHV model reproduces the predictions of 
quantum mechanics in this instance. 

For other values of ri 7^ 12/, with / odd, an identical contradiction applies to a 
graph consisting of r = {n — 12) (mod 24) < 24 unconnected nodes and a ring of 
size n — r. So for any value of n we have a ring of size n — r = 12{2k + 1) for some 
integer k, and we have proven that in order to reproduce all quantum correlations 
we must have 



d>2{2k + l) -1 =4A; + 1 = 4 



n-12 
24 
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as desired. 



4.3.4 Universal model 



To end this chapter, we lay out a communication model that actually does reproduce 
all quantum correlations for Pauli measurements on all graph states. We know that 
this model cannot be site invariant and that it must also have a communication 
distance that scales linearly with the number of qubits. For these reasons we impart 
our model with two properties. One property, motivated by our understanding of 
site invariance, is that there is a special node capable of breaking symmetry. The 
second property, motivated by what we know about communication distance, is that 
we allow communication to span the entire graph. 

Our universal model involves two rounds of communication. During the first 
round, everyone sends their measurement choice to a single party. During the second 
round of communication, this special party sends a classical bit to each qubit needing 
to change their LHV output. The algorithm allowing this special party to decide 
which qubits must change their LHV result is explained below. 

There are 2" possibly non-unique submeasurements, for an n-qubit graph state, 
which need to be correct in the end. These correspond to the 2" binary vectors e 
that specify which measurement results of the global measurement, M, to include in 
the submeasurement. However, we need not worry about random submeasurements 
and only focus on those definite submeasurements satisfying Eq. (4.42). In fact, 
Eq. (4.42) can be used in a simply way to find all definite submeasurements. To see 
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this, rewrite the condition to be a definite submeasurement as 

n 

[r2(M)],e, =2 J]r,fc[ri(M)]fcefc 

k=l 

n 



fc=i 



^ re' =2 0, 

where we have defined Tjk = Tjk[ri{M)]k — [r2{M)]j6jk. We also have, by multiplying 
these equations by [ri{M)]j, that definite submeasurements must satisfy 

n 

^([ri(M)],r,fc[ri(M)], - [ri(M)],[r2(M)],5,,)e, =, 



k=l 



^ T'e^=,0, 



(4.77) 



where F^^ = [ri{M)]jTjk[ri{M)]k — [ri{M)]j[r2{M)]j6jk is a symmetrized version of 

By rewriting the condition to be a definite submeasurement in this way what we 
find is that definite submeasurements are binary vectors in the modulo 2 nuUspace 
of some matrix F. This means that definite submeasurements form a subspace. If 
we can arrange so that all definite submeasurements corresponding to a basis of this 
nuUspace are correct, then we claim that all definite submeasurements will be correct. 
For, consider two binary vectors in a basis for the modulo 2 nuUspace of F, e and 
/. The quantum mechanical prediction for one of these submeasurements, say the 
one corresponding to e, is denoted by (— 1)**" and can be written from Eqs. (4.42) 
and (4.40) as i to the power 



Y,[n{M)],e,[Tjk[ri{M)]kek - [r2(M)],-e,-^ 

n n 

J2e,Y,{[n{M)]jT,k[ri{M)]k - [n(M)],[r2(M)],-5,-,)efc (4.78) 



i=i fc=i 



: J2 e.q.e, = ere^. 
j,k=i 
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From this, we find that the quantum mechanical outcome for the submeasurement 



given by e + / is 2 to the power 



T 



e + / P e + / . (4.79) 



We want to know the value of Eq. (4.79) modulo 4. We can do this with the help 
of a basic fact of modular arithmetic. That is, a coefficient times an even number 
modulo 4 is the same as that coefficient modulo 2 times the even number, i.e., 

c • (even number) =4 c ■ (even number). (4.80) 



The reason this is relevant is because, since e + / describes a definite submeasure- 
ment, the entries of T%[e + f) are even by Eq. (4.77). Hence we have that this 
submeasurement result is i to the power 

(e + /)r(^T7)'' = er(^T7)^ + /r(^T7)^. (4.8i) 

Now, once again, eF^ has even entries since e corresponds to a definite submeasure- 
ment and F^^ = Ff-. Therefore the exponent of i is really, 

er(e + fY + frfe + fY = ePe^ + ft'f + e(r + irf)f 

^ ^ ^ ^ (4.82) 

= Se + Sf + 2eFV^ =4 Se + Sf, 

where we used the even entries of T^f^ to say that 2eF''/^ =4 0. Thus if our model 
correctly predicts Se and Sf, it also correctly predicts Se+f- 

So it only remains to fix the LHV table for those submeasurements corresponding 
to a basis for the modulo 2 nuUspace of F. Any basis will do, so we choose a basis with 
the property that each basis vector contains a measurement that is not contained 
in any other basis vector. Such a basis exists and can be obtained by writing the 
vectors of any basis as the rows of a matrix, and then row reducing that matrix to 
obtain a basis with the desired property. Having such a basis is nice because the 
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submeasurement corresponding to each basis vector can be corrected, if need be, by 
having the measurement unique to that basis vector change its LHV result. 

We now have a universal communication model for Pauli measurements on sta- 
bilizer states. First, each party should communicate their Pauli measurement to 
a special party. This party then calculates a basis for the modulo 2 nullspace of 
r as described above. For each basis vector whose corresponding submeasurement 
disagrees with the quantum prediction, this special party sends one bit of classical 
information to a party who can fix it. With this done, all sub measurements are 
guaranteed to agree with quantum mechanics. 

4.4 Concluding remarks 

This chapter explores the degree of nonlocality present in stabilizer states by consid- 
ering communication-assisted LHV models. These models can be classified into two 
types, depending on whether or not they are site invariant, and can be parameterized 
by the allowed distance of communication, where the distance between two qubits 
is defined as the minimal number of edges between the corresponding nodes in the 
stabilizer-state graph. Interestingly, a simple nearest-neighbor site-invariant com- 
munication model is capable of yielding the global quantum mechanical correlation 
for any measurement of Pauli products on any stabilizer state, but the submeasure- 
ments of these global measurements are much harder to reproduce. To replicate the 
predictions of quantum mechanics for all submeasurements, it is necessary to violate 
site invariance and the communication distance must scale as n/6 or faster in the 
number n of qubits. 

We also looked at specific graphs of interest in measurement-based quantum com- 
putation. In particular, for the two-dimensional cluster state, site invariance must be 
broken, while for the one-dimensional cluster state a successful site-invariant model 
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exists provided communication span the entire graph. Nearest-neighbor communi- 
cation is insufficient to reproduce all submeasurement correlations for both types 
of graphs, suggesting that site invariance, rather than communication distance, is 
a more relevant quantity in the context of measurement-based quantum computa- 
tion. Indeed, a site-invariant model can be made to reproduce all submeasurement 
correlations in a Pauli product measurement for the one-dimensional cluster state, 
but not for the two-dimensional cluster state. Finally, a successful model exists for 
all submeasurements if the restrictions of both site invariance and limited distance 
communication are relaxed. 

The obvious next step with communication models is to go beyond Pauli mea- 
surements on stabilizer states. Including one additional measurement is problematic 
because probability-preserving tables do exist, and so it is likely that any communi- 
cation strategy will have to be probabilistic in nature. An alternative is to consider 
Pauli measurements on a non-stabilizer state for which Pauli measurements are uni- 
versal. The ffist step there would be to see if probability-preserving tables exist for 
Pauli measurements on that state. Keeping in mind the results for non-Pauli mea- 
surements on stabilizer states, it seems unlikely that probability-preserving tables 
exist for Pauli measurements on these non-stabilizer states. However, if they do ex- 
ist, then examining communication strategies for such correlations would hopefully 
result in a plethora of new and interesting results. If probability-preserving tables do 
not exist in either case, then we must begin looking at probabilistic communication 
models. 
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Conclusion 



Fully harnessing the information processing power inherent in quantum systems re- 
quires an understanding of entanglement and its role in quantum information pro- 
cessing tasks. One perspective on the role entanglement plays in empowering quan- 
tum computation derives from measurement-based quantum computation, a unique 
method by which universal quantum computation can be achieved. A significant 
implication of this scheme is that the entanglement resource required for the compu- 
tation is independent of the particular algorithm being implemented; every algorithm 
begins with the same fiducial state, usually a stabilizer state, so that algorithms dif- 
fer only in the particular measurements made. The measurement outcomes from 
the computation are correlated in such a way that is impossible to reproduce with 
classical correlations alone. 

The central theme in this work is the study of correlations, which are usually 
imagined as resulting from projective measurements on quantum states. According 
to measurement-based quantum computation, these quantum correlations somehow 
capture all the computational power a quantum computer can possess. The difficulty 
now is to understand how these correlations relate to computational power. We take 
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the point of view that one should "subtract off" classical correlations by modeling 
measurement results with necessarily imperfect local-hidden-variable tables. The 
focus then shifts to how one should quantify the difference between classical and 
quantum correlations based off the performance of these LHV tables. There are 
indeed many ways one could do this. 

We quantify the difference between classical and quantum correlations by supple- 
menting LHV tables with classical communication and analyzing the success or fail- 
ure of a communication model with regards to various properties of the model. The 
idea of using classical communication to help understand entanglement properties of 
quantum states is a familiar one [37, 69]. In fact, one perspective is that quantum 
mechanics is a nonlocal theory; that is, measurement results are predetermined, but 
measurements on one system can effect another system instantaneously. The role of 
classical communication is to make up for this "spooky action at a distance." 

The idea of using classical communication and shared randomness to reproduce 
quantum correlations goes back, in the bipartite scenario, to Toner and Bacon [92], 
and, in the multipartite scenario, to Tessier et al. [9U]. One big difference between 
these two cases arises from the lack of submeasurements in bipartite systems. Re- 
sults in this work suggest that submeasurements give rise to correlations that are 
much more difficult to reproduce than global measurements alone. In the bipar- 
tite scenario all measurement correlations can be reproduced with just one bit of 
communication. This includes not only Pauli measurements, but all projective mea- 
surements. Therefore, in the attempt to generalize our results beyond Pauli measure- 
ments, it is reasonable to suppose that a similarly simple strategy may work for global 
measurements on multipartite states. Indeed, extending our simple communication 
strategy beyond Pauli measurements would strengthen the claim that submeasure- 
ments should receive more attention than global measurements. Submeasurement 
correlations, then, would hold the key to gaining more understanding of the power 
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of measurement-based quantum computation. 

Perhaps a reasonable approach to extending our results beyond Pauli measure- 
ments is to use Bell inequalities with auxiliary communication [6] . This is essentially 
the main idea behind our proofs in Subsecs. 4.3.2 and 4.3.3. Extending the results 
in Chapter 3 could help in this endeavor. Since one expects the communication 
required to reproduce quantum correlations for all submeasurements to grow expo- 
nentially, Bell inequalities with auxiliary communication seem more likely than a 
direct approach to yield interesting results. Perhaps Bell inequalities with auxiliary 
communication could even lead to an entanglement measure capable of capturing 
the computational power inherent in a quantum state. 
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Appendix A 

Graphical description of 
controlled-Z gates 



The transformation rules given in Subsec. 2.3.5 suffice to describe the effect of any 
local Clifford operation on a stabilizer state. In order to complete the set of trans- 
formation rules for the Clifford group, we include transformation rules for the '"Z 
gate here. In the interest of brevity, we consider only transformations of reduced 
stabilizer graphs. Transformation rules for general stabilizer graphs can be derived 
by ffist using equivalence rules El and E2 in Sec. 2.3.6 to convert to an equivalent 
reduced graph and then applying the transformation rules below. 

T(viii). Applying Z between two solid nodes complements the edge between them. 

T(ix). Applying '^Z between a hollow node and a solid node complements the edges 
between the solid node and the neighbors of the hollow node. 

Flip the solid node's sign if the two nodes were initially connected and the 
hollow node did not have a sign or if the two nodes were not connected and 
the hollow did have a sign. 
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T(x). Applying Z between two hollow application nodes complements an edge pro- 
vided its endpoints neighbor distinct application nodes. 

Nodes that neighbor both application nodes flip their signs. If an applica- 
tion node initially had a sign, flip the signs of nodes connected to the other 
application node. 

Transformation rule T(viii) is trivial since the Z gate simply commutes into 
layer two of the reduced-graph-form circuit. The circuit identities needed to prove 
rules T(ix) and T(x) are given in Fig. A.l. 
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Figure A.l: Circuit identities for transforming a reduced stabilizer graph under Z 
gates. In (a)-(c) the identity is illustrated for the case of four neighbors connected 
to a qubit of interest; equivalent expressions hold for other numbers of neighbors. 
Identities (a) and (b) give the rule T(ix). The effect of a Z gate on the hollow qubit 
in the third layer of (a) or (b) is to deposit an additional Z gate on both the hollow 
and solid qubits in the third layer of the final circuit. The identity in (c) relies on (b); 
this identity is extended in (d) to a demonstration of the rule T(x). If there is a Z 
gate on the lower hollow qubit in the third layer of (c), the result in the third layer 
of the final circuit is to put a Z gate on that qubit and on the neighbors of the other 
hollow qubit. Translated to (d), this means that a Z gate in the third layer on either 
hollow qubit leads in the third layer of the final circuit to additional Z gates on that 
qubit and on the neighbors of the other hollow qubit. 
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Proof that the equivalence rules 
are complete 



Having described a set of rules in Subsec. 2.3.6 for converting between equivalent 
stabilizer graphs, we show in this appendix that, in each case, the aforementioned 
rules generate the entire equivalence class of stabilizer graphs. The proof is divided 
into three parts. The first part shows how to use rules El and E2 to convert an 
arbitrary stabilizer graph to an equivalent graph in reduced form. The second part 
explains how testing the equivalence of a pair of reduced stabilizer graphs can be 
simplified, using rules E(i) and E(ii), to a special form. Finally, the third part proves 
that the graphs on the two sides of such a simplified equivalence test are equivalent 
only if they are trivially identical. Taken as a whole, this proves that the set of 
equivalence rules given in Subsec. 2.3.6 is sufficient to convert, reversibly, between 
any two equivalent stabilizer graphs and thus that they generate all stabilizer graphs 
that are equivalent to the same stabilizer state. Similarly, considering only the final 
two parts of the proof shows that rules E(i) and E(ii) are sufficient to generate all 
reduced stabilizer graphs. 
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B.l Converting stabilizer graphs to reduced form 

Two features identify a stabilizer graph as reduced. In a reduced graph, hollow nodes 
never have loops, and hollow nodes are never connected to each other. Equivalence 
rule El can be used to convert looped nodes from hollow to solid. Applying rule El in 
this sort of situation can cause hollow nodes to acquire loops, but each application 
fills one hollow node of the graph, so the procedure will terminate in at most a 
number of repetitions equal to the number of hollow nodes in the graph. Similarly, 
connected hollow nodes in the resulting graph can be converted to solid nodes using 
the appropriate case of rule E2. Once again, the process is guaranteed to terminate 
because the number of hollow nodes to which the rule might be applied decreases by 
two with each application of the rule. Concomitantly, the conversion of a stabilizer 
graph to an equivalent reduced graph never increases the number of hollow nodes. 

B.2 Simplifying reduced-graph equivalence 
testing 

Equivalence testing for pairs of reduced graphs is facilitated by simplifying the equiv- 
alence such that nodes that are hollow only in the first graph never connect to nodes 
that are hollow only in the second. This simplification can be accomplished by iter- 
ating the following process. Choose a pair of nodes, a and b, that are connected in 
either graph and such that a is hollow in one graph and b is hollow in the other. To 
the graph in which they are connected, apply the relevant reduced equivalence rule 
to the selected nodes. Among other things, the equivalence operation reverses the 
fill of the two nodes it is applied to. Since one node is hollow and the other solid, 
this preserves the total number of hollow nodes while yielding a node that is hollow 
in both graphs. Because it is applied only to unpaired hollow nodes, subsequent uses 
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of this equivalence rule do not disturb the newly paired hollow node. Consequently, 
this process also terminates in at most a number of repetitions equal to the number 
of hollow nodes in each of the graphs. 

B.3 Trivial evaluation of simplified reduced-graph 
equivalence tests 

The two reduced graphs composing a simplified reduced-graph equivalence test are 
equivalent, i.e., correspond to the same state, if and only if the graphs are identical. 
To see why this is so, we return to graph-form quantum circuits. In addition to the 
standard restrictions for reduced graphs, the circuits corresponding to the graphs in 
a simplified equivalence test have the following property: if in one of the circuits, 
a qubit with a terminal H participates in a, Z gate with a second qubit, which 
cannot have a terminal H, then in the other circuit, it cannot be true that the 
second qubit has a terminal H and the first does not. We prove the triviality of 
simplified reduced-graph equivalence testing by considering an arbitrary simplified 
reduced-graph equivalence test and showing that the two graphs must be identical 
if they are to correspond to the same state. 

In terms of unitaries, an arbitrary graph-form circuit equality can be written as 

geHi jeSi heZi yeCi k m 1 ^ 

= n ^^ n ^ n ^'^ n %n^'^io)''" 

geHr j&Sr heZr iGCr k 

where C lists the pairs of qubits participating in Z gates and Ti., Z, and S are 
sets enumerating the qubits to which H, Z, and S gates are applied respectively. 
The total number of qubits is denoted by n and the subscripts / and r discriminate 
between the circuits on the left and right hand sides of the equation. In terms of 
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these sets, a reduced-graph-form circuit satisfies the constraints 7i Ci S = and 
{a, b} ^ C for all a,b ^ H. The circuits in simplified tests also satisfy {a, b} ^ Q, Cr 
for all a EHinHr and 6 G 7i/ fl Ti^, where H denotes the complement of H, i.e., the 
set of qubits to which H is not applied. 

Suppose now that the two graphs have hollow nodes at different locations; i.e., 
at least one of the sets. Hi fl Hr and HinHr, is not empty. For specificity, let's say 
that Hi n Hr is not empty. In the language of circuits, this means that there exists 
a qubit a that has a terminal H on the left side of Eq. (B.l), but not on the right 
side. Since qubit a is part of a circuit for a reduced graph, it does not participate in 
Z gates with other qubits that possess terminal H gates. Consequently, on the left 
side of Eq. (B.l), the Z gates involving qubit a can all be moved to the end of the 
circuit where they become X gates with a as the target. Doing this and transferring 
the '"X gates to the other side yields. 



geHi jeSi heZi jeCi s.t. a^7 k 

h&Ni(a) geHr j€Sr heZ,- 76Cr k 



where A/; (a) denotes the set of qubits that participate in '"Z gates with qubit a on 
the left hand side of Eq. (B.l). 

Because the original graph equality, Eq. (B.l), was simplified, J\fi{a) HHr = 
and, by assumption, a ^ Hr, so the X gates and the terminal Hadamards on the 
right side of Eq. (B.2) do not act on the same qubits. Thus we can commute the X 
gates past the terminal Hadamards. Moreover, we can then move the X gates to 
the beginning of the circuit where they have no effect and can therefore be dropped. 
During this migration, however, they generate a complicated menagerie of phases. 
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The resulting expression for the right side of Eq. (B.2) is 

(\lSr(") / \lZr(") 

/ heZr \feM{a) 



11% n % n ^^ n^'^io)' 

7eC,. 5eAfir{a) ceAfi{a)rWr{a) k 



(B.3) 



where Mr {a) is defined similarly to J\fi{a), 

Mria) = {{p,q}\p e M{a),q e Kia),p ^ q} , (B.4) 

and 1 represents an indicator function, i.e., l^^(a) equals 1 if a G ^^ and otherwise. 

It might appear that we have made things substantially worse by this rearrange- 
ment, but in an important sense, Eq. (B.3) is now very simple with regard to qubit a: 
H is applied to qubit a followed by a sequence of unitary gates all of which are di- 
agonal in the standard basis. This implies that there are an equal number of terms 
in the resultant state where qubit a is in the state |0) and the state |1). On the left 
side of Eq. (B.2), however, the only gate remaining on qubit a is the identity or an 
X, depending on whether a E Zi; thus qubit a is separable and is either in state 
|0) or |1) depending on whether a & Zi. Consequently, our initial assumption that 
Til 7^ TYj. is incompatible with satisfying the equality. 

The preceding discussion shows that two graphs composing a simplified equiva- 
lence test are equivalent only if they have hollow nodes in exactly the same locations. 
Given this constraint, the terminal Hadamards can be canceled from both sides of 
Eq. (B.l), giving 

h&Zi jdSi -fGCi k 

h^Zr j&Sr "l&Cr k 

The state after the initial Hadamards is an equally weighted superposition of all the 
states in the standard basis. The subsequent unitaries are diagonal in the standard 
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Figure B.l: An example of the process described in Appendix B.3. A simplified 
reduced-graph equivalence test is shown in (a). The test is simplified because a node 
that is hollow in only one graph is never connected to a node that is hollow only in 
the other. The test is also false, since such an equivalence is satisfied if and only if 
the two graphs are identical. By switching to circuit notation, it is possible to see 
why these graphs correspond to different states. Translating the graph equality in (a) 
into circuit notation yields the circuit equality in (b). The fact that the graphs are 
reduced implies that qubits with terminal H gates, such as qubit 3 in the left hand 
circuit, are not connected to other qubits with terminal H gates. This allows us to 
pull out the '^Z gates acting on qubit 3 on the left side and transfer the resulting 
'"X gates to the other side of the equation, yielding the equality in (c). That the 
graph equality was simplified guarantees that H gates do not prevent us from then 
pushing the X gates to the beginning of the right hand circuit in (c), as shown 
in (d). Though this is not generally the case, qubit 3 winds up separable on both 
sides of the equation, allowing us to verify the inequivalence of the prepared states 
since |0) ^ (|0) -i|l))/v^. 



basis, so they put various phases in front of the terms in the equal superposition. 
Since a unitary is fully described by its action on a complete set of basis states. 
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Figure B.2: An example sequence of graphs representing parts 1 through 3 of the 
equivalence checking procedure in Appendix B. The graph on the left hand side of (a) 
is equivalent to the reduced graph on the left hand side of (b) by the application of E2 
to the pair of nodes {1,5}. Likewise, an application of El to node 3 in the graph 
on the right hand side of (a) yields the reduced graph on the right hand side of (b). 
The graph of (c) results from the application of E(ii) to the pair of nodes {1,2} in 
the graph on the left of (b). An application of E(i) to the pair of nodes {3,4} in 
the graph on the right hand side of (b) also yields the graph in (c), verifying, as per 
Appendix B.3, that the graphs in (a) represent the same state. 

demanding equality term by term in Eq. (B.5) amounts to requiring that 



n^^n^n%=n^^n^ n 



Cr, 



^7 ' 



(B.6) 



h£Zi 



jeSi jeCi heZr jeSr 76C,. 

which is only satisfied when Zi = Zr, Si = Sr, and Ci = Cr- Thus, after simplification, 
the equivalence of pairs of reduced graphs is trivial to evaluate, since equivalence 
requires that the two graphs be identical. 

An example which illustrates the circuit manipulations described algebraically 
in the text of this section is given in Fig. B.l. An example of the entire process of 
testing graph equivalence is given in Fig. B.2. 

As mentioned above, our proof provides a constructive procedure for testing the 
equivalence of stabilizer graphs. Moreover, it shows that equivalence rules El and 
E2 given in Subsec. 2.3.6 are sufficient to convert between any equivalent stabilizer 
graphs, and that the rules, E(i) and E(ii), for reduced graphs are sufficient to convert 
between any two equivalent reduced stabilizer graphs. Since the conversion of an 
arbitrary stabilizer graph to reduced form never increases the number of hollow 
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nodes and the rules E(i) and E(ii) preserve the number of hollow nodes, we conclude 
that the reduced stabilizer graphs for a stabilizer state are those that have the least 
number of hollow nodes. 
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Proof of the Pauli-product 
measurement rule 



In this Appendix we derive the measurement rule given in Subsec. 2.3.7 for trans- 
forming a stabihzer graph under a measurement M, where M is any tensor product 
of I and Z Pauh operators. The proof proceeds in three stages. In the first, we 
determine the effect of M , considered as a Clifford unitary operation, on a stabilizer 
state \il)). This enables us, in the second stage, to find the action of the measure- 
ment projector P^ = [/ + (— l)™'M]/2 on |^) and thus to determine whether the 
measurement is certain or random. The post-measurement state is then found via 
a simple circuit identity in the last stage. Notice that, in determining the effect of 
M on I?/)), we must retain the overall phase, since, in the second stage, the overall 
phase becomes a relative phase in the superposition Wip) + {—l)™'M\%p)]/2. 
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Figure CI: Circuit identities used to determine the effect of a Z operator on a 
stabilizer state. Identity (a) follows from the fact that Z and controUed-sign gates 
commute and from the identity HZH = X. A similar identity holds if the node 
has a loop or a sign since Z commutes with itself and with S. The first equality 
in identity (b) is easily verified in the standard basis, and the second equality is an 
application of (a). Identity (c) follows from the equalities shown; the second equality 
uses the identity in (b). 



C.l Action of Z on a stabilizer state 

To begin, consider the effect of a unitary M on the n-qubit stabilizer state \ip), where 
M is known to be a tensor product consisting of only the operators / and Z. The 
action of I is trivial, so we can focus on determining the action of Z. As illustrated 
in Fig. C.l (a), applying Z to a solid node is equivalent to applying an X operator to 
the same node prior to all other Clifford gates in the circuit. Similarly, the action of 
Z on a hollow node can be reexpressed using the circuit identity in Fig. C.l(c). This 
identity shows that an identical state is obtained by adding a leading X operator to 
each neighbor of the hollow node and, if the hollow node has a sign, introducing an 
overall phase of —1. For properly simplified measurements, hollow measured nodes 
are only neighbored by solid measured nodes, so only members of A^s? the measured 
solid nodes, collect leading X operators. The number of X operators collected by 
each member of TVIs is 1 plus the number of neighbors it has in the set Ain- Since 
X"^ = I, the net result is that X operators are added only to members of A^se, the 
set of solid measured nodes with an even number of hollow measured neighbors. 
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Summarizing, the stabilizer state M\iIj) can be obtained from the state |0)®" by 
first applying an X to each member of TWse and then applying the Clifford gates 
needed to obtain {ip) from the initial state |0)®". In addition, an overall phase of 
(—1)^ must be applied where b = {Ain fl Z\ is the number of measured hollow nodes 
with a sign. That is. 



M|^) = Mf/|0)®" = (-l)^f/( Y[ XJ|0)®", 

VjeXsE / 



(C.l) 



where U is shorthand for the sequence of Clifford gates for preparing {ip) = f/|0)®"' 
as indicated by the stabilizer graph. 



C.2 Certain and random outcomes 

The second stage of our proof applies this result to find the action of Pm on {ip) 
and the probability of obtaining measurement outcome (—1)™, which is given by 
{ip\Pm\'ip) ■ We have immediately that 

\ jGXsE / 



which gives 



\ Vie-MsE / / 



l^™|^) = oU + (-l) (0| I 11 ^.■llO)^"| • (C.3) 

Vje-MsE 

Since (0|X|0) = 0, this means that measurements are classified into two types: if 
A^SE = 0, the outcome probabilities are 1, for m =2 b, and 0, for m 7^2 b, but if 
■M.SE ^ ^ , ITT- has a 50% chance of being either or 1. 
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(a) |o>4if} 



^W}*y^{BW- |0)^.|. cp [Wy \0)^H}*jj{Wy 



n 



(b) |o>THn-M-rff 



\o)-{hU^{W^ 



(c) 






i|0>^ 



|o)- 












|0>^i]- 




^Wh 



Z"'- 



-^^ 




Figure C.2: Circuit identities used to determine the post-measurement state. Iden- 
tities (a) and (b) utilize basic circuit identities for pushing a controUed-X gate past 
a controlled-Z gate. For the identity in (c), the first three qubits represent the ele- 
ments of A^sE, with the chosen qubit p placed on top. The dashed box delineates 
the Clifford operations that put the qubits in A^se into an appropriate cat state, 
i.e., those operations in Eq. (C.5) that are applied before [/, whereas the gates after 
the dashed box are the relevant portion of the gates in f/, i.e., the gates that create 
the original stabilizer state. For illustration, V contains no S* or Z gates except that 
the sign of the chosen node is indicated by Z'^ . The first equality in (c) is trivial. 
The second uses the identity in Fig. C.l(b) to push Z'"+^ to the right, turning it into 
j^m+h g^j-^^ depositing Z"^^^ on each neighbor of p. The next equality eliminates the 
controlled-Z gates that initially connected p to other nodes using the identity in (a) 
for connections to nodes in A^se and using (b) for connections to nodes outside of 
A^SE- The final simple equality returns the circuit to graph form and produces a 
possible overall phase which is ignored. 
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C.3 Post-measurement state 

When the measurement outcome is (—1)™", the post-measurement state is {ip') = 
Prn\4') / ^/wl^^riW} ■ When A^sE = 0, the outcome m =2 6 occurs with certainty, 
and the post-measurement state is the same as the original stabihzer state {ip)- 

When A^sE 7^ 0, the situation is more comphcated. In this case (iplPmlip) = 1/2, 
so 



^^ V jeMsE J 



^^ \ jeMsE 

Thus the post-measurement state is obtained using the Chfford circuit that created 
the original stabilizer state, but applied to an initial state that, for the qubits in 
A^sE, is changed to a cat state, i.e., an equal superposition of all Os and all Is, with 
the sign between the two terms in the superposition given by (—1)™+^. To construct 
a graph for the post-measurement state, we introduce a standard quantum circuit for 
making the cat state from lO)*^" and use circuit identities to put the overall quantum 
circuit into graph form. Thus we write 

\^')=u( n ^')( n %^z;^^' n H,\or, (c.s) 

\ieMsE\{p} / VfceA4sE\{p} / jeXsE 

where p denotes the chosen node. Figure C.2 translates these Clifford operations into 
circuit notation and develops the identities needed to determine the post-measurement 

state \ip'). 

Steps 1-3 in Sec. 2.3.7 follow directly from Fig. C.2(c), while step 4 follows from 
applying an S gate to both sides of Fig. C.2(c) and using T3. Step 1 describes the 
edge complementation that occurs in the third equality in the figure. Step 2 correctly 
describes the node signs indicated in the final circuit. Finally, step 3 expresses the 
fact that neighbors of the chosen node are now those nodes in A^se, and the chosen 
node is hollow. 
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Proof of Theorem 1 



Before we can proceed with the proof of Theorem 1 there is one definition to get 
out of the way. The correlation range of n random variables -R(ai), . . . , -R(a„) is the 
range of possible values for the expectation value of their product. This interval is 
denoted as C(ai, . . . , an) so that R{ai) ■ ■ ■ R{an) = R{p) -^ p E C(ai, . . . , a„). The 
idea of correlation ranges is crucial for our proof. 

We prove Theorem 1 by induction on the number of random variables. It is 
convenient that the case n = 2 is already covered. Thus we assume that the theorem 
is true for n variables and we wish to find all the possible p such that 

R{ai) ■ ■ ■ R{an)R{an+i) = R{p) 

^ Vj R{p)R{a,) = n,^, R{ak) (D-1) 

^ C{p, ttj) n C(ai, . . . , ttj-i, ttj+i, . . . , ttn+l) 7^ 0. 

The inductive hypothesis now lets us find both these correlation ranges to see if they 
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intersect. We can immediately write down the second correlation range. 

C(ai, . . . , flj-l, flj+l, . . . , Ctn+l) = 



(D.2) 



Since the sequence in Eq. (D.l) is valid for any j, we go ahead and find the correlation 
ranges for some j 7^ n + 1. 



Case 1: p < 0. 



. There are now three conditions 



In this case, C{p, aj) = |p + a^l — 1, 1 — (a^ —p) 
which, taken together, are equivalent to the correlation ranges overlapping. The first 
is that the correlation ranges be increasing, i.e., that they are actually ranges, which 
is really just a technical requirement equivalent to p G [—1, 1]. The second is that 
the lowest value in C{p, aj) be below the highest value in C{{ak}k^j), 

n 

\p + aj\- I < (n- 1) -^ak + an+1. (D.3) 

The third is that the highest value in C{p, aj) be above the lowest value in C{{ak}kj^j), 

l-aj+p>^ak-{n-l). (D.4) 

We want to show that these are equivalent to p G [—1, 1] and 

n 

^Ofc -n <p< r2-^afc + a„+i. (D.5) 

k fe=l 

[^] First we show that C{p, aj) being a range gives us that p G [—1, 1]. For, 
—p — flj — 1 < 1 — (flj — p) or 



ij + p — 1 < 1 — {aj — p) 



(D.6) 



^ p> -1. 
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Furthermore, p < 1 since it is negative. Now we show that Eqs. (D.3) and (D.4) 
imply Eq. (D.5). From Eq. (D.4) we get 

1 - aj + p > 22 o-k - {n - I) <^ p > 2J a-k - n, (D.7) 

k^j k 

which gives us the lower bound in Eq. (D.5). The upper bound in Eq. (D.5) is 
positive by the fact that a^ G [0, 1]. Therefore, since p is negative in this case, we 
get the entirety of Eq. (D.5). 

[<^] First, it is easy to see that p G [—1,1] implies that the correlation ranges are 
increasing from the backwards arrow in Eq. (D.6). Now, as we have seen, Eq. (D.4) 
is equivalent to p > J2k ^k — n. To show Eq. (D.3), we have to consider two cases. 

li p + ttj > 0, consider the upper bound in Eq. (D.5). 

n n 

p <n-^ak + an+i =^ p + ttj - 1 < (n - 1) -^a^ + a„+i. (D.8) 

fc=l k^j 

If p + ttj < 0, consider the lower bound. 

P>Y.k(^k-n>Y,f^ak- 2(a„+i + aj) - n 



(D.9) 



^ -(p + aj) -1 <{n-l) - YTk^j ^fc + "ri+i- 
And we have the result. 



Case 2: p > 0. 

Now, C{p, aj) = 
to 



p + aj — 1,1 — \aj —p\ 



and the ranges overlapping is equivalent 



p + ttj - 1 < (n - 1) - ^ ttfc + a„+i (D.IO) 

ky^j 



l-|aj-p| >5^a,.-(n-l). (D.ll) 

kj^j 
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We want to show that these equations, along with the correlation range being a valid 
range, is equivalent to Eq. (D.5) when supplemented with p G [—1,1]. 

[^] As before, we begin by showing that C{p, aj) being a range gives us p G [—1,1]. 
For, 

p + aj — 1 < 1 — {aj — p) or 

p + aj-l<l-{p-aj) (D.12) 

^ p <l. 

Again, p being non-negative takes care of the lower bound. There are now two cases 
to consider. 

If p < a„+i, then from Eq. (D.ll) we get 

1 — aj + p> \^ o-k — (n — 1) ^ p > y^ Ofc — n. (D.13) 

kjtj k 

But then, 

(D.14) 

<n- YJk=l (^k-P + fln+l- 

This means that p < n — Y12=i ^k + dn+i- 

If p > a„+i, then we want to consider Eq. (D.IO). 

n n 

p + aj — 1 < {n — 1) — y^ ttfc + a„+i -^ p < n — N^ a^ + a„+i. (D.15) 

k=/=j k=l 

Now the other inequality follows from this. 

n 

> y^ Ofc — a„+i + p — n > \J cik— n— p^p> y^ ajt — n. (D.16) 

k=\ k k 
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[<^] The correlation ranges that yield the overlapping conditions in Eqs. (D.IO) 
and (D.ll) are valid ranges by Eq. (D.12). Also, we have already seen that Eq. (D.IO) 
is equivalent to p < n — J2t=i ^k + o-n+i- So, to show Eq. (D.ll), consider the two 
cases. 

If p < flj, consider the lower bound in Eq. (D.5). 
p > 2jctfc -n^l+p-aj> 2jctA: - {n- 1). (D.17) 

k ky^j 

li p > ttj, then consider the upper bound. 

p <n- Yl^j ak - aj + a„+i 

< ^ - YlljLj «fc - ctn+i + a-j (D- 18) 

^l-p + aj> J2k^j Ofc - (n - 1). 
Thus we have proven the result in this case as well. 
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